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INTTRCOUCTICN 


This  report  describes  the  results  achieved  by  the  Regis  College  Research 
Center  during  a 26-nionth  contract  with  the  Air  Force  Geophysics  Laboratory 
(APGL) . The  scientific  disciplines  under  investigation  were  varied  and  all  of 
than  paralleled  some  of  the  efforts  undertaken  by  the  Meteorology  Division, 
APGL.  llTe  Research  Center  werked  on  three  general  areas:  atmospheric  analysis, 
mesoscale  forecasting,  and  satellite  meteorology. 

In  atmospheric  analysis  we  performed  a number  of  modeling  ejqjeriments 
to  explore  flux-profile  relationships  in  the  atmospheric  surface  layer.  In 
mesoscale  forecasting  we  concentrated  on  the  development  of  autonated  proce- 
dures for  specifying  the  type  and  intensity  of  surface  weather  phencmena  and 
the  development  of  indirect  methods  for  specifying  and  predicting  landing 
zone  slant  visual  range.  In  satellite  studies,  we  used  data  acquired  by  the 
Man-ocrnputer  Interactive  Data  Access  System  (McIDAS)  to  perform  analyses  in 
the  following  areas:  the  physical  properties  of  clouds,  the  objective  extrac- 
tion of  satellite  sensed  vertical  soundings,  the  use  of  imagery  data  for 
invert Lgating  automated  short  range  forecasting  techniques,  and  the  objective 
specification  of  cloud  type  fron  the  Defense  fteteorological  Satellite  Program 
(CMSP)  very  high  resolution  satellite  imagery  by  advanced  statistical  methods. 
This  last  effort  in  CMSP  imagery  was  conducted  in  collaboration  with  the  Air 
Force  Global  Weather  Central  (APOC:)  and  was  undertaken  because  of  AFQCs 
need  for  assistance  in  this  area. 


1.0  ATMOSPHERIC  ANALYSIS 

1.1  Error3  of  Est±nates  in  the  Atmospheric  Surface  Laver 

1.1.1  Introduction 

Many,  if  not  all,  arpirical  formulas  in  the  physical  sciences  are 
determined  by  sene  optimization  procedure  applied  to  sets  of  measurements.  A 
frequently  used  procedure,  for  exairple,  is  the  method  of  least  squares  in  \diich 
the  best  formula  is  the  one  rendering  the  sum  of  squares  of  the  differences 
between  estimated  and  measured  values  a minimum  for  a given  set  of  measurements. 
The  best  formula  obtained  in  this  manner  is  thus  subject  to  variations  due  to 
changes  in  the  set  of  measurements  anployed  in  the  c^imization  and  to  changes 
in  the  variable  optimized.  It  is  therefore  no  mere  idle  curiosity  either  to 
ask  how  good  a best-fit  arpirical  formula  is  or  to  examine  hew  sensitive  such 
a formula  is  to  differences  in  the  optimization  procedure.  The  answers  to  these 
questions  are  important  not  only  in  determining  the  accuracy  of  the  formulas 
but  also  in  assessing  the  utility  of  the  formulas. 

Businger,  wyngaard,  Izuni  and  Bradley  (1971)  proposed  the  following  set 
of  fon.^ulas  for  the  relationship  between  the  vertical  gradients  of  the  mean 
wind  and  potential  temperature  and  the  eddy  flirces  within  the  atmospheric 
surface  layer: 

(1)  l^)se  stratification 

di  = 35  ^ (1  + 15  j;)  (1.1) 

= .74  0*  z-Ml  + 9 (1.2) 

dz  L 


> 1 
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(2)  stable  stratification 


dz 


u* 

.35 


z"^(l  + 4.7  J) 


(1.3) 


^ = e*z”^  (.74  + 4.7  ^ ) (1.4) 

dz  ^ 


in  which  V is  wind  speed,  0 potential  barperature,  and  z the  height  above 
the  ground  at  vAiich  V and  0 are  referred.  Three  ndcrcmeteorological 
parameters  u*,  fi*  and  L are  defined,  respectively, by 


u* 


= ( Iq)^  = ( -uw  )^ 


(1.5) 


0*  = 


Ho 


/ .35  u* 


wS" 

.35  u* 


(1.6) 


3 - 

L = - li*— 

.35  q m 


(1.7) 


vrfiere  u*  is  the  friction  velocity,  0*  the  scaling  tenperature  and  L is 
the  Ntinin-ObuJthDV  length.  9 is  the  mean  potential  tarperature  of  the  atmospheric 
surface  layer. 

S#  and  0w  are  the  covariance  between  the  longitudinal  and  vertical 
oorponents  of  wind  and  the  ooveuriance  between  the  potential  tenperature  and 
the  vertical  oorponent  of  wind,  respectively. 

These  formulas  were  derived  in  accordance  with  the  principle  of  least 
squares  using  the  results  of  measurements  taken  in  the  Kansas  1968  Field  Program 
of  the  Air  Force  Cambridge  Research  Laboratories.  A carplete  report  on  the 
measuranents  taken  and  the  data  oollected  in  the  field  program  has  been  prepared 
by  Izuni  (1971)  in  the  Kansas  1968  Field  Program  E)ata  Report  (hereafter  referred 
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to  as  the  Report) . The  veracity  and  significance  of  the  formulas  rest  both  cai 
the  careful  measuranents  and  the  meticulous  analyses  of  the  data,  and  are 
highly  appreciated  by  the  scientific  comtunity  (Panofsky  (1974) ) . 

Ironically,  it  is  this  recognition  of  excellence  that  has  led  us  to  raise 
the  question  concerning  the  utility  of  these  formulas.  For,  in  almost  any 
other  perceivable  situation,  in  vAiich  some  applications  of  the  formulas  are 
desired,  we  would  not  expect  to  have  the  sophistication  of  instnmentation  and 
the  accuracy  of  measuranents  as  those  that  were  present  in  the  Kansas  1968  Field 
Program.  It  is  more  likely  that  the  formulcis  would  be  utilized  in  conjunction 
with  data  that  are  obtained  in  a less  refined  manner.  A few  inevitable  questions 
then  are:  Can  these  formulas  still  be  used  with  the  same  degree  of  confidence? 
What  reservations  should  be  kept  in  applying  these  formulas? 

The  work  reported  belcw  is  our  attarpt  to  answer  these  questions  without 
the  benefit  of  another  field  program.  By  depriving  ourselves  of  seme  of  the 
measurements  ctotained  in  the  field  program  and  enployed  in  the  derivation  of 
Bq.  (1) , we  ccaistructed  a hypothetical  set  of  measurements  that  might  be 
c±>tained  with  less  ccmplete  yet  practical  instinmentation.  We  then  sought 
the  best-fit  formulas,  in  those  that  would  give  the  leeist  square  errors  of 
estimate  in  a slightly  different  but  equally  sound  manner  as  that  employed  by 
Businger  et  al.  Cotparisons  were  made  among  various  choices  of  values  of 
parameters  in  these  formulas  by  calculating  the  square  errors  of  estimate  as 
defined  here  using  the  data  available  in  the  field  program  as  appeared  in  the 
Report. 

1.1.2  Experiments 

While  adhering  to  the  same  aissiitption  about  the  surface  layer  and  to  the 


same  notican  about  optimality  as  those  atopted  by  Businger  et  al.,  vge  have  chosen 
variations  to  two  major  decisicais  they  bode  in  the  course  of  arriving  at  Eq.  (1) . 

These  deviatiems  resulted  in  the  use  of  different  data  for  estimating  values  of  ^ 

the  parameters  and  different  objective  functicais  for  minimizing  the  difference 
betMeen  estimated  and  measured  values.  We  shall  elaborate  upon  these  differences 
in  this  section. 

On  the  bcisis  of  analyses  by  Haugen,  Kaimal  and  Bradley  (1971) , Businger  et 
al.  assone  that  the  heat  and  mcmentun  fluxes  ranain  constant  with  height  in  the 
atmospheric  surface  layer  and,  therefore,  parameters  u*,  0*  and  L in  Eg. (1) 
are  constant  in  each  sanple  case.  However,  in  order  to  reconcile  the  mecisure- 
ments  of  the  Reynolds  stress  by  drag  plates  at  the  ground  and  by  the  sonic 
anemcmeters  on  the  tower,  they  recognized  the  trend  of  variation  of  the  measured 
uw  with  height  between  5.66  m-  and  22.63  m-level  and,  took  the  latter  into 
account  in  estiitating  the  surface  value  of  the  stress.  The  friction  velocity 
u*  characteristic  of  the  surface  layer  was  Ccilculated  fron  this  surface  value 
of  the  stress.  On  the  other  hand,  they  thought  it  best  to  average  the  values 
of  heat  flux  measured  at  three  different  levels  (5.66,  11.32  and  22.63  m)  in 
defining  the  value  of  the  surface  layer  heat  flux.  " 

We  imagine  that  we  set  up  the  required  instrumentation  a posteriori  to  the 
assurpticxi  that  the  fluxes  remain  constant  with  height.  It  would  therefore  be  ^ 

natural,  for  reasons  of  econaiy,  to  measure  each  flux  at  cxie  level.  This  we  do  | 

at  the  lowest  level  available  in  the  Report.  Thus,  for  the  momentun  flux  we 
take  the  average  of  two  drag-plate  roeasuranents,  corrected  for  the  difference  * ^ 

between  the  plates  and  the  sonic  anancmeters.  It  is  the  same  as  that  used  by  li 

Businger  et  cil.  For  the  heat  flux,  however,  ve  take  the  value  of  wS"  measured  I 
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at  z = 5.66  m as  the  value  r^resentative  of  the  surface  layer,  in  contra- 
distinction to  that  of  Businger  et  al.  Likevd.se,  while  Businger  et  al.  ncoH 
the  value  of  0 at  z = 4m  as  the  estimate  of  0 in  the  expression  of  L,  we 
enploy  the  vailue  of  0 at  z = 2m  for  0 . 

In  minimizing  the  difference  between  estimate  and  c^aservation,  Businger  et 
al.  chose  for  the  objective  fianction  of  minimization  the  sun  of  squares  of 
differences  between  the  estimated  and  observed  vertical  gradients.  Since  the 
vertical  gradients  were  never  measured  directly  in  practice,  they  resorted  to 
an  alogarithm  to  obtain  the  so-called  observed  valxies  (Businger  et  al.  (1971)) . j 

As  shcvm  in  ^^pendix  1.1,  since  the  assixned  vertical-gradient  equations  may  j 

be  readily  integrated  with  respect  to  height  to  yield  the  difference  between  any  j 

tWD  levels  above  ground,  we  choose  for  the  objective  function  of  minimization  | 

the  weighted  sun  of  squares  of  differences  between  estimated  and  observed  I 

I 

differences  between  two  consecutive  levels.  The  weights  are  chosen  to  be 
prcportional  with  the  difference  in  logarithm  of  the  height.  Thus,  all  the 

sublayers  between  2m-  and  32m-levels  considered  in  this  study,  from  2m  to  4m,  i 

from  4m  to  &n,  from  8m  to  l&n,  and  from  16m  to  32m,  have  the  same  weight. 

1.1.3  Opti  'ization  of  the  gradient  equations  * 

We  assune  that  the  vertical-gradient  equations  are  given  by  ( 


(1)  l^)se  stratificaticMi 


(2)  Stable  stratification 


dz 


Z ^ (By  + ^yZ  ) 


(3.3) 


^ z-MBjj  + 6hZ  ) (3.4) 

in  which  (kv,  Yy,  ny) » (^H'  '^H'  "h>  ^^P®®  stratification  and  (By,  6y) 

(Bjj,  6^)  under  stable  stratification  are  to  be  determined  such  that  quantities 

such  cis 


N 5 
E E 
n=l  ]c=2 


K - Vi> 


Wk  - Vk-ll] 


2 

n 


vrfiere  symbol  - is  used  to  denote  the  estimate,  are  minimized.  Here,  integer 
index  k derxjtes  the  level  in  the  vertical,  with  k=l  referring  to  the  lowest, 
z = 2m  and  k=5  to  the  highest,  z = 32m.  Each  value  of  integer  index  n identi- 
fies one  saitple  case  vAiicdi  is  a 15  minute  period  during  vAiich  a set  of  concurrent 
measuranents  on  V,  ©,  iw  and  is  available  (see  Tables  1.1  and  1.2  of  the  Report) . 
Vfe  shall  call  these  the  universal  ccsistants  of  the  flux-profile  relationships. 

Excluding  runs  46  axKi  47  (as  called  in  the  Ri^rt)  in  view  of  missing  data 
of  uw  and  we  at  z = 22.63  m,  we  find  altogether  68  cases  with  lapse  stratifica- 
tion and  42  cetses  with  stable  stratificaticai.  Each  grocp  has  been  treated  as 
one  sanple  groc?*  characteristic  of  the  respective  stratificatica'.,  after  the 
initial  trial  of  dividing  the  groi^js  into  subgrxxps  of  ^3fxraximately  twaity 
cases  fadled  to  reveal  any  sign  of  (consistency  in  the  optinal  values  of  the 
universed.  constants  among  different  subgroups. 

The  cptimed.  sets  of  the  universeil  constants  are  sought  within  the  following 


(Confines: 


1 

(1)  ky  and  1^  are  positive  real  nimbers  with  two  decimal  places, 

(2)  Yy  and  Yjj  are  real  integers, 

(3)  n^  and  njj  are  real  integers, 

(4)  6y,  gjj  and  cire  real  nutnbers  with  one  decimal  place. 

It  is  readily  recognized  that  all  the  values  proposed  by  Businger  et  al.  fall 
within  the  above  restrictions.  That  these  restrictions  have  a sufficiently  fine 
resolution  for  the  given  set  of  measurements  may  be  seen  from  the  vad-ues  of  root- 

mean-square  errors  of  estimate  tabulated  in  Tables  1.1a,  1.1b  ani  1.1c.  i 

1 

Table  1.1  is  presented  to  demonstrate  the  structure  of  topography  of  the  j 

root-mean-square  errors  of  estimate.  Each  of  the  tables  gives  a three-dimen-  j 

I 

sionail  cross-section  on  a hyperplane  defined  by  n^  = const.  The  entry  in  the 

Icist  oolvinn  under  (2.32) , vAiere  the  overbar  signifies  averatge,  gives  the  value  | 

of  the  c*)jective  function  of  optimization  at  a given  set  of  (k^,  Yy)  v^iich 

defines  the  coordinates  of  tiie  valley  of  the  tc^ogr^iiy.  Each  entry  in  other  ^ 

oolunns  is  the  root-mean-square  error  of  estimate  of  the  wind-speed  difference 
between  the  two  levels  bounding  the  sublayer.  The  anallest  value  in  each  oolixnn 
is  underlined  for  easy  recognition.  ; 

It  appears  fron  these  tables  that  there  are  a nuonber  of  equ^Qly  acc^Jtable 
choices  for  the  cptimed  set  of  universal  constants.  By  using  the  values  of 

mean  square  error  of  estimate,  hcvever,  it  is  possible  to  discriminate  further  j 

4 / ICy  • 34 1 Yv  ~ • I 

f 

The  finer  structure  cf  topograEhy  of  the  average  square  of  error  of  p 

estimate  in  the  vicinity  of  the  optimal  set  is  given  in  Table  1.2.  A cotpariscni  ^ 
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of  the  vcilues  of  the  avereige  squcire  errors  between  this  set  and  that  of 
Businger  et  ad.,  shows  that  the  difference  is  insignificant.  Similar  procedures 
au?e  follcwed  in  the  determination  of  the  optimad  sets  for  other  categories.  Ttie 
results  for  la^ise  stratification  are  suimarized  in  Table  1.3a  and  those  for  stable 
stratification  in  Table  1.4a.  Acoonpanying  each  set  is  the  reoGKmended  set  of  the 
universal  constants  vhich  is  obtained  vhen  the  optimadity  of  the  cA>jective 
function  in  the  optimal  set  is  given  vp  for  uniformity  of  values  of  constants. 

The  values  of  the  root-mean-square  errors  of  estimate  are  included  in  the  table 
to  shew  the  effect  of  such  a trade-off. 

Tables  1.5  and  1.6  are  pr^ared  to  exhibit  finer  details  of  these  selected  sets 
of  the  universad  oonstamts  as  well  as  of  the  sets  proposed  by  Businger  et  ad. 
Superscripts  r,  o and  b are  entered  in  ^^aropriate  places  to  designate  the  rows 
referring  to  the  recanoended  sets,  the  optimad  sets,  and  the  sets  of  Businger 
et  al,  re^jectively.  It  is  cleau:  fron  these  tables  that,  meaisured  in  terms  of 
the  objective  function  enployed  in  this  study,  the  sets  of  Businger  et  ad.  differ 
little  fron  the  other  two  in  both  wind  speed  and  potentiad  temperature  under 
laipse  stratification.  Hiis,  hcMever,  is  not  the  case  under  stable  stratificaticai. 

In  both  wind  speed  amd  potential  temperature  they  are  associated  with  considerably 
larger  errors  of  estimate,  paurticularly  in  the  vpper  two  sublayers. 

1.1.4  Uses  of  the  gradient  equations 

There  aure  two  practiced  uses  to  vAiich  the  verticad  gradient  equations  can 
be  put.  One  is  to  estimate  wind  speed  and/or  potentiad  tenperature  within  the 
surface  layer  at  a level  away  fron  the  level  of  measurement.  ‘iTiis  may  be 
akxxitplished  by  integrating  the  gradient  equations  if  there  are  measurements  or 
estimates  of  the  micronetearological  paurameters  together  with  the  meaisurenent  of 


» 


i 

I 


tv, 

f 
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wind  and/or  taipearature  at  one  level,  llje  other  is  to  estimate  the  values  of 
eddy  fluxes  from  multi-level  measurements  of  wind  and  tenperature.  Ihis  is  done 
by  inverting  the  integrals  of  the  gradient  equations,  "nie  application  for  such 
purposes  of  the  gradient  equations  c^Ttained  in  the  previous  section  are  di  .ctmigsgH 
in  the  following. 

1.1. 4.1  Estimation  of  wind  and  tenperature  with  the  aid  of  micrcroeteorological 
measurements. 

A cotparison  of  errors  of  estimate,  vAiich  are  c±>tained  with  the  use  of  the 
data  in  the  R^ort,  is  made  using  three  different  methods  of  estimating  wind  and 
terperature  at  various  levels  fran  the  measurements  taken  at  z = 2m. 

Ihe  micranetearological  method  makes  use  of  the  Vcilue  of  u*  and  0* 
calculated  fran  the  measurements  of  the  mcmentun  flux  at  the  ground  and  of  the 
heat  flux  at  z = 5.66  m.  Hie  value  at  level  z^,  F(Z]^) , is  estimated  by 

F<V  - F(zi)  + ^ f dz  (4.1) 

vtere  zj  = 2m  and  ^ is  the  verticeil  gradient  defined  by  caie  of  Eqs  (3.1)  'v 
dz 

(3.4)  with  I je  oorrespcffiding  optimal  set  of  the  universal  constants. 

The  method  of  linear  regressicai  estimates  F(zj^)  by 

F(zjt)  = + B,5F(zi)  (4.2) 

vrfiere  are  determined  by  linear  regression.  The  method  of  polynomial- 

in-logarithn  assunes  that 

F(z,^)  = F(zi)  + P4(ln  ^ ) (4.3) 


J 


' I 
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in  \*u.ch  (x)  is  the  fourth-order  polynanied  in  x.  "Hie  veilues  of  edl  oonstcints 
in  the  polynanieLL  are  determined  in  acxxsrdance  vdth  the  least-squares  principle, 
with  the  aid  of  the  same  set  of  data,  eis  were  eind  in  tlie  method  of  linear 
regression,  that  is  used  in  the  micrometeorological  method. 

Ihe  results  are  surmarized  in  Table  1.7.  It  shews  that,  at  leeist,  with  the 
set  of  measurennents  avciilable  there  is  no  evidence  suggesting  that  the  micro- 
meteorologiccil  method  is  si?)erior  to  the  other  two,  cis  far  as  the  accuracy  of 
estimation  is  oonoemed.  Cn  the  contrary,  the  errors  in  stable  stratification 
cure  considerably  larger  than  those  of  the  other  two  methods.  Uiis  may  be  due 
to  the  fact  that  the  thickness  of  the  surface  layer  cam  be  much  less  than  32 
meters  under  stable  stratification  and,  consequently,  the  flux-profile 
relationships  used  cu:e  no  longer  valid  in  the  layer  above  z = 8m. 

1.1. 4. 2 Estimation  of  eddy  fluxes  with  the  aiid  of  multi-level  measurenents 
of  wind  eind  tenperature 

With  the  use  of  definiticxis  of  friction  velocity  u*  and  scaling  tanpera- 
ture  0*  given,  respectively,  by 

22  = u*2  (4.4) 

p 

and 

— = -kvU*0*  (4.5) 

pc^ 

where  p is  the  eiir  density  and  Cp  the  specific  heat  of  air  at  constant 
pressure,  the  problem  of  estimating  the  momentun  flux  Tq  and  the  heat  flux 
in  the  surfeKse  layer  is  reduced  to  that  of  estimating  u*  and  0** 
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In  designing  the  schemes  of  estimating  u*  and  0*  from  multi-level 
measurements  of  wind  and  temperature  that  are  described  below,  we  have  followed 
the  sequential  procedure  in  vAiich  the  parameters  are  estimated  one  at  each  st^ 
uniquely  and  the  value  of  a peurameter  estimated  at  a later  st^  depends  on  the 
values  of  those  estimated  at  earlier  steps. 

(a)  Scheme  I - This  method  assixnes  no  knowledge  of  the  roughness  length 
Zq  and  requires  measurements  of  wind  at  three  different  levels  and  of  tempera- 
ture at  two  different  levels.  Let  these  be  denoted  by  = V(zi) , i = 1,  2,  3 
and  Qj  = 0(zj),  j = 1,2,  respectively.  Hie  estimation  procedure  runs  as  follows: 
St^  1:  Obtcdn  ratio  defined  by 


Rl  = ^3  - V2 
V2  - Vi 


St^  2:  Estimate  the  Monin-Obukhov  length  L from  with  identity 

^ lv(zi,  Z2;  L) 

vhere  Z^(zi.  Z2;  L)is  the  integral  of  the  gradient  for  wind  with  respect  to 
height  frcm  level  zj^  to  level  Z2  with  the  value  of  the  Monin-Obukhov  length  L. 
Step  3:  Estimate  u*  frcm  V]^,  V2,  V3  and  L in  accordance  with  the  least-squares 
principle.  It  may  be  readily  shown  that  u*  is  given  by 

; _ Z27  L)  (V2  - Vi)  + Iv(22/  Z3?  L)  (V3  - V2) 

u*  ~ 

Iy  ^^2.^  ^2^  ^v  ^^2^  ^3^ 


f ^ 

I • 


t 

r'' 


L 
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Sti^  4;  Estimate  0*  fron  9]^  and  ©2  and  L with  identity 


= ^2  ~ Ql 

l0(Zj/  Z2;  fi) 

where  I0  is  the  integral  of  the  vertical  gradient  for  potential  temperature 
between  Zj  and  z^. 

St^  5:  Estimate  Zq  fron  Vj  and  u*  with  identity 

k^V  j — ( Zq  t z j ; L) 

l*ider  stable  stratification,  since  both  integrals  and  I0  are  linear 
functions  of  L,  the  equation  in  step  2 may  be  rewritten  as  an  e:q>licit  equation 
for  L are!  the  values  of  and  I0  can  be  readily  determined  with  the  knowledge 
of  L.  The  only  transcendental  relation  arises  in  Step  5. 

Under  l^se  stratification,  however,  either  tabular  or  grc4>hical  repre- 
sentations of  Iv  as  a function  of  L with  Zj  and  Z2  as  parameters  are  required 
to  obtain  L.  We  have  found  the  tabular  r^resentation  with  variable  ’jicrement 
in  L seirves  the  purpose  quite  adequately. 

It  is  readily  shown  that  the  proposed  vertical  gradient  of  wind  restricts 
the  value  of  R]^  within  certain  ranges  that  depend  on  the  levels  of  meacure- 
ments,  regardless  of  vailues  of  the  universal  constants.  For  example,  with  the 
choice  of  wind  measurements  at  z = 2,  4 and  8m,  it  is  expected  that 
.8409  < Rj  = 1 under  l^se  stratification 

and 

1 < Rj  < 2 under  stable  stratification 
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Consequently,  whenever  the  ctoserved  winds  give  rise  to  values  of  outside 
the  prescrihed  ranges,  the  scheme  fails  to  yield  logically  satisfactory  estimates. 

The  conflicts  between  estimate  and  c±»servatiOTi  may  be  conveniently  divided 
into  two  kinds.  A conflict  of  the  first  kind  occurs  vhen  the  ctoserved  vcdue  of 
Rl  is  either  less  than  the  lower  limit  under  lapse  stratificaticn  or  greater 
than  the  v?3per  limit  under  stable  stratification.  In  either  case  the  nearest 
possible  estimate  calls  for  L = 0 v^iich  invalidates  the  proposed  gradient 
equation.  On  the  other  hand,  a conflict  of  the  second  kind  cAjtciins  vhen  the 
observed  value  of  is  either  greater  than  the  L^jper  limit  under  lapse  strati- 
fication or  less  than  the  lower  limit  under  stable  stratification.  In  either 
case  the  closest  estimate  calls  for  L ->  »>,  a situaticai  vrfiich  is  likely  to  cause 
inconsistencies  among  u*,  6*  and  L. 

When  we  applied  this  schane  to  the  measiaranents  in  the  lowest  levels  of 
the  R^xart,  we  found  that  29  and  7 out  of  the  68  cases  under  lapse  stratifi- 
cation had  conflicts  of  the  first  and  second  kinds,  respectively.  The  dis- 
covery v«s  so  discouraging  as  to  make  us  look  for  Scheme  II  v^ch  is  described 
below. 

(b)  Scheme  1*  - This  method  cissvines  that  the  roughness  length  Zq  is  known  and 
requires  measurements  of  both  wind  and  tenperature  at  two  levels. 

Step  1:  Obtain  ratio  R2  defined  by 


Step  2:  Estimate  L from  R2  with 

R-  = ^2> 
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step  3:  Estimate  u*  according  to 


Iv(2o,  zi;  L)Vi  + Iv(Zo,  Z2;  1)^2 

1^2  Zj^;  L)  + Iv^(Zq,  Z2;  L) 

Step  4:  Estimate  0*  according  to 
9^=  ^2  - Qi 

As  in  Scheme  I,  the  proposed  gradient  equations  of  wind  limit  the  value  of 
R2  within  certain  ranges.  With  the  choice  of  z = 2 and  4m  and  Zq  = 2aa  for 
instance,  it  may  be  shewn  that 

1.0736  < R2  ^ 1.1505  under  l^jse  stratification 

1.1505  < R2  < 2.0101  under  stable  stratification 

Similcu:  conflicts  as  those  of  may  be  shewn  bo  exist  between  estimated  and 
observed  vcilues  of  R2. 

When  we  applied  this  scheme  to  the  data  of  the  R^ort  in  the  lowest  levels, 
we  found  4 conflicts  of  the  second  kind  in  the  68  ceises  of  lapse  stratification 
and  10  conflicts  of  the  second  kind  in  the  42  cases  of  stable  stratification. 

There  was  not  a Ccise  of  conflict  of  the  first  kind  under  either  stratification. 

In  the  results  presented  belcw,  we  opted  to  ignore  the  inconsistencies  associated 
with  the  conflict  of  the  second  kind  and  obtained  both  u*  and  0*  with  estimate 
L -►  ® . This  cissured  us  of  finite  values  of  estimated  mcmentun  and  heat  fluxes. 

Figures  1.1  - 1.4  present  the  scatter  diagrams  of  u*  and  0*  in  which  the  values 


20 


estimated  with  Schane  II  are  plotted  in  ordinates  against  those  calculated  fron 
the  measurotnents  of  uw  and  w0.  The  roughness  length  Zq  = 2.44an  was  mentioned 
by  Businger  et  al.  and  is  believed  to  be  representative  of  the  site  of  the 
measurements.  The  scatter  diagrams  obtained  with  two  other  values  of  Zq,  2.00 
and  5.00cm  are  included  to  shew  the  sensitivity  of  the  estimates  to  changes  in 
the  value  of  Zq.  The  influence  of  the  accuracy  of  on  that  of  the  estimates 
on  u*  and  0*  with  Schane  II  is  quite  obvious  fron  these  figures. 

1.1.4. 3 Estimation  of  wind  and  terperature  with  the  aid  of  micrciteteorological 
estimates. 

The  twD  tables  in  Table  1.8  present  the  root-mean-square  errors  of  estimate 
obtained  vAien  the  vertical  gradients  are  integrated  to  estimate  wind  speed  and 

potentiaLL  tanperature.  Each  table  lists  the  errors  obtained  with  the  set  of  j 

I 

the  universal  constants  found  optimal  or  reconnended  in  this  study  and  those  ! 

with  the  set  preposed  by  Businger  et  al.  Along  with  each  set  of  the  constants  j 

are  two  sets  of  values  of  the  micrcmeteorological  parameters:  one  measured  in  j 

the  lowest  levels  and  the  other  estimated  fron  the  wind-  and  taiperature- 

measuraixsnts  at  z = 2m  and  4m.  The  integration  of  a vertical  gradient  starts  ^ 

at  z = 4m  ariJ  the  roughness  length  Zq  = 2.44cm  is  used  throughout.  I 

We  see  in  the  results  under  lapse  stratification  what  we  have  anticipated.  i 

In  the  first  place  the  errors  with  measured  values  of  the  parameters  are  * 

analler  than  the  errors  with  estimated  values,  with  a few  exceptions  of  anall  |. 

differences.  All  categories  shew  an  increase  with  height  of  the  magnitude  of  : 

the  error.  There  is  little  difference,  as  we  found  earlier  in  Section  1.1.3,  i' 

between  the  two  sets  of  the  universal  constants,  either  with  measured  or  with 
estimated  values  of  the  paraneters. 
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The  results  under  stable  stratification,  on  the  other  hand,  sean  to  defy 
our  expectations  and  to  elude  our  rationalizati '■n  on  one  important  count.  On 


the  one  hand,  they  are  consistent  with  those  under  lapse  stratification  in 
showing  that  the  errors  increase  with  height  in  all  categories.  They  are  also 
consistent  with  the  results  found  in  Section  1.1.3  in  that  the  errors  with  the  use 
of  the  reoanmended  set  of  the  constants  are  snaller  than  the  errors  with  the 
use  of  the  set  of  Businger  et  al.  But,  they  also  shew  that,  at  all  levels  and 
in  all  sets  of  the  constants,  the  errors  produced  with  the  use  of  estimated 
values  of  the  parameters  are  significantly  atialler  than  those  resulting  when 
measured  values  of  the  parameters  are  used.  Since  these  universal  constants 
have  been  obtained  by  optimization  of  the  verticcil  gradients  or  vertical 

differences  with  the  aid  of  the  measured  values  of  the  pararreters,  such  a j 

showing  raises  a doubt  as  to  vhether  the  optimization  procedure  was  ^iprcpriate  j 

for  the  purpose  of  the  estimaticsi  on  hand.  It  may  also  be  due  to  the  simple 

» 

fact  that  there  is  a cerponent  of  considerable  magnitude  in  these  measure-  j 

ments  that  has  not  been  accounted  for  at  all  by  the  proposed  form  of  the 
vertical  gradients . 

I 

1.1.5  Coiclusions  | 

I 

We  have  carried  out  a nunber  of  exercises  using  the  data  collected  in  the 

Kansas  1968  Field  Program  Data  Report  to  explore  the  flux-profile  relationships  < 

. h 

in  the  atmospAieric  surface  layer.  ' 

The  experiment  on  optimization  (Section  1.1.3)  shews  that  the  use  of  an  . i 

equivalent  but  slightly  different  objective  function  of  optimization  and  less  ■ 

odtplete  data  for  the  micrometeorological  peirameters  result  in  a set  of  the 

universal  constants  different  from  that  proposed  by  Businger  et  al.  There  is. 
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hcvever,  no  significant  difference  in  root-mean-square  errors  between  these 
sets  under  lapse  stratification.  The  only  serious  discrepancy  arises  in  the 
case  of  potential  tenperature  under  stable  stratification. 

The  ccnparison  of  accuracy  between  the  micrcmeteorological  method  and 
purely  statistical  methods  (Section  1.1. 4.1)  fails  to  exhibit  the  advantage  of  the 
micrcmeteorological  measuranents  so  far  as  estimation  of  wind  and  tenperature 
at  other  levels  within  tie  surface  layer  is  concerned. 

In  designing  a scheme  of  estimating  micrcmeteorological  parameters  from 
wind  and  taiperature  measurements  we  have  found  the  two-fold  inportance  of  the 
roughness  length  Zq.  (Section  1.1. 4. 2).  In  the  first  place,  the  assurtption  of  its  ^ 

knowledge  allows  a schane  v^ch  requires  one  fewer  ctoservation  on  wind  speed  ! 

and  yet  is  more  tolerant  of  the  otoservations  that  cannot  be  accounted  for  by 
the  proposed  form  of  the  vertical  gradient.  Secondly,  the  accuracy  of 
estimates  o«i  friction  velocity  u*  and  scaling  taiperature  0*  is  found  to  be 

I 

sensitively  affected  by  the  accuracy  of  Zq.  ' 

We  find  it  difficult  to  understand  v*iat  we  observe  in  the  estimation  of  ! 

wind  and  taiperature  at  upper  levels  from  those  at  lower  levels  with  the  aid  of  ^ 

estimated  mJ  crrmeteorological  parameters  under  stable  stratification  (Section  ( 

1.1. 4. 3) . We  speculate  that  the  results  are  due  to  the  fact  that  the  top  of  the 
surface  layer  under  stable  stratification  is  on  many  occasions  lower  than  sane 
of  the  upper  levels,  say  l&n  or  32m.  ^ ^ 

f 

I 

I . 
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V* 


1.1 


IntegrcLl 


z~Ml  + dz 
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The  integrcil 

22  , 

(Zj,  z^j  a)  = ^ z“Ml  + az)“^/’^  dz 

where  Zj,  Z2  > 0,  a = 0 and  n a positive  integer  may  be  transformed  into 

r^2  nv^~^  , 

In  (Yi'  y2^  = yi  (yn  - 1)  ^ 

in  vrtiich  yi  = (1  + azi)^/^  , y2  = (1  + 


For  n = 1,  2,  3 and  4,  the  integral  may  be  expressed  in  terms  of  elementary 
functions. 


I (z  , z ; a)  = In  _2. 
i i ^ zi 


InLLfi 

1 + az, 


^2' 


13(2^,  Z2' 


(1  + az2)**  - 1 

^ (1  + azj)^  - 1 


(1  + az2)^  + 1 

(1  + azj)^  + 1 


(1  + az2) 1/3  - 1 {(1  t az2)"/3  + (1  + az2)i/3  t 1}^ 


(1  + azj) 


+ ^ tan- 
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TABiz;  1.1  - The  root-nean-square  errors  of  estimate 
wind  speed  differaioe,  lapse  oonditicai 
(n\/sec) 


a)  = 3 


kv 

7v 

(2,4) 

(4,8) 

(8,16) 

(16,32) 

(2,32) 

.32 

14 

.067 

.060 

.093 

.097 

.081 

.33 

11 

.064 

.058 

.089 

.093 

.078 

.34 

10 

.065 

.058 

.086 

.092 

.076 

.35 

8 

.066 

.059 

.082 

.089 

.075 

.36 

7 

.•071 

.061 

.080 

.088 

.075 

.37 

6 

.076 

.063 

.077 

.087 

.076 

.38 

5 

.082 

.067 

.075 

.086 

.078 

.39 

4 

.087 

.070 

.074 

.086 

.080 

.40 

3 

.092 

.075 

.074 

.089 

.083 

.41 

3 

.105 

.080 

.071 

.086 

.086 

.42 

2 

.108 

.084 

.074 

.093 

.091 

.43 

2 

.120 

.090 

.071 

.089 

.094 

.44 

2 

.132 

.098 

.071 

.087 

.100 

27 


I 


L 


I 


i 

TABLE  1.1  (ocxitinued)  f 

c)  n = 5 

i 

, I Layer  


7v 

(2,4) 

(4,8) 

(8,16) 

(16,32) 

2,32) 

.32 

42 

.067 

.057 

.079 

.092 

.075 

.33 

34 

.068 

.057 

.077 

.091 

.074 

.34 

27 

.069 

.058 

.075 

.091 

.074 

.35 

22 

.073 

.060 

.074 

.090 

.075 

.36 

17 

.075 

.062 

.073 

.091 

.076 

.37 

14 

.081 

.065 

.072 

.091 

.078 

.38 

11 

.085 

.069 

.072 

.092 

.080 

.39 

9 

.092 

.073 

.071 

.092 

.082 

.40 

8 

.101 

.078 

.069 

.092 

.086 

.41 

6 

.106 

.082 

.070 

.093 

.089 

.42 

5 

.114 

.087 

.070 

.094 

.093 

.43 

4 

.122 

.092 

.071 

.096 

.097 

.44 

3 

.129 

.096 

.072 

.099 

.101 

Sie  mean  square  error  of  estimate 
wind  speed  difference,  lapse  condition 

(10~®  m^/sec^) 


rV  = 4 

K 


.33 

.34 

.35 

7523 

6063 

5770 

6954 

5786 

5523 

6511 

5613 

5592 

6175 

5528 

5730 

5928 

5515 

5926 

5759 

5564 

6170 

5655 

5665 

6453 

5608 


5810 


6770 


TABLE  1.3b  - Reocninended  sets  of  the  universal  constants 
l^)se  stratification 


Wind  Speed 

K % 

^v 

RMSE 

(nv'sec) 

Terperature 

’^h 

Yh 

RMSE 

(°K) 

.34 

4 

18 

.075 

1.00 

18 

.085 

.35 

4 

16 

.076 

1.00 

16 

.086  i 

TABI£  1.4a  - Optimcil  sets  of  the  universal  ocaistants 
stable  stratification 


Wind  Speed 


K 

^v 

«v 

RMSE 

(itv/sec) 

kh 

6h 

RMSE 

(°K) 

.34 

mn 

3.0 

.200 

1.00 

1.3 

1.1 

.184 

.35 

3.0 

.199 

1.00 

1.4 

1.1 

.184 

Tenperature 


TABI£  1.4b  - Reoatmended  sets  of  the  universal  constants 
stable  stratification 


kv 

&v 

«v 

RMSE 

(nv'sec) 

kh 

eh 

^h 

RMSE 

(OR) 

.34 

3.1 

.203 

1.00 

1.0 

.192 

.35 

3.2 

.206 

1.00 

1.0 

mi 

.193 

31 


TABLE  1.6a  - The  rootMtean-square  errors  of  estimate  bn/sec) 
wind  speed  differences,  stable  stratification 


TABLE  1.6b  - The  root-mean-square  errors  of  estimate  (°K) 

potential  tenperature  differences,  stable  stratification 

= 1.00 

Layer 


K 

(6h»'5h) 

(2.4) 

(4.8) 

(8.16) 

(16.32) 

(2.32) 

( .0,  1.2)r 

.089 

.105 

.189 

.306 

.192 

.34 

(1.3,  1.1)0 

.107 

.114 

.171 

.286 

.184 

(1.0,  1.2)1^ 

.089 

.106 

.189 

.308 

.193 

.35 

(1.4,  1.1)0 

.113 

.119 

.169 

.283 

.184 

(.74,  4.7)*5 

.152 

.321 

.677 

1.447 

.818 

r 


!• 

I 

I 

I 
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TABLE  1.7  - Ttie  root-fiean-squcire  errors  of  estimate 


a)  Wind  speed  (nv'sec) 


Strati- 

fication 

Method 

2 

4 

Level  (m) 

8 

16 

32 

- 

Micrcitet 

0 

.067 

.104 

.148 

.178 

• 

Lapse 

Lin.  Reg. 

II 

.059 

.125 

.186 

.250 

Pol.  Log. 

II 

.158 

.318 

.445 

.583 

Micrcmet 

tl 

.217 

.676 

1.563 

3.475 

Stable 

Lin.  Reg. 

II 

.049 

.112 

.181 

.303 

1 

1 ^ 

Pol.  Log. 

If 

.132 

.261 

.359 

.442 

I 

1 

b)  Potential  tarperature  (°K) 


Strati- 

fication 

Method 

2 

4 

Level  (m) 

8 

16 

32 

Micrtitiet 

0 

.253 

.345 

.465 

.565 

Lapse 

Lin.  Reg. 

II 

.226 

.346 

.429 

.490 

Pol.  Log. 

II 

.210 

.311 

.395 

.448 

Micronet 

II 

.128 

.359 

.823 

1.825 

Stable 

Lin.  Reg. 

II 

.131 

.312 

.592 

1.016 

Pol.  Log. 

II 

.213 

.292 

.349 

.383 

3. 


■ 


TABLE  1.8  (continued) 


b)  Potential  tenperature  (°K) 


stratification 

universal 

constants 

micro 

parameters 

8 

Level  (m) 
16 

32 

n^  = 2 
kv  = .34 

measured 

.105 

.226 

.332 

lapse 

kh  = 1.00 

Yh  = 20 

estimated 

.143 

.218 

.289 

rV  = 2 
kv  = .35 

measured 

.195 

.280 

.421 

kh  = 1/.74 
Yh=  9 

estimated 

.253 

.426 

.572 

II 

measured 

.236 

.707 

1.713 

®h  " 

«h  = 1.2 

estimated 

.039 

.156 

.393 

stable 

kv  = -35 

measured 

1.185 

3.594 

8.526 

Bh  = .74 

6h=  4.7 

1 

1 

estimated 

.460 

1.414 

3.430 

[■ 
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1.2 

Tayi 

1.2.1  Introduction 

Contour  maps  of  forecast  errors  are  produced  to  exhibit  the  geographical 
distribution  of  the  difference  between  prediction  and  observation  in  the  course 
of  studying  the  APGW2  boundary  layer  model.  The  maps  are  used  primarily  to  help 
us  detect  any  extraordinary  feature  or  consistent  pattern  vrfiich  may  be  associa- 
ted with  any  of  the  parameters  (such  as  time,  position,  forecast  category,  model, 
etc.)  used  to  identify  the  forecasts.  They  serve  this  purpose  very  well  either 
individually  or  in  small  groups  by  permitting  easy  visual  ccnparisons . 

As  the  number  of  the  maps  increased  with  increases  in  both  the  nixtiber  of 
synoptic  cases  and  the  number  of  models,  however,  it  became  obvious  to  us  that 
seme  systematic  method  for  categorization  and  registration  had  to  be  devised 
in  order  that  we  could  retain  correctly  and  use  properly  the  information  con- 
tained in  the  individual  maps.  What  is  needed  are  measures  of  those  attributes 
in  the  m^)s  that  can  be  efficiently  oiployed  to  characterize  the  individual 
maps  and  to  discriminate  their  various  properties.  As  such,  the  task  of  neeting 
these  objectives  falls  under  the  category  of  pattern  recognition. 

To  put  the  vrork  reported  here  in  proper  perspective  it  must  be  pointed  out 
that  in  designing  the  intervals  and  symbols  of  the  contour  maps  we  did  not 
anticipate  the  undertaking  of  this  study  and,  cis  a result,  we  did  not  have  in 
the  contour  maps  all  the  qualities  that  are  desired  for  the  purpose  of  this 
study.  Vfe  have  sinply  taken  vrfiatever  were  available  to  us  and  tried  to  see  how 
v^ll  MS  could  do  with  them. 

For  the  purpose  of  understanding  the  rationale  and  technique  of  the  methods 
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developed,  it  is  not  necessary  to  know  v^t  constitutes  each  of  the  models  or 
vrfiat  distincticxi  there  are  among  different  levels  mentioned  in  the  presentation 
below.  Nevertheless,  for  the  purpose  of  balance  and  reference,  we  present  a 
brief  description  of  the  models  considered  and  the  nature  of  the  contour  maps 
analyzed. 

Six  prediction  models  of  the  planetary  boundary  layer  are  considered.  They 

are: 

(1)  Medel  "O":  the  APGWC  operational  model, 

(2)  Ntodel  "6":  the  APGWC  model  adopted  for  use  at  the  AFGL  ooiputing  facilities, 

(3)  Model  "C" : fiodel  "6"  with  a new  method  of  estimating  eddy  ^ffusivity,  j 

(4)  Model  "H":  Model  "C"  with  virtual  tarperature  used  in  the  hydrostatic  • j 

relation, 

(5)  Model  "V":  Model  "C"  with  a new  vertical  differencing,  and 

(6)  Model  "R":  the  persistence  forecast  in  which  the  value  observed  24  hours  ago 

f 

I 

is  taken  as  the  predicted  vcilue. 

i 

All  models  predict  the  valiae  of  ciir  tenperature  at  each  of  the  27x29  <^id  | 

points  on  each  of  the  8 levels  in  the  dctnain  of  predicticm.  The  forecast  error  | 

at  each  grid  point  is  defined  to  be  the  difference  between  the  predicted  and  ^ 

the  observed  vcilues.  The  latter  is  ciotained  by  the  objective  analysis  and  is  j 

f 

ocrmon  to  all  models.  , [ 

The  contour  lines  are  drawn  using  bilinear  interpolation  on  the  27x29  f 

grid-point  values  of  the  foreccist  error  on  each  level.  The  interval  of  5°C  j 

for  tenperature  is  chosen  for  convenience  and  econony  in  the  prediction  and  f' 

inspection  of  the  maps.  ' 

Figures  1.1  and  1.2  are  samples  of  the  contour  plots  of  errors  of  tenperature.  ' 
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self-e>^lanatx)ry  once  it  is  known  v^iat  infonnation  they  contain.  Each  ' 

contains  tlie  model  nunber  and  three  statistics  for  the  entire  grid  in  the  right 
hand  colimn.  Below  that  are  the  values  and  locations  of  the  maximtin  and  minimun. 

Next  are  the  contour  symbols  and  belcw  that  the  contour  incranents.  Tie  symbols 
and  incranents  are  the  same  for  all  the  plots  in  each  variable.  The  range  of 
forecast,  the  hour  and  date  of  the  forecast,  and  the  level  in  the  vertical 
appear  on  the  top  of  the  plot  under  the  header. 

1.2.2  Technique  of  analysis 

I 

1.2. 2.1  General  contour  features  and  qualities  i 

1 

Concentric  positive  or  negative  contours  do  not  necessarily  indicate  ' 

contours  of  different  absolute  value.  In  some  cases  the  inner  contour  is  of 
the  same  absolute  value  (as  measured  fron  the  zero  level  contour)  as  the 

adjacent  outer  ccaitour.  This  occurrence  was  labeled  with  the  term  "volcano",  | 

while  the  more  conventional  occurrence  of  a higher  value  for  the  inner  contour 
was  labeled  as  a "peak"  or  "hill". 

Cotparative  estimates  of  relative  height  were  determined  primarily  by  j 

counting  ccrcentric  oontours.  This  was  known  to  be  not  entirely  accurate,  but  j 

in  the  few  cases  checked  from  the  table-printout,  the  correlation  was  suffic- 
iently high  to  enable  meaningful  conclusions  to  be  drawn.  * 

Relative  gradients  were  also  determined  by  counting  concentric  oontours, 

as  well  as  examining  contour  spacing.  Thus, if  two  first-increment  oontours  ! 

* 

were  nearly  the  same,  and  each  had  a second  inner  ccaitour,  but  one  had  also  a 1^, 

third  inner  contour,  then  this  latter  structure  would  be  described  as  having  r** 

a higher  relative  gradient.  Again,  the  lack  of  proper  contour  symbol  identi- 
ficaticxi  caused  a certain  degree  of  uncertainty. 
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1 . 2 . 2 . 2 Measures 


Hie  analysis  was  separated  into  two  modes.  The  first  was  a measure  of 
deviatic*!  of  the  individucil  model  forecasts  from  the  observation.  The  second 
was  a ccrparison  of  patterns  among  the  several  models. 

The  deviation  measure  was  defined  so  as  to  indicate  the  presence  of  a ' j 

specified  range  of  error  in  a particular  region  of  tfie  datiain.  The  range  vAiich 
was  selected  corresponded  to  values  sixdi  that  no  second  increment  contour  would 
be  evident.  A quadrant  was  selected  as  the  unit  region.  Since  such  a measure 
corresponded  to  the  absence  of  the  second  concentric  contour,  initial  recordings 

reflected  this.  Each  quadrant  was  later  ccrnpared  with  the  tabular  printout  to  i 

! 

verify  the  presence  of  "volcanoes"  (See  preceding  for  term  definition) . Iheir 
frequency  of  occurrence  was  subsequently  recorded.  i 

This  measure  of  deviation  (D-S)  was  intended  to  provide  an  overall  indi- 
cation of  "goodness"  of  prediction  in  a given  quadrant  because  it  measures  j 

directly  the  amount  of  error  present  in  the  measured  quantity.  | 

Ihe  discrimination  of  patterns  is  acocitplished  by  first  ocmparing  specific  j 

contours  or  contour  grovpings  as  to  similarity  of  location.  Next,  the  peaks  I 

are  compared  as  to  height,  quantity,  location,  and  gradient.  Different  | 

categories  describe  each  of  thefee  attributes.  Thirdly,  included  in  each  | 

category  is  an  allcvranoe  for  certain  of  these  qualities  to  vary,  and  a notation  t 

I' 

/ I 

is  provided  for  recording  this. 

Hius  the  categorizatiOTi  scheme  consists  of  three  parts.  The  first  is  a I 

contour  similarity  ocmparison  primarily  as  to  location.  Second  is  a peak 

similarity  cxmparison  as  to  height,  quantity,  location,  and  gradient.  Hurd  is  > 

the  gray  sc:ale  comparison  that  delineates  vdiich  category  has  more  of  a given 
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contour  level  or  higher  peaks  (larger  gradients) . This  categorization  schane 
is  sunmarized  in  Similarity  Coiparison  Categorization  Table  1.1. 


1.2.3  Results  of  analysis 

Table  1.2  presents  the  fraction  of  the  so-called  "good"  forecasts  of  each  of 
the  six  models.  A "good"  forecast  is  one  that  has  no  second  contour  on  either 
side  of  the  zero  ccaitour,  that  is,  one  with  no  error  greater  than  or  equal  to 
10°C  in  magnitude.  Ihe  fraction  is  the  ratio  of  the  nanber  of  quadrants  with 
good  forecasts  to  the  total  number  of  quadrants  analyzed.  The  figures  in  the 
table  were  obtained  after  volcanoes  and  hills  in  the  contour  maps  had  been 
distinguished  by  referring  than  to  the  corresponding  tabular  printouts  of  the 
grid-point  values  of  the  error. 

Analysis  of  the  oorparative  structures  of  the  contour  plots  among  the 
Vcirious  model  grotpings  vqs  carried  out  utilizing  the  categorizaticn  detailed 
earlier  in  Section  1.2. 2. 2.  The  number  and  distribution  of  forecasts  and  quadrants 
analyzed  were  those  listed  in  Table  1.2.  Various  qualities  were  examined  and  the 
results  are  discussed  belcw. 

1.2. 3.1  ff>dels  "C",  "H"  and  "V" 

Ttve  ocrrparisons  among  models  "C",  "H",  and  "V"  showed  a significant  lack 
of  differences.  In  only  four  quadrants  at  one  level,  specifically  a 12  hour 
night  forecast  at  300  meters,  it  was  not  possible  to  assign  the  S-NI  descriptor 
to  the  full  grouping.  Oatibining  this  information  with  the  closeness  of  the 
fractions  of  "good"  forecasts  ctotained  in  Table  1.2,  we  conclude  that  the  three 
models  are  insignificantly  different. 

It  is  also  clear  from  Table  1.2  that  this  group  out-performed  any  other 
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model  in  all  but  one  category,  the  24-hour  night  forecast.  In  the  latter  model  > 

"R"  is  seen  to  be  best  of  all  and,  model  "7"  is  as  good  as  the  groip. 

1.2. 3. 2 Model  6 

The  contours  of  model  "6"  v^e  found  to  show  a bias  toward  positive  for 

both  12-hour  night  forecast  and  the  24-hour  day  forecast  vhen  corpared  with  the 

grotp  of  "C" , "H" , and  "V" . Ootparisons  with  model  "0"  for  these  forecasts  also 

showed  the  same  bias.  The  biais  is  found  primarily  at  levels  above  the  surface.  j 

At  the  svirface,  the  bias  was  sometimes  discernible,  but  only  infrequently.  In 

order  to  c±»serve  the  bias  properly  it  was  generally  necessary  to  have  the  (S-G)  ! 

I 

descriptor  applicable.  However,  the  corparisons  at  surface  were  often  made  . | 

with  the  (S-C)  or  (S-NI)  descriptors.  | 

I 

No  bias  was  (deserved  in  the  12- hour  day  forecast.  The  bias  in  the  24-hour 
night  forecast  was  at  best  erratic,  none  discernible  at  surface  and  most  foiand 
at  900  and  1600  meters.  However,  the  usage  of  the  (S-G)  oorparison  was  quite  | 

widespread  in  this  forecast  category.  ! 

I 

Presence  or  absence  of  the  bias  appears  to  correspond  well  with  the  values 
of  the  fraction  of  "good"  forecasts  in  Table  1.2.  Only  in  the  12-hour  day  fore-  j 

cast  is  it  found  that  tlie  performance  of  model  "6"  is  as  good  eis  the  other  | 

models,  while  it  is  significantly  inferior  to  than  in  all  of  the  other  cate-  ' 

I 

gories.  , 

1.2. 3. 3 Similarity  corparisons 

i. 

ConpariscaTS  of  models  "O"  and  "6"  with  the  group  of  "C",  "H",  and  "V"  ^ 

were  made  cm  the  bcisis  of  contour  pattern  by  counting  the  nunber  of  quadrants  'f 

in  which  carpeurison  S-C  or  better  was  applicable.  The  results  are  simnarized 


46 


I 


in  Table  1.3.  Distinctions  are  drawn  among  the  cases  if  an  entire  forecast 
ocmtaining  16  quadrants  falls  in  the  category  (forecasts) , if  an  entire  level 
with  4 quadrants  on  it  is  in  the  category  or,  if  neither  of  the  above  is  true 
(quadrants) . 

It  is  obvious  that  the  contrast  in  the  counts  between  models  "O"  and  "6" 
attplifies  vdiat  has  been  found  in  the  deviation  analysis.  Exc^t  for  the  12  hour 
day  forecast  in  which  model  "6"  exhibits  no  bias,  all  categories  show  that 
model  "O"  is  a great  deal  similar  in  contour  pattern  to  model  "6". 

When  the  statistics  on  the  two  different  attributes  are  in  agreement,  they 
reinforce  our  confidence  on  the  visual  inpressions.  On  the  other  hand,  as  in 
the  oatparisons  of  models  "7"  and  "R",  it  was  quite  often  that  no  unanimous 
inference  could  be  drawn  from  the  two  analyses. 
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Figure  1.1  Contour  maps  of  tanperature  forecast  errors,  12  hr  forecasts 
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Figure  1.2  Contour  maps  of  temperature  forecast  errors,  24  hr  forecasts 
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2.0  AUTOMATED  WEATHER  SYSTEM  DEVELOPMENT 

2.1  Introduction 

The  advent  of  low-cost,  reliable  microprocessing  technology  cov5)led 
with  the  pressing  need  within  the  U.S.  Air  Force  for  a modernized  c^jerational 
veather  support  capability  has  spawned  a major  program  of  research  and 
develcpnent.  The  overall  goal  of  the  program  is  the  development  of  a fully- 
autcmated  airfield  weather  observation  and  forecasting  system  vdiich  will 
fulfill  the  v^ther  support  requirements  of  today's  and  tonorrow's  Air  Force 
v^le  utilizing  considerably  fewer  people  than  does  the  present  system. 

Integral  to  this  goal  is  the  solution  of  several  difficult  and  heretofore 
unsolved  autonated  atmospheric  sensing  problems;  namely,  present  weather 
determination,  cloud  ceiling  height,  cloixi  cover  and  glidescopje  slant  visual 
range.  Mditionally,  the  future  needs  of  the  Air  Force  include  objective  and 
autonated  methodologies  (algorithms,  models,  etc.)  to  generate  terminal  area 
weather  predictions  of  aviation-critical  elements  (cloud  ceiling,  runway  and 
slant  visual  range  and  landing  zone  wind  conditions)  for  forecast  intervals 
ranging  fron  a few  minutes  to  2-3  hours.  Within  the  period  of  performance  of 
this  contTc-  ':t  effort,  pjarticular  attention  has  been  paid  to  the  development  of 
automated  procedures  to  specify  the  type  and  intensity  of  surface  weather 
phencmena  (i.e.,  present  veather  determination)  and  to  the  development  of 
indirect  methods  to  specify  and  predict  landing  zone  slant  visual  range. 

These  studies  will  be  described  in  subsequent  sections  of  this  Chapter. 

An  R&D  undertaking  of  this  type  required  the  collection  of  atmospheric 
data  representing  widely  varying  weather  conditions.  A Weather  Test  Facility 
(WTF)  was  established  by  the  Air  Force  Geophysics  Laboratory  on  Otis  AFB, 

Capie  Cod,  Massachusetts,  for  tliis  purpose  in  late  1975.  It  is  ideally  situated. 
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being  in  an  area  frequented  by  considerable  low  cloudiness,  coastal  fog  ard 
a full  range  of  precipitation  forms  ard  is  located  on  an  open  ard  flat  area 
well  removed  from  the  airfield's  runways.  It  is  described  in  detail  in  the 
next  section. 

2.2  Weather  Test  Facility  Site  Description 

The  Weather  Test  Facility  (WTF)  at  Otis  AFB,  MA,consists  of  an  array  of 
instnmented  meteorological  towers  ard  grourd  stations  distributed  in  the 
manner  shown  in  Figure  2.0.  While  there  are  several  one-of-a-kird  sensors 
deployed  in  the  WTF,  there  are  certain  other  sensors  vAiich  are  deployed  at 
multiple  locations.  These  include  Climatronics  Mark  I wind  sensors,  BG&G 

forward  scatter  meters  (FSM)  and  BG&G  terperature-devipoint  sensors  vdiich  are  ' 

i 

placed  on  the  towers  at  levels  approximately  50  ft  apart.  The  central  tower  j 

cdtplex  consists  of  three  towers  extending  to  20  ft,  100  ft,  200  ft  (denoted  ' 

by  C,  B,  and  A in  Figure  2.0) . Hiey  are  configured  to  simulate  an  aircraft  ^ 

approach  zone  and  to  provide  continuous  weather  observations  at  Category  I and  ! 

II  decision  height  levels  (200  ft  and  100  ft  respectively)  essential  to  the 
evaluation  of  glidescope  slant  visual  range  and  low-level  wind  shear  techniques.  | 

TV®  additional  100  ft  tcwers  (labelled  P and  Q in  Figure  2.0)  are  placed  about  J 

1500  ft  either  side  of  the  simulated  approach  zone  oorplex  for  the  purpose  of  ' 

f 

providing  predictor  data  for  the  offset  tower  techniques  for  glidescope  slant  ^ 

1 

visual  range  specification  and  prediction. 

Table  2.1  lists  the  various  sensors  deployed  in  the  WTF,  the  atmospheric  i 

variable  they  are  designed  to  sense,  their  respective  location (s)  aixi  the  [ 

nutiber  of  sensors  at  each  location.  The  placement  cf  weather  sensors  at  ’ ■ 

F 

\ 

location  X (the  ground  instrumentation  site  in  Figure  2.0)  is  for  the  purpose 
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THE  AFGL/LYU  WEATHER  TEST  FACILITY  LOCATED  AT  OTIS  AFB 


of  interccnparisons  among  sensors  designed  to  measure  the  same  atmosjAieric 
variable  sixdi  as  extinction  coefficient,  wind  speed  and  direction,  etc.  A 
short  description  of  the  sensing  principle  enployed  in  each  instrvment  follows. 


Hie  BG&G  Model  207  Forward  Scatter  Visibility  Meter  (FSM)  is  a single- 
ended  device  which  measures  the  amount  of  light  scattered  forward  into  the 
detector  through  a sairpling  volute  of  0.05  m^.  It  has  an  effective  visual 
range  of  0.06  km  to  greater  than  16km. 

The  transnissoieter  measures  the  percentage  of  trananitted  light 
detected  fron  a projector  light  source,  a specified  distance  away.  Three  of 
the  four  units  in  the  WTF  are  standard  AN/GM3-10  units  used  the  Air  Vfeather 
Service.  The  fourth,  located  at  X,  has  a solid  state  modification  kit 
installed. 

The  Iirpulsphysik  Videograph  is  another  single-ended  visibility  device 
which  measures  the  amount  of  light  scattered  backward  into  the  detector  fron  a 
scattering  volune  which  extends  over  a range  frcm  3 to  25  m fron  the  detector. 

The  MRI  Integrating  Nephelctneter  measures  the  reduction  to  visibility 
due  to  suspended  particulates  (haze,  stoke,  and  dust)  by  drawdng  dried  air 
samples  through  a chamber  where  the  scattered  light  from  a pulsed  flash  lamp 
is  detected  and  compared  to  calibrated  reference  scattering  fron  clean, 
filtered  air. 

The  output  of  each  of  the  visibility  devices  (FSM,  transmisscmeter, 
videograjh  and  nefhelcmeter)  is  converted  to  units  of  extinction  coefficient 
for  interccmparison  purposes.  Extinction  coefficient  is  prcportional  to 
visibility  through  Allard's  law  for  night  time  conditions  and  Koschnieder's 
law  for  daytime  conditions. 

Three  different  Eppley  radiometers  are  utilized  at  the  Otis  WIF.  The 
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pyrananeter  cxmsists  of  a vare-wDiond  thermcpile  vrfiich  measures  total  sun  and 
s]cy  incaiiing  radiation.  The  normal  incidence  pyrhelianeter,  mounted  on  an 
equatorial  nount  to  track  the  sun,  measures  direct  solar  intensity.  The  so- 
called  infrared  radiometer  (pyrgecmeter)  measures  unidirectional  terrestrial 
long-wave  radiation. 

The  liminance  meter  is  a laboratory  version  of  a simple  light  meter 
vAiich  measxires  backgroiand  luninance  in  a narrow  field  of  view.  As  d^loyed 
at  Otis,  it  is  mounted  at  the  100  ft  level  of  tcwer  A and  pointed  towards  the 
horizon  along  the  simulated  glide  scope. 

The  Climatronics  xMark  I cup  and  vane  wind  set  measures  wind  speed  by  a 
light  chopping  principle  and  wind  direction  using  a non-contacting  transducer. 

It  features  a light  weight  construction,  lew  power  consumption  and  has  a very 
lew  tlireshold  speed  (0.43  knots) . 

The  J-Tec  vortex  anenoteter  measures  wind  speed  by  measuring  the  spacing 
between  vortices  formed  ty  the  movement  of  air  past  an  obstruction  ty  means  of 
an  ultrasonic  beam.  Wind  direction  is  obtained  frem  a moving  vane  using  a 
potentiometer . 

The  "tosemount  orthogonal  airspeed  sensor  measures  airspeed  vectors  by 
developing  differential  pressiares  after  air  flow  through  pressure  sensings  ports 
located  along  the  cylindrical  body  of  the  device.  It  has  a no-^ioving-parts 
design  and  has  a heating  elanent  built  into  the  sensing  head  to  insure  an  all- 
weather  c^iability. 

The  Gill  prepvane  is  an  instrument  vMch  records  azimuth  angle  for  wind 
direction  and  the  rotational  speed  of  the  propellor  for  wind  speed. 

The  BG&G  temperature  set  utilizes  platinum  resistance  thermometers  to  an 
accuracy  of  - 0.3C.  Ihe  dewpoint  set  uses  a thermo-electrically-cooled 
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TABIE  2.1 

Otis  WTF  Weather  Sensor  Inventory 


Sensor  Type 

Variable 

Location  (Number) 

BG&G  FSM 

Extinction  Coeff. 

A(5), 

B(3),  C(2), 

P{3), 

Q(3),  X(l) 

Trcinsmisscmeter 

Atm.  Transmission 

A(l), 

B(l),  X(2) 

Videografdi 

Extinction  Coeff. 

X(l) 

Nefheloneter 

Extinction  Coeff. 

X(l) 

Pyrancmeter 

Solar  radiation 

X(2) 

Pyrhelioneter 

Direct  sun  radiation 

X(l) 

Pyrgeoneter 

Terrestrial  rcidiation 

X(l) 

Luninance  meter 

Background  luninance 

X(l) 

Climatronics  Cup  and  Vane 

Wind  Speed  + Direction 

A(4), 

P(2),  X(l) 

J-Tec  Vortex 

Wind  Speed  + Direction 

X(l) 

Rosemount  Orthogonal 
Airspeed 

Wind  Speed  + Direction 

X(l) 

Gill  Prc^jvane 

Wind  Speed  + Direction 

X(l) 

BG&G  Thermistor 

Dry  bulb  tamperature 

A(3), 

X(l) 

BG&G  Cooled  Mirror 

Dewpoint  temperature 

A(3), 

X(l) 

Tipping  bucket  reiingage 

Precipitation  Amount 

X(l) 

r 

i 

I 

1 
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detection  mirror  held  at  the  dewpoint  (or  frost  point)  and  a condensate-  ' 

detecting  optical  system  to  precisely  determine  the  dewpoint  tenperature. 

The  Belfort  tipping-buc)cet  rain  guage  records  precipitation  rate  and 
amount  using  a 0.005  inch  bucket.  The  accimulation  of  bucket  tips  over  a 
specified  period  corresponds,  therefore,  to  the  precipitation  amount. 

The  weather  data  vAiich  are  continuously  and  autonatically  gathered  at 
the  Otis  WTF  is  done  so  means  of  a Doric  200-channel  data  logger  system. 

Each  sensor  is  polled  in  the  system  5 times  per  minute  (ten  times  in  the  case 
of  winds)  and  records  on  7-track  magnetic  tape  the  voltage  output  received 
fron  each  device.  One  magnetic  tape  is  written  per  day;thus  weekend  recordings 
cire  paced  at  one  scan  per  minute  to  ensure  the  tape  lasting  frcm  Friday 

afternoon  until  Monday  morning.  The  tapes  are  routinely  transmitted,  by  mail,  I 

to  APGL  at  Hansccm  AFB  vrfiere  detailed  processing  is  conducted  on  the  CDC  6600  1 

to  monitor  the  daily  performance  of  the  various  sensors,  to  prepare  data  for  j 

subsequent  analysis  and  to  assist  in  the  analysis  of  significant  weather  events.  i 

I 

I 

2.3  Data  Reduction  and  Analysis  ' 

i 

fi  lengthy  CDC  6600  carputer  program  called  "REQUOTE"  was  developed  to  i 

convert  raw  signal  voltages  fron  the  nvxnerous  weather  sensing  instnments  to 

laseful  meteorological  data.  The  program  has  evolved  through  several  stages  and  ^ 

t 

has  recently  achieved  expansion  to  include  its  sixteenth  instrument  group  ' i, 

f. 

additicxi.  : 

A major  revision  of  the  original  logic  was  made  necessary  in  early  1977  t . 

vAien  a new  data  recording  systan  replaced  the  original  design.  Repeatedly  I 

since  its  inception,  aciditions  and  revisions  have  been  necessary  as  new 

. i 

instruments  have  been  added  to  the  experiments  or  prc*)latis  develcped  in 
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existing  ccnponents. 

Each  voltage  received  from  the  various  sensors  is  converted  into 
iteteorological  form  and  averaged  over  the  minute's  sairpling.  These  one  minute 
averages  are  printed  out  in  sequenced  pages  as  shewn  in  Figure  2.1  to 
Figure  2.11.  One  page  r^resents  one  hour's  worth  of  data  and  the  groupings 
of  sensors  per  page  are  arranged  for  ease  in  evaluating  sensor  performance, 
technique  develqpnent  and  evaluation  and  in  the  analysis  of  significant  weather 
events.  Printed  sensor  output  is  returned  to  the  Otis  WF  site  personnel  for 
additic^ial  sensor  evaluation  and  monitoring. 

t 

Figure  2.1  shews  the  day-of-year  and  time  of  day  (GMT)  and  the  voltages  I 

recorded  (-99.99  indicates  no  voltage  present  from  that  sensor) . Letters  refer 
to  the  specific  towers  or  ground  locations  shewn  in  Figure  2.0  and  the  numbers,  j 

the  height  above  ground  except  for  numbers  46  to  55  vAiich  are  other  identifica-  | 

tion  nunbers. 

I 

Figure  2.2  contains  the  time  span  covered  by  the  input  data  tape,  the  , 

sunrise-sunset  times  for  the  days  involved  and  the  forward  scatter  meters  used  t 

for  the  slant  visual  range  algorithm  discussed  in  Section  2.5.  Also  included 

1 

are  the  start  and  end  times  of  reduced  visibility  episodes  vAiich  occurred  ’ 

during  this  tape's  time  frame  and  a table  of  the  units  referred  to  in  subsequent  ! 

pages.  f 

' i, 

Figure  2.3  contains  the  extinction  coefficients  calculated  frem  the  r 

r 

18  forward  scatter  meters  and  the  four  trananissaneters  plus  the  ctoserved  and  I 

i- 

estixnated  averages  of  certain  scatter  meters  used  in  slant  visual  range  and 
in  the  determination  of  the  most  prctoable  extinction  coefficient  calculated 
frem  a Markov  prediction  model  being  evciluated  for  short  range  visibility 
prediction. 
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Figure  2.4  details  the  wind  speed  in  neters  per  secxDnd  and  the  wind  | 

direction  of  seven  wind  sets,  the  temperature  and  the  dewpoint  ^read  in 
degrees  Celsius  and  the  rainfall  amount  indicator.  The  rain  code  contains  the 
acconulated  nonber  of  .005-inch  rain  bucket  tips  fron  the  beginning  of  the 
tape  to  the  left  of  the  decimal  point.  The  nunber  of  tips  per  scan  for  a 
given  minute  is  depicted  in  the  five  positions  to  the  right  of  the  decimal 

* 

point. 

Figure  2.5  contains  the  data  from  the  rantiaining  wind  sets  and  the 
observations  fron  the  ranaining  visibility  sensors  and  radiation  meters. 

i 

Figure  2.6  through  Figure  2.10  contain  intercoiparison  data  of  the 
hourly  averaged  wind  output.  Horizontal  and  vertical  profiles  of  wind 
oonditions  can  be  analyzed  fron  these  data  as  well  as  the  relative  performance 
characteristics  of  the  various  wind  sensors  at  location  X. 

Figure  2.11  depicts  the  output  of  the  scatter  meters  used  each  minute 
to  calculate  the  obseirved  and  estimated  values  of  slant  visual  range  using  the 
methods  discussed  in  Section  2.5. 

The  FBQOOTE  program  also  evaluates  the  FSM  data  to  identify  episodes  of 
reduced  visibility  according  to  prescribed  criteria.  In  addition  to  develcping 
^)ecialized  data  listings  for  these  periods,  correction  factors  for  the 
transmisscmeters  are  calculated  and  printed  out  to  aid  in  evaluating  the 
instrunent  drift  during  reduced  visibility  episodes.  This  is  a cotmon  problem 
with  trananisscmeters  which  must  be  accounted  for  prior  to  their  utilization 
in  analysis  and  prediction  studies. 

A separate  subroutine  was  added  to  REX3U0TE  in  mid-1978  to  evaluate  the 
slant  visual  range  estimation  and  short  range  prediction  procedures  previously 
developed.  A control  forecast  technique  which  conputes  the  most  probable 
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extinction  coefficient  using  a well-tested  Markov  model  is  part  of  the  evalua- 
tion. Both  techniques  assume  normality  in  converting  the  deterministic 
predictions  to  exceedance  prc±>abilities. 

Recently,  routines  have  been  added  to  inclixJe  scatter  plots  of  the 
Rosarount,  Gill  and  J-Tec  wind  direction  and  speed  output  against  the  Clitna- 
tronics  wind  set.  Finally,  RBQUQTE  generates  a specialized  data  file  (stored 
on  permanent  file  in  the  CDC  6600  systan)  containing  the  sensor  outputs 
required  by  the  ootputer  program  being  used  in  the  present  weather  autaiation 
studies  discussed  in  Section  2.4. 

The  evolving  nature  of  the  research  requirements  and  sensor  configura- 
tions of  the  Otis  WTF  put  heavy  demands  on  and  add  to  the  ooiplexity  of  a basic 
data  processiip  program  like  REQUOTE.  The  resulting  problems  have  been  further 
corplicated  by  seme  major  setbacks  in  the  data  logging  system  at  Otis.  Special 
data  editing  routines  had  to  be  developed  to  monitor  each  sensor's  output 
voltage  to  ensure  its  validity  before  inclusion  in  subsequent  averaging  and 
conversion  rxjutines.  Criteria  for  data  acceptability  underwent  evolution  as 
experience  was  gained  with  the  data  logging  system.  Other  routines  had  to  be 
developed  i.  ' filter  excessive  noise  from  the  useful  data  stream.  Data  loss 
rates  became  excessive  in  mid-1977  when  the  6600  system  could  not  successfully 
read  magnetic  tapes  being  generated  at  Otis.  After  numerous  attatpts  to 
circumvent  the  problem  with  software  modifications  to  RBQUOTE,  the  problem  was 
traced  to  the  tape  drive  at  the  Otis  WTF  which  was  replaced.  Subsequent 
routine  processing  of  the  collected  data  has  generally  been  without  incidence 
and  timely. 
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2.4  Autonation  of  visual  Present  Weather  Observations 


The  autonated  array  of  weather  instnments  described  previously  provides 
an  ideal  setting  for  the  evaluation  of  techniques  for  discrimination  of  present 
weather  elements.  A decision-tree  approach  was  chosen  to  determine  the  extent 
to  vArLch  an  objective  solution  to  the  problon  ooiiLd  be  made.  Primary  consid- 
eration was  given  to  extracting  present-weather  information  from  a selected 
subset  of  the  weather  instnments  cited  in  Section  2.2. 

The  decision  tree  develcped  under  this  contract  is  shown  in  generalized 
form  in  Figure  2.12.  It  consists  of  two  basic  programs,  (a) , an  ctostructicai 
to  visibility  determination  and,  (b) , a present  weather  determination.  One  of 
the  more  obvious  decision  points  in  the  first  of  these  programs  is  the  test  to 
determine  if  a restriction  bo  visibility  exists.  The  forward-scatter-meter 
provides  this  information.  The  next  test  determines  the  hunidity  of  the  air. 

If  low,  then  the  program  branches  to  the  "dry"  obstructions  and  with  ir^nit 
fron  the  nep^ielcmeter,  snow  conditions  and  wind  gustiness  discriminates  between 
stoke  (K) , haze  (H)  , dust  (D) , or  blowing  snow,  sand  or  dust  (BS,  BN,  or  BD)  . 
If  the  hunidity  is  high,  then  the  program  leads  to  the  fog  branch  and  deter- 
mines ground  fog  (GF) , fog  (F)  or  ice  fog  (IF)  through  examinations  of  the 
vertical  variations  of  visibility  and  tarperature. 

The  present  weather  determination  is  triggered  by  a positive  report 
fron  the  tipping-bucket  rain  guage.  The  frequency  of  tips  is  related  bo 
intensiby  and  is  used  to  separate  the  precipitation  into  three  classes, 
light  (-) , moderate  ( ) , and  heavy  (+) . Tarperatirre  tests  separate  freezing 
precipitation  frcm  liquid  and  differences  between  the  backscatter  meter  and 
forward-scatber-meter  discriminate  snow  (S)  fron  raiin  (R) , freezing  rciin  (ZR) 
eind  freezing  drizzle  (ZL) . Ehpirical  evidence  which  shows  that  the 
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transmissoneter  and  FSM  respond  similarly  in  drizzle  (anall  drops)  but  depart  ' i 

significantly  in  rain  (large  drcps)  is  used  to  separate  rain  (R)  from  ] 

drizzle  (L)  observations.  • ! 

Hie  model's  algorithms  produce  minute-by-minute  observations  of  weather.  ) 

These  in  turn  are  being  evaluated  as  part  of  an  AFGL  research  study  against 
the  operational  weather  observations  routinely  made  by  FAA  control  tower 
personnel  at  Otis  AFB.  The  control  bcwer  is  located  approximately  one  mile 
east  of  the  Weather  Test  Facility.  Figure  2.13  is  an  example  of  one  hour's 
data  considered  by  the  decision  tree,  the  present  weather  it  deduced  and  the 
cdiparative  hunan  ctoservations  for  the  same  period  on  8 June  1978. 

Column  1 denotes  the  time  of  observation  to  the  nearest  minute,  Oolixtins 
2 through  6 are  one-minute  averages  of  extinction  coefficient  (x  10““^  m”^)  j 

measured  by  the  FSM  at  five  levels  of  tcwer  A (200  feet,  150  feet,  100,  50  feet  | 

and  at  ground  level) . Column  7 is  extinction  coefficient  obtained  from  the  ■ 

transnissaneter.  Colimns  8 through  17  shew  the  one-minute  averages  of  j 

tanperature  (xlO  degrees  C)  at  a given  level  of  tcwer  A (200,  150,  100,  50  feet 
and  ground  level)  followed  by  the  tarperature-dewpoint  difference  at  the  same 
height.  Column  18  is  a one-minute  count  of  the  number  of  tips  made  by  the 

I 

rain  guage  in  response  to  precipitation.  Each  count  represents  .005  inches  of 
precipitation.  Thus  the  intensity  is  easily  obtained.  The  next  oolirm,  19,  ! 

shews  the  extinction  coefficient  c±)tained  from  the  nephelcmeter . Colum  20  ^ 

displays  the  extinction  coefficient  observation  of  the  backscatter  meter. 

I 

Colvxm  21  shews  the  algorithm's  determination  of  present  weather  and  obstruc- 
tion  to  vision.  The  last  column  shows  the  present  weather  observations  taken 
by  the  FAA  ctoserver  during  this  hour.  The  symbol  T refers  to  a thunderstorm 
occurrence.  Other  symbols  shewn  here  denote  rain  (R) , fog  (F)  and  drizzle  (L) 
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with  the  intensity  vciriations  mentioned  above. 

Preliminary  results  of  the  cottparisons  between  the  present  weather 
program  and  the  FAA  observations  have  been  very  encouraging.  The  technique 
has  proved  to  be  very  responsive  to  rapidly  changing  weather  events  as  illus- 
trated by  the  variations  of  precipitation  during  this  particular  hour.  Visi- 
bility variations  are  as  quickly  observed.  The  importance  of  quickly  deter- 
mining significant  changes  in  weather  for  aviation  needs  cannot  be  overerpha- 
sized. 

2.5  Remote  Ibwer  Slant  Visual  Range  Studies 
2.5.1  Background 

As  the  cperational  minima  for  low  visibility  approaches  and  landings  are 
reduced  at  certain  airfields  due  to  improved  autcmated  landing  systans  being 
installed,  the  requiranents  for  accurate  and  timely  measurements  of  airfield 
visibility  are  beocming  more  urgent.  Currently,  information  sujplied  to  the 
pilot  consists  of  hunan  observer  estimates  of  prevailing  visitjility  and  instru- 
mented measurements  of  runway  visual  range  (RVR)  at  one  or  more  locations 
adjacent  to  the  runway.  Although  these  observations  supply  essential  informa- 
tion, experience  heis  shewn  that  the  surface  RVR  measurements  and  general-area 
visibility  ctoservations  are  often  not  representative  of  the  conditions 
encountered  by  the  pilot  as  he  proceeds  along  the  glide  path  through  the  deci- 
sion height  (DH)  to  touchdown. 

The  development  of  an  effective  technique  for  the  measurement  of  slant 
visual  range  (SVR)  ahead  of  the  aircraft  alor>g  and  belo/  the  glide  slope  path 
has  been  the  primary  objective  of  this  research.  Field  experiments  at  the 
Otis  Weather  Test  Facility  (WTF)  have  been  conducted  using  an  instnmented 
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tx:v^r  aj^roach.  Sunmarizmg , tJie  study  seeks  to  determine  the  accuracy  with 
vrfiich  SVR  in  a simulated  approach  zone  could  be  specified  at  selected  time 
intervals  frcm  time  zero  to  plus  ten  (10)  minutes  using  various  ccmbinations 
of  forward  scatter  meter  (FSM)  measurements  made  at  time  zero  at  a tower 
1500  ft  offset  frcm  the  simulated  approach  zone. 

2.5.2  Initial  data  analysis 

Five  episodes  of  reduced  visibility,  totalling  9480  minutes,  were 
selected  to  serve  as  a data  base  for  this  study.  The  March  1976  episode  was 
used  as  a dependent  set  with  the  other  four  episodes  ccrprising  the  independent 
set.  Coastal  advection  fog  occurred  in  the  March  1976,  May  1976  and  July  1976 
episodes.  The  August  1977  episode  was  a radiation  fog  episode. 

The  dependent  data  set  was  exhaustively  examined  in  order  to  define  the  I 

approach  to  be  taken.  A program  SCATTER  was  developed  to  evaluate  the  linear 
relationships  among  those  FSMs  positioned  at  identical  heights  on  towers  A and 

I 

Q above  ground.  The  examination  of  the  statistics  frcm  SCATTER  revealed 
systenatic  errors  among  the  FSMs,  the  result  of  a gradual  accumulation  of  a 
residual  film  on  the  lens  of  the  FSMs.  Correction  factors  were  calculated  and 
applied  to  the  d^sendent  data  set  thereby  minimizing  this  systematic  error. 

A similar  process  was  applied  to  the  independent  data  set.  , 

A tool  which  proved  very  helpful  in  editing  and  validating  the  basic  ' ' 

data  sets  was  the  Interactive  Graphics  device  available  through  the  AFGL 
CDC  6600  system.  The  graphics  software  package,  provides  options  to  display 
edit,  on  a CRT,  user  selected  subsets  of  sensor  measuranents  collected  at 
WTF.  Ihe  software  can  display  up  to  360  data  points  (typically  6 hr  of 
data  in  our  case)  of  one  to  four  sensors  at  any  one  time.  As  simple  instructions 
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are  displayed  in  the  CRT,  the  liser  interacts  with  the  CDC  6600  system  using 
either  a light  pen  or  keyboard.  A hard  ccpy  pen  plot  option  is  also  available. 

The  data  generated  at  the  Otis  WTF  was  frequently  examined  using  the 
interactive  graphics  device.  In  particular,  portions  of  the  five  data  sets 
were  edited  and  hard  copies  made  through  this  device.  An  indepth  description 
of  the  procedures  vdiich  ccnprise  the  software  package  can  be  found  in 
Appendix  2A. 

As  illustrated  in  Figure  2.0  (Section  2.2),  the  Otis  WTF  includes 
transnisscmeters  mounted  along  a slant  path  of  about  12  degrees  from  the 
horizontal  between  towers  A and  C.  They  were  so  positioned  to  provide  SVR 
measinrements  appropriate  to  a 100  ft  decision  height.  The  initial  formulation 
for  these  studies  was  to  incorporate  the  slant  transmisscmeter  measurements  into 
the  predictand  definition.  Unfortunately,  the  slanted  position  of  the  trans- 
misscmeters  often  resulted  in  extinction  coefficient  values  vdiich  vere  non- 
representative  of  actual  conditions  during  periods  of  precipitation.  Liquid 
water  on  the  detec:tor  surface  and  that  which  accumulated  at  the  base  of  the 
detector  hood  undoubtedly  contaminated  the  measuranents . A program  was  written, 
called  TRANSCDR,  to  evaluate  the  transmisscmeter  data  against  the  tcwer-mounted 
FUM  output  during  selected  episodes.  An  analysis  of  its  output  confirmed  the 
bias  towc'rds  high  extinction  coefficient  values  frcm  the  transmisscmeters. 

2.5.3  Analysis  technique  formulation 

An  alternative,  and  equally  valid,  approach  bo  defining  the  predictand 
for  the  study  was  then  formulated.  An  unweighted  average  of  the  FSM 
measurements  at  tcwer  A up  to  the  100  ft  level  was  conputed  bo  represent  the 
SVR  in  the  simulated  approach  zone: 
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^ ^^100  ^0 

where  siabscripts  100,  50,  and  10  refer  to  the  height  of  the  FSM  in  feet. 

Specific  conbinations  of  Q^qO'  ^50'  ^10  ^10  selected  as  predic- 

tors of  A.  The  series  of  predictors  were  chosen  to  investigate  specification 
accuracy  based  on  the  nutnber  and  spacing  of  the  instruments  as  follows: 

1)  A = (Q^qq  + Q^q  + Qio)/3.  + 6 

2)  A = Kj^Q5q  + K2Aj^q  + 6 

3)  A = Aj^q  + 6 

Seme  sirtplifying  ass\jmptions  were  made  to  establish,  in  advance,  the 
specification  algorithms.  Method  (i)  assumes  horizontal  hemogeneity  of  fog 
density  over  the  1500  ft  distance  between  towers  A and  Q.  Methcd  (2)  is  a 
prediction  algorithm  with  estimates  of  coefficients  and  K2  obtained  through 
nwltiple  linear  regression  techniques  using  the  dependent  data  set.  Method  (3) 
is  merely  a representation  of  the  current  operational  methodology  vrfiich  uses 
the  runway  dsual  range  (RVR)  to  aid  the  pilot  in  his  landing  decision.  For 
the  purposes  of  this  study,  it  represents  a control  technique  against  vdrLch 
the  two  offset  tower  methods  (1  and  2)  will  be  evaluated.  Methods  (1)  through 
(3)  provide  detenru nistic  solutions  for  SVR  specification  and  short  range 
prediction.  Transformation  to  exceedance  probability  form  was  pursued  throi>gh 
two  approaches.  The  first  model  was  defined  as  follows: 

SVRrp  — SVR|j| 

^ RMSE 
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vtere  p = equivalent  normal  variate 

SVRp  = SVR  threshold  (daytime:  *5  mile  : SVRji  = 46) 

SVR|j,  = SVR  at  t = <(1 

RMSE  = Root  Mean  Square  Error 

RMSE  was  enpirically  determined  to  be  a value  of  20  in  extinction  coefficient 
units.  The  second  model  was  defined  as  follows: 


vAiere  p = equivalent  normal  variate  ^ 

f^  = SVRp  equivalent  normal  variate  . ; 

fo  = SVRjji  equivalent  normal  variate  ; 

R = hour-to-hour  correlation  coefficient  of  equivalent  normal  variates 
Note  that  the  dencminator  is  the  standard  error  of  estimate  (SE) . A set  of 
mutually  exclusive  and  exhaustive  equations,  fined  tuned  for  time  of  day  and  f 

season,  cure  used  to  relate  SVRp  (extinction  coefficient)  to  f^  (equivalent  * 

normal  variate) . 

2.5.4  Results 

The  three  specification/prediction  methods  were  applied  to  each  of  the 
four  ind^jendent  data  episodes  on  a minute-by-minute  basis.  Since  fog  density 

i 

tends  to  be  horizontally  uniform  over  distances  of  several  kilcroeters,  oaipari-  ^ 

t 

sons  of  the  specified  SVR,  by  the  two  tcwar  methods,  with  the  ctoserved  SVR  (A) 

reveal  nearly  a one-to-one  correspondence.  Even  though  the  time  variations  of  ■, 

SVR  are  significant  over  periods  as  short  as  several  minutes,  the  offset  tower 

estimates  closely  track  the  c±)served  conditions  1500  ft  away. 

Method  2 vdiich  uses  data  fran  the  50  ft  level  of  tower  Q in  oonbination 
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witli  the  KVR  value  (A^^q)  yields  results  equally  as  good  as  Method  1 vdiich 
draws  on  data  to  the  100  ft  level  of  tower  Q.  The  information  content  in  the 
layer,  surface  to  50  ft,  sufficiently  identifies  the  intensity  of  the  vertical 
gradient  cf  extinction  coefficient  such  that  the  increment  frcm  50  to  100  ft 
is  superfluous.  Meanv^le  both  methods  iirprove  upon  the  control  forecast 
technique  (Method  3)  a significant  amount. 

A convenient  way  to  cotpare  the  deterministic  solutions  of  the  three 
methods  is  to  evaluate  their  skill  relative  to  preselected  thresholds  of 
extinction  coefficient  v^ch  have  operational  significance.  The  threshold 
used  was  12  kra~^  which  corresponds  to  daytime  visibilities  of  one-quarter  mile 
and  nighttime  of  one-half  mile  respectively.  By  suntnarizing  the  results  into 
2x2  contingency  tables,  for  each  threshold  separately,  one  can  ccnpute  the 
Probability  of  Detection  (PoD)  for  each  method  and  for  each  independent 
episode.  Table  2.2  surmarizes  the  PoD's  for  SVR  specification  (0  lag)  and 
predictions  ranging  frcm  2 to  10  min.  With  the  exception  of  the  August  1977 
radiation  fog  episode,  the  PoD  exceeds  90%  for  both  Method  1 and  Method  2 
for  the  0 lag  and  2 min  prediction.  The  high  tenporal  variability  of  extinc- 
tion coeff  "ient  beccmes  clearly  evident  by  10  minutes  after  observation  time. 
By  contrast,  however,  the  control  forecast  technique  (Method  3)  has  consider- 
ably less  skill  in  the  advection  episodes.  Since  it  relies  exclusively  on 
the  surface  RVR,  it  consistently  overestimates  SVR  in  radiation  fog  episodes 
v^iich  explains  its  high  PoD. 

Comparison  of  the  two  methods  of  transforming  the  deterministic  solu- 
tions to  exceedance  probabilities  was  achieved  through  the  use  of  the  p-score 
statistic  V(^ch  calculates  the  mean-squared-probability  error.  It  has  the 
characteristic  that  a value  of  zero  relates  to  perfect  estimation  and  a value 
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of  1 to  totally  iirperfect  estimation.  Based  on  the  independent  episodes 
simnarized  in  Table  2.3,  the  second  method,  which  considers  the  equivalent 
normal  variates  of  SVR,  generally  resulted  in  lower  (i.e.,  better)  p-scores, 
especially  for  prediction  times  of  15  minutes  or  greater. 


t. 

ii' 
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TABLE  2.2 


THRESHOIiD  = 

12  km"l 

PRCBABILITY 

LAG 

OF  EffTTECTION  (%) 

(MINUTES) 

MEOTCO 

EPISCOE 

0 

2 

5 

10 

1 

1 APR  76 

AD  FOG 

94 

90 

83 

71 

JULY  76 
RAIN/AD  Fog 

97 

96 

91 

83 

MAY  76 

AD  FOG 

97 

97 

97 

94 

AUG  77 

RAD  FOG 

80 

80 

80 

79 

2 

1 APR  76 

AD  FOG 

91 

88 

80 

69 

JULY  76 
RAIN/AD  Fog 

91 

90 

86 

79 

MAY  76 

AD  FOG 

97 

96 

96 

94 

AUG  77 

RAD  FOG 

78 

78 

78 

78 

3 

1 APR  76 

AD  FOG 

13 

13 

13 

12 

JULY  76 
RAIN/AD  Fog 

56 

56 

56 

54 

MAY  76 

AD  FOG 

87 

87 

87 

87 

AUG  77 

RAD  FOG 

100 

100 

100 

98 
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TABI£  2.3 


IHRESHOLD: 

5 km“^  (day) 

P 

SCORE 

12  km~^  (night) 

LAG 

(MINUTES) 

METOCO 

EPISODE 

2 

5 

15 

30 

60 

1 

APRIL  76 

.065 

.105 

.211 

.259 

.439 

MAY  76 

.006 

.010 

.031 

.058 

.111 

JULY  76 

.014 

.024 

.052 

.071 

.115 

AUGUST  77 

.049 

.053 

.071 

.106 

.187 

2 

APRIL  76 

.066 

.108 

.175 

.225 

.315 

MAY  76 

.007 

.010 

.031 

.056 

.108 

JULY  76 

.013 

.023 

.048 

.063 

.091 

AUGUST  77 

.049 

.053 

.068 

.095 

.150 

i 
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2A  APPENDIX:  Interactive  Graphics  Program  Description 
2A-1  Introduction 

Graphics  is  an  interactive  capability  available  at  Terminal  #2  (Bldg 
1105B,  Roan  156)  of  the  AFGL  CDC  6600  Systan.  Section  17  of  the  AFGL  User's 
Guide  (revised  Jun  1977)  provides  detailed  description  of  program  setup, 
terminology  and  internal  display  control  routine  description. 

The  graphics  systan  is  used  for  LYU  purposes  to  display  a maximum  of 
360  points  - six  hours  of  data  - of  one  to  four  sensors.  This  systan  is 
primarily  used  to  edit  the  data  collected  at  Otis  Vfeather  Test  Facility 
according  to  the  user's  specification.  The  operator  can  change  many  of  the 
values  of  the  variables  in  the  program.  The  user  can  do  this  either  by 
punching  data  in  on  the  available  typewriter,  or  activating  a series  of 

instructions  by  means  of  the  light  pen  - a pen  placed  to  the  right  of  the  j 

console . 

2A-1.1  There  are  four  program  decks  that  are  used  in  conjunction  with  ' 

graphics  - Files,  Ready,  Two  Graphics  decks.  All  of  the  data  for  the  Otis  ■ 

i 

Project  is  either  stored  on  an  archived  tape  or  processed  directly  off  a 

Cap  Tape.  In  either  case,  each  record  contains  a one  minute  average  of  a | 

I 

maximum  of  150  weather  parameters  (Exco,  Speed,  Tatp,  Neph,  etc.)  and  the  time.  j 

Files  is  a program  that  converts  the  data  you  wish  to  edit  into  six  hour  ^ 

blocks  of  one  sensor  to  be  irput  to  graphics.  Three  permanent  files  are  ‘ 

[•' 

t created  in  this  program  - two  random  - mass  storage  files  arranged  in  six-  i 

hour  blocks  (edited  and  unedited)  and  a file  with  the  start  date  of  the  data,  ‘ 

i' 

start  date  of  the  first  block  and  the  number  of  blocks. 

J* 

• 2A-1.2  Ready  is  a program  that  converts  the  six  hour  block  data  of  one  sensor 

into  one  minute  averages  of  one  hundred  fifty  sensors  with  time.  This  new 
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edited  file  can  new  be  lised  as  ii^t  to  any  statistical  analysis  program. 

2A-1.3  Tbe  two  graphics  programs  are  used  to  put  the  graphics  overlays  into 
the  systan.  Ihe  first  program  created  a permanent  file  with  any  changes  the 
user  might  vgant  to  add  to  the  program.  Semetimes  labels  might  be  changed  or 
default  values  may  be  altered.  This  is  done  by  the  delete  and  insert  state- 
ments. The  second  program  attaches  all  four  files  necessary  for  running 
gr^hics.  These  are  the  edited  and  unedited  files  of  data,  the  ident  file 
containing  time  and  nutiber  of  blocs,  and  the  graphics  file  containing  the  user 
changes  and  the  original  graphics  overlays. 

2A-2  Overlay  (0,0) 

The  graphics  program  operates  on  a sysban  of  overlays.  Each  overlay 
deals  with  a specific  area  of  the  program.  There  are  10  overlays  labelled 
from  0 to  9.  The  Overlay  (0,0)  contains  all  the  logical,  conmon,  dimension 
and  subroutines  that  open,  read,  write,  and  close  mass  storage,  set  the  origin 
for  the  screen,  and  set  the  labelling  for  the  hard  copy.  This  overlay  iitmedi- 
ately  goes  into  the  next  overlay  without  the  help  of  a light  button.  A light 
button  is  an  instruction  on  the  screen  that  can  be  activated  by  the  light  pen. 

2A-3  Overlay  (1,0) 

Overlay  (1,0)  is  only  used  at  the  beginning  of  the  gr^^cs  program. 

In  tliis  overlay,  the  start  of  the  data,  the  first  block  time,  and  the  nunber 
of  blocks  for  a specific  case  are  read  into  the  program  from  a card.  The 
beginning  of  the  data  is  put  \ap  on  the  screen.  A cross  always  appears  on 
the  screen  where  the  origin  is  set.  The  origins  (x,  y values)  of  the  four 
possible  graphs  are  placed  in  two  arrays.  The  screen  is  approximately 
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17"  X 17"  using  the  Icwer  left  hand  comer  as  the  origin.  Four  graphs  must 
be  fitted  into  this  frame  of  reference. 

In  the  block  data  statanents,  max  and  min  (default  values)  are  given 
for  batp,  direction,  speed  and  exoo.  At  this  point  in  the  program,  they  can 
be  altered.  Several  statanents  will  flash  on  the  screen.  The  instructions 
are  self-evident.  Type  XXX. XXX.  This  tells  you  to  type  in  two  numbers  with 
a maximum  of  three  whole  nimbers  before  the  decimal  point.  This  is  ths 
instruction  for  the  temp,  direction  and  the  exco  values.  The  instruction  for 
speed  is  XX.XX.  - 2 numbers  for  the  max  and  min  values  of  2 vAx5le  numbers 
before  the  decimal  point. 

There  is  only  one  safeguard  in  typing  in  instructions.  If  wrong  num- 
bers are  punched  in  the  correct  format,  the  program  accepts  the  incorrect  data.  | 

Howsver,  if  either  more  or  less  than  the  correct  number  of  characters  are  { 

punched  in,  the  program  will  reject  it  and  you  will  be  allcwed  to  type  in  ' 

» 

again.  This  overlay  also  specifies  vihether  you  wish  to  plot  the  Exco  values  j 

on  a log  or  linear  scale.  This  is  actually  semi-log  because  the  X-axis  is 
always  linear  (time  every  minute  up  to  six  hours) . 

2A-4  Over  .Lay  (2,0) 

Overlay  (1,0)  iimvediately  goes  to  the  third  Overlay  (2,0),  DATO  (DATA)  . ^ 

This  overlay  sets  up  the  data  used  for  the  first  block  and  all  subsequent  ' ^ 

i. 

I blocks.  There  are  four  possible  options.  These  options  are  activated  by 

placing  the  light  pen  on  the  instruction  at  the  botbcm  of  the  screen  that  says  * 

. i' 

one,  two,  three  or  four. 

In  Option  One,  you  type  in  8 digits  to  specify  4 sensors  you  wish  to 
work  with.  In  the  program,  the  numerical  values  are  translated  into  sensor 
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names.  The  data  for  6 hours  of  each  of  four  sensors  specified  is  read  off  of 
a randan  access  file  that  was  created  before  the  program  was  executed.  The 
index  of  the  file  must  be  corputed  in  the  program  by  the  block  number.  If 
the  index  is  missing,  the  data  are  filled  with  missing  data  9999.  The  labels 
of  increments  for  plotting  are  ccnputed  by  subtracting  the  minimum  value  fron 
the  maximum  value  and  dividing  by  seventeen.  The  time  window  for  the  current 
time  is  put  on  the  screen.  The  time  array  is  corputed  with  a time  associated 
with  every  360  minutes. 

The  randan  access  file  is  arranged  in  the  following  manner.  Each 
record  contains  360  minutes  of  data  of  one  sensor.  Ihere  is  an  overlap  of  60 
minutes  on  each  record.  The  first  hour  of  one  record  and  the  last  hour  of 

the  previous  record  are  the  same.  Therefore  you  really  only  get  five  new  I 

hours  of  the  data  on  the  screen.  j 

In  subsequent  uses,  any  of  the  four  options  can  be  used  depending  on 
the  user's  needs.  Option  One  can  be  activated  if  the  user  wants  to  advance 
a six  hour  period  with  the  same  four  sensors.  The  time  array  is  advanced 
5 hours  because  one  hour  is  overlapped  fron  previous  frame.  The  new  data  are 
read  fron  the  randan  access  file  first  updating  the  block  by  one. 

2A-4.1  Option  TVo  can  be  activated  if  the  user  wants  to  look  at  four 

I 

different  sensors  at  the  same  time  as  was  used  previously.  The  number  of  , 

sensors  are  typed  in.  There  is  added  instruction  of  U or  E.  The  program  is  ^ 

given  twa  identical  randan  mass  storage  files  created  to  be  input  to  graphics.  i 

Oie  of  these  files  U remains  untouched  through  graphics.  The  other  file  ' 

E can  be  edited,  in  the  course  of  using  graphics,  according  to  the  user's 
specifications.  In  this  option  the  scale  for  plotting  can  be  changed  fron 
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linear  to  logarithmic  or  vice  versa. 


2A-4.2  Cation  Three  is  activated  if,  at  any  time,  the  user  wants  to  reset 
the  time  frame  to  any  block  of  data  in  the  rardcm  access  file.  If  the  time 
that  is  typed  in  is  not  the  beginning  of  a block,  the  program  cycles  to  the 
nearest  block.  The  user  types  in  the  8 digits,  e.g.,  02034147  (2  for  each 
sensor) , he  wishes  to  look  at.  Then  he  types  in  10  digits  to  specify  the 
start  time  he  is  interested  in,  2 each  for  month,  day,  year,  hour,  and  minute, 

(e.g.,  1027770830) . The  program  cotputes  the  nearest  block  index  and  reads 
the  appropriate  data.  Then  it  sets  the  correct  time  array  based  on  the  current 
time. 

2A-4.3  Cation  Four  is  activated  if  the  user  wants  to  display  data  (maximum 

six  hours)  within  two  blocks.  The  user  again  specifies  the  four  sensors 

(8  digits)  he  wishes  to  use.  He  types  in  whether  he  wishes  to  use  the  unedited 

or  edited  file.  He  then  types  in  the  initial  time  (10  characters)  and  the  > 

final  time  (10  characters) . The  maximum  span  of  time  frcm  initial  to  final 

time  is  six  hours.  Any  time  span  less  than  this  vrould  be  acceptable.  The 

data  and  t.  ne  arrays  are  appropriately  set.  \"Jhen  this  overlay  is  ccmpleted, 

the  program  goes  immediately  into  the  Eras  Overlay  described  below.  Two  light 

battons  are  set  in  the  Datn  Overlay  - (1)  Eras,  (2)  End.  Since  the  Datn 

Overlay  goes  imnediately  to  the  Eras  Overlay,  these  light  buttons  are  [• 

superfluous. 

2A-5  Overlay  (3,0) 

The  Eras  Overlay  erases  the  screen.  The  program  puts  up  several  light 
buttons.  A light  button  activates  an  overlay.  You  can  activate  a light  button 
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by  means  of  tJie  light  pen.  Hitting  one  of  the  names  of  the  overlays  on  the 
side  of  the  screen  puts  the  operator  into  that  overlay.  The  light  buttons  for 
the  Eras  Overlay  are  Quad,  Datn,  Zoon,  End,  and  Draw.  Eras  also  puts  i:p  the 
time  for  the  beginning  of  the  data  and  either  the  current  time  or  initial  and 
final  time. 

2A-6  Overlay  (4,0) 

Quad  is  an  overlay  that  sets  up  the  sensors  to  be  plotted.  There  are 
four  light  buttons  displayed  on  the  right  hand  side  of  the  screen  - Draw,  End, 

Zocm,  and  Delt.  On  the  left  hand  side  of  the  screen  is  a message  - "Which 

I 

Quadrant".  Underneath,  there  are  two  colonns,  the  left  hand  column  has  four  I 

rows,  each  row  containing  the  name  of  one  of  the  four  sensors.  The  right  hand 
column  contains  the  same  four  rcws  with  a + beside  each  name.  The  fifth  row 

I 

on  this  side  contains  the  ward  "GO".  Each  of  these  words  in  both  columns  is 

! 

activated  by  a light  pen.  There  is  a distinct  difference  between  the  two  ' 

columns.  When  the  message,  "Pick  Quad",  is  flashed  at  the  bottan  of  the  ! 

screen,  the  user  must  make  a decision  as  to  which  column  to  activate.  If  he 

only  wants  to  plot  the  six-hour  block  of  one,  two,  three,  or  four  sensors  he 

has  picked  in  data,  he  may  pick  the  sensors  fron  the  right  hand  column  and 

then  hit  go.  This  will  immediately  erase  the  screen  and  plot  the  sensor  (s) 

the  user  has  selected  (go  to  Draw  Overlay) . j 

t 

However,  if  he  wants  to  do  sonething  else,  he  uses  the  left  hand  column. 

He  activates  the  sensors  he  wishes  to  use.  Then  he  uses  the  light  pen  to 

i 

activate  the  overlay  on  the  right  hand  side  of  the  screen.  He  has  three  real  ; 

choices.  He  may  activate  Zocm  vAiich  will  blew  up  vertically  and  horizontally 
a anall  scale  of  time.  He  may  activate  Delt  vAiich  will  allow  the  user  to  edit 
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the  data  he  is  plotting.  His  third  choice  is  Draw  vdiich  just  plots  the  sensor 
he  picks.  He  can  also  end  the  program  at  this  point  with  the  End  Overlay. 


2A-7  Overlay  (5,0) 

Draw  is  the  overlay  that  plots  the  six  hour  block  (360  points)  of  fron 
one  to  four  sensors,  log  or  linearly.  Several  light  buttons  are  put  up  on  the 
right  hand  side  of  the  screen.  They  are  Quad,  Datn,  Zocm,  Delt,  Stor,  Hard, 
Eras,  and  End.  He  may  decide  to  go  to  any  of  these  overlays  fron  Draw.  Draw 
puts  on  the  screen  whatever  is  in  the  data  array. 

2A-8  Overlay  (6,0) 

Zoan  is  the  overlay  that  sets  up  the  sensor  that  is  specified  for  a 
zooning  routine.  The  y - axis  - in  units  of  the  sensor  is  expanded  to  fill 
the  entire  screen.  The  x - axis  - time  scale  - is  typed  in  during  the  execu- 
tion of  the  Zocm  Overlay.  The  instructions  are  shewn  at  the  bottom  of  the 
screen.  Type  in  XX. XX  (format)  to  specify  hour  to  zocm  (beginning  of  hour 
to  zocm)  . Type  in  XX. XX  (end  hour) . There  is  a line  to  type  in  beginning  and 
ending  hours.  The  time  array  and  data  array  will  be  set  up  for  plotting  in 
the  zocmin  array.  After  all  the  instructions  have  been  satisfied,  this  over- 
lay will  go  immediately  to  the  Draw  Overlay.  Frcm  here  the  user  can  go  to  the 
editing  overlay  (Delt) . Many  times  it  is  easier  to  edit  a sensor  in  the  Zocm 
Overlay  because  it  is  blcwn  up  in  such  a way  that  the  bad  data  are  easily  seen 
and  accurately  deleted. 

2A-9  Overlay  (7.0) 

Delt  is  the  overlay  vrfiich  edits  the  data.  Seven  light  buttons  are  put 
up  on  the  right  hand  side  of  the  screen.  From  this  overlay,  the  user  can  go  to 


93 


I E)atn,  Stor,  Delt,  Hard,  End,  Draw,  or  Eras.  The  overlays  Quad  and  Draw  or  f 

! Zoon-Draw  is  always  activated  first.  Fran  one  to  four  plots  are  already 

(xi  the  screen  at  this  time.  In  the  Delt  Overlay,  the  names  of  the  sensors  on 
! the  side  of  the  plots  are  light  buttons  and  can  be  activated  by  the  light  pen. 

The  first  instruction  is  "Pick  Quad".  The  user  must  activate  by  the  light  pen 
the  sensor  he  is  going  to  edit.  Ihe  next  instructicgi  is  "Move  Cross  in 
(Xiadrant".  Ihe  following  instruction  is  "First  Point".  Ihis  means  that  the 

user  by  means  of  the  light  pen  moves  the  cross  to  the  point  on  tlie  screen  j 

I 

v^iere  he  wishes  to  begin  editing.  This  point  will  not  be  transferred  into 
internal  storage  in  the  program  until  the  instruction  "First  Point"  is  acti- 
vated the  light  pen.  Immediately  the  next  instruction  "Last  Point"  will  I 

flash  CXI  the  screen.  The  user  now  moves  the  cross  to  the  point  in  the  data 
vAiere  he  wishes  to  end  his  editing.  This  point  is  transferred  into  storctge 

after  the  instrucrtion  "Last  Point"  is  hit  with  the  light  pen.  | 

The  user  may  not  move  the  Icist  point  to  the  left  of  the  first  point. 

If  done,  an  instructicxi  will  flash  on  the  screen  "Pick  Point  Again".  The 

user  must  nove  the  cross  to  picdc  a good  secxmd  point  and  hit  the  instructicxi  I 

"Pick  Point  Again"  with  the  light  pen.  j 

dice  the  span  of  editing  has  been  selected,  the  user  c»n  decide  viiich  | 

type  of  editing  he  wishes  to  use.  The  instruction  will  read  "Which  Opticxi"  f 

' t 

follcMad  by  "Delt",  "Subs",  and  "Intp".  Each  of  these  opticxis  is  a li^t 

f 

button,  activated  ty  the  light  pen.  The  "Subs"  cspticxi  has  a further  instruc- 

• K- 

ticxi  "Type  in  XXX. X"  To  Specify  Sub  Value".  When  the  specific  Vcilue  is  typed  [■' 

in  hy  the  user,  that  value  is  plotted  in  the  ^lan  specified  fcx:  editing.  The  f 

I 

"Delt"  option  is  used  to  delete  erronecxis  data  (bad  instrument) . The  "Delt"  1 

option  puts  9999  or  missing  data  into  the  edit  span.  The  9999.  vedue  goes 
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into  the  data  array  but  sane  default  value  is  put  on  the  screen.  Ihe  flag  is 
either  a value  above  the  high  value  for  the  sensor  or  below  the  low  value. 

It  is  easily  distinguishable  fran  true  data.  The  "Delt"  option  must  be 
before  the  "Intp"  option  can  be  utilized. 

The  "Intp"  option  works  on  missing  data.  Interpolation  is  a routine 
vAiich  develops  a "best  fit  curve"  to  estimate  values  for  missing  data  points 
by  using  good  data  points  on  either  side  of  the  missing  data.  If  not  enough 
good  data  points  (at  least  10  on  either  side)  are  present  to  use  in  the  inter- 
polation routine,  the  user  will  be  given  instructions  to  expand  the  range  by 
picking  his  first  and  last  points  again.  The  interpolation  routine  is  lvigp»H 
on  a systems  routine  "SPIU  FIT"  vrfiich  was  adapted  to  our  graphics  routine. 

In  all  of  these  options,  after  our  editing  span  and  our  option  has 
been  selected,  the  screen  is  conpletely  ercised  of  all  plots.  In  a few 
seconds,  the  plots  reappear  with  the  corrections  added  to  our  selected  plot. 
Delt  can  be  activated  over  and  over  again  to  edit  the  data  as  much  eis  nec- 
essary. If  the  editing  is  done  in  the  Zoom  Overlay,  it  will  be  retained  and 
shown  when  four  plots  are  put  on  the  screen. 


I 

2A-10  Overlay  (8,0) 

Hcurd  is  the  overlay  that  puts  the  data  that  are  plotted  on  the  screen  | 

into  the  plotting  queue.  You  can  have  one,  two,  three,  or  four  plots  on  one  ' | 

graph.  The  user  can  put  three  exco  plots  up  on  the  screen  at  the  same  time.  j 

I 

All  the  gra{^  will  have  their  own  scale  depending  on  the  type  of  sensor.  !' 

Clearly  labelled  symbols  will  distinguish  one  sensor  from  another.  Ihe  time 

t 

series  will  be  reeil  time  for  hard  copy.  However,  cai  the  screen,  the  user  will  j 

see  hours  beyond  24  (26,  27) . This  is  just  a convenience  used  for  the  program 


\ 
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1 

i 


f 


and  they  indicate  a change  of  day.  Ihe  time  scale  for  hard  copy  will  only  | 

be  the  first  five  hours  even  though  six  hours  of  data  are  presented  on  the 

screen.  As  we  stated  earlier,  the  sixth  hour  is  a dc^licated  hour.  The  I 

Overlay  Draw  must  be  activated  before  Hard.  (Draw  puts  the  plots  on  the  | 

I 

screen) . Once  Hard  is  activated,  the  sensor  names  on  the  side  of  the  plots 

beocme  light  buttons.  The  first  instruction  is  "Select  Tah^is  or  Go".  This 

label  is  eilso  a light  button.  The  user  then  must  pick  with  the  light  pen  the  i 

sensors  to  be  plotted  and  hit  "Go"  can  the  bottom  of  the  screen.  When  the 

message  at  the  bottom  of  the  screen  disappears,  the  Hard  Overlay  has  been 

corpleted.  The  plot  has  been  put  into  the  plotting  queue  and  will  be  plotted 

at  Central  Site  sane  time  later.  The  user  can  usually  pick  v;p  his  plots  the 

next  day.  1 

i 

2A-11  Overlay  (9,0)  ' 

Star  is  the  overlay  that  chcinges  the  data  array  to  oanply  with  the  ( 

editing  procedinre  acoonplished  during  the  execution  of  the  program.  TVo  data 
arrays  are  read  in  at  the  beginning  of  the  program.  For  storing  puzpx>ses,  one 
file  is  labelled  as  unedited  data;  the  other  eis  edited.  Whenever  ein  editing  | 

routine  (Delt  Overlay)  is  activated,  the  data  that  cure  edited  are  stored  in  an  I 

I 

intermediate  curay.  When  Stor  Overlay  is  activated,  the  data  in  the  inters  j 

mediate  array  are  transferred  to  the  edited  file.  In  the  beginning  of  the  , [ 

program,  both  files  are  identical.  After  the  Stor  routine  is  ooipleted,  the  r 

edited  file  is  edtered.  [ 

• 1“.' 

IT 

2A-12  Overlay  (10,0)  \y 

End  is  an  overlay  that  exits  the  program.  All  grEqphic's  jobs  must  be  f 

completed  by  activating  the  End  overlay.  This  gives  a normal  ending  to  your 
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• 1 1 

: I 

i i 

program  and  assures  that  all  your  editing,  hard  copy  etc.  is  saved.  f ‘ 

2A-13  There  are  a few  limitations  to  graphics  and  the  user  can  easily  1 

ocrpensate  for  these.  It  is  better  to  keep  the  time  on  the  initial  oraitrol 
card  set  relatively  low  (25  sec)  than  to  try  to  do  all  your  editing  in  one 

I 

exercise  of  the  program.  Graphics  takes  a large  amount  of  core  and  supersedes 
anything  else  on  the  oarputer  system.  Since  grajrfu.cs  has  top  priority  in  the 
system,  once  you  put  your  program  in  the  ccmputer,  it  almost  imtediabely  ocmes 
ip  on  the  screen.  However,  if  the  system  is  overloaded,  graphics  is  the  first 

to  be  diirped.  If  graphics  is  not  exited  the  normal  way  (End  Overlay)  but  is  j 

either  pulled  off  the  systan  by  an  operator  or  a "com  lost",  all  the  editing  ! 

or  hard  ceding  the  user  has  done  is  irretrievably  lost.  It  is  better  to  | 

load  the  program  several  times  and  do  fairly  small  sections  at  one  time. 

I 

t 

I 


i 

i 


I 


r 

1-' 
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2A-14  We  origincil  list  set  i?)  for  grajAiics  is  as  follows 


1. 

FSM  A15 

24. 

SPD  AlO 

47. 

PYRH  I 

2. 

FSM  AlO 

25. 

DIR  AlO 

48. 

PYFH  II 

3. 

FSM  AOS 

26. 

SPD  AOS 

49. 

BKG 

4. 

FSM  AOl 

27. 

DIR  AOS 

50. 

RAIN 

5. 

FSM  BIO 

28. 

SPD  XOl 

6. 

FSM  BOS 

29. 

DIR  XOl 

7. 

FSM  BOl 

30. 

W1P  A15 

8. 

FSM  COl 

31. 

EffT  AlS 

9. 

FSM  PIO 

32. 

IMP  AlO 

10. 

FSM  P05 

33. 

DPT  AlS 

11. 

FSM  POl 

34. 

IMP  AOS 

12. 

FSl  QIO 

35. 

DPT  AOS 

13. 

FSM  QOS 

36. 

IMP  XOl 

14. 

Fa4  QOl 

37. 

DPT  XOl 

15. 

FSM  XOl 

38. 

TRAN  AB 

16. 

FSM  YOl 

39. 

TRAN  BC 

17. 

FSM  ZOl 

40. 

TRAN  KL 

18. 

SPD  PIO 

41. 

TRAN 

19. 

DIR  PIO 

42. 

NEPH 

20. 

SPD  P05 

43. 

VID 

21. 

DIR  P05 

44. 

PYRG 

22. 

SPD  A15 

45. 

PYR  AI 

23. 

DIR  A15 

46. 

PYR  All 
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3.0 


SATELLITE  STUDIES 


3 . 1 Introduction 

The  satellite  program  is  a fair  ranging  one  that  extends  from  the  operation 
and  si^port  of  interactive  computer  hardvrare  to  advanced  concepts  for  the 
ajplications  of  satellite  imagery  data  in  operational  forecasting.  This  chapter 
will  address  the  work  done  on  the  McIDAS  (^to-ccmputer  interactive  tota  Aocess 
Systan)  and  the  use  to  vhich  it  is  being  put  in  the  research  investigations. 
Research  work  has  been  done  in  the  following  areas:  the  physiccil  properties  of 
clouds;  the  objective  specification  of  cloud  type  from  very  high  resolution 
satellite  imagery  by  advanced  statistical  techniques;  the  objective  extraction 
of  satellite  sensed  vertical  soundings;  and  the  collection,  pirooessing  and  analy- 
sis of  imagery  data  for  investigation  of  automated  short-range  forecasting 
techniques . 

3.2  McIDAS  Support 
3.2.1  McIDAS  operations 

McIDAS  was  installed  at  AFGL  in  Jime  1974  by  the  University  of 
Wisconsin.  This  system,  vhich  blends  the  precision  of  a digital  computer  with 
the  convenience  of  a TV  display  and  the-  ingenuity  and  imagination  of  the 
researcher,  provides  a very  powerful  tool  for  the  management,  processing  and 
analysis  of  satellite  data  (Figure  3.1) . In  1975,  a real  tiite  data  acquisition 
capability  was  added  to  the  systan.  The  24-ft  antenna  with  its  cissociated  1687 
mhz  receiver  and  bit  synchronizer,  danodulator,  and  interface  to  the  McIDAS 
oorputer  allowed  the  reception  of  the  visible  and  IR  data  from  too  GOES  earth 
synchronous  satellites,  located  over  the  equator  at  135°E  and  75%.  Coverage 
extended  from  l^jalein  Island  in  the  Pacific  to  Ascension  Island  in  the  Atlantic. 
A 12-db  signal  to  noise  ratio  was  eu:hieved  in  this  euxjuisiticsi  systan  with  a 
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MCIDAS  operator's  console. 


mean  error  rate  of  less  than  1 per  100,000  characters. 

The  software  was  then  i;5x3ated  to  give  a real-time  system  capability.  The 
McIDAS  capability  to  ingest  and  process  satellite  data  includes  the  fitting  of 
orbital  parameters  by  known  lantinarks  to  achieve  a high  degree  of  positional 
accuracy.  Tests  shewed  that  the  positional  error  could  be  reduced  to  close  to 
i 1 picture  elanent  (pixel)  in  the  acquired  image.  The  average  error  measured  on 
235  sets  of  pixels  with  4-mile  resolution  was  4.7  miles.  This  navigational 
capability  permits  the  location  of  the  data  in  areas  of  coiplete  cloud  cover  or 
over  landmark  free  oceans. 

The  hardware  v^ch  makes  up  the  McIDAS  facility  consists  of  the  data 

acquisition  system  previously  mentioned,  a Harris  6024/5  corputer  with  65K  of 

memory  (24-bit  word  length) , 5 disk  drives  with  a total  of  130  megabytes  storage  j 

capacity,  2 tape  drives  for  9-track  800  bpi  magnetic  tape,  line  printer,  card  I 

recorder,  and  TV  console  for  eperator  interaction  with  the  data.  In  December  | 

1977  a second  TV  terminal  was  added  hy  the  University  of  Wisconsin.  This  update  j 

also  provided  a Write-Randem  Raster-Read  Manory  (WRRRM)  to  generate  graphics  and  ' 

plotting  presentations  overlayed  or  independent  of  the  satellite  images,  as  | 

desired.  Fol  owing  this  hardware  and  software  modification,  data  from  the  FAA 

high  speed  weather  circuit  were  prograititied  into  McIDAS  so  that  surface  and  upper  | 

( 

air  data  for  North  J'merica  could  be  stored  and  analyzed.  The  c^aerating  system 

used  by  McIDAS  is  a Harris  Disk  Monitor  System  modified  by  the  University  of  f 

Wisconsin. 

I 

IXiring  the  period  when  support  was  given  to  operations  in  the  western  y 

Pacific  Ocean,  data  were  required  frem  the  Pacific  GCES  satellite  for  the  Kwajalein 
area,  and  in  addition,  the  short-range  forecasting  program  needed  data  from  the 
Atlantic  satellite  for  the  New  England  area.  This  could  be  acoerplished  by 
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switching  the  antenna  fran  one  satellite  to  the  other  during  their  transnission  * • 

periods.  The  eastern  GOES-2  trananits  fron  00  to  18  minutes  and  from  30  to  48 
minutes  after  every  hour.  As  long  as  data  are  not  required  south  of  40°S  latitude, 
both  transnissions  can  be  received.  Manual  r^ositioning  of  the  antenna  requires 
about  three  minutes  and  curtails  the  extent  of  the  data  vrfiich  can  be  collected. 

Dual  antenna  cotmand  units  were  installed  to  give  push-button  selection  of 
antenna  positicxi.  This  change  redxiced  the  repositioning  time  fran  3 minutes  to 
40  seconds  and  has  eliminated  the  operational  prt±>lejn  of  supporting  missions  with 
data  fran  t»o  satellites  simultaneously. 

3.2.2  McIDAS  software 

The  McIDAS  software  is  oorposed  of  two  parts.  First  is  the  system 
software  that  consists  of  data  files,  control  of  the  video  display,  data  ingest  J 

and  control  of  data  applications  software  program  files.  The  second  is  the  i 

data  applications  software,  vrfiich  includes  cloud  altitude  algorithns,  winds  frxan  j 

cloud  motions,  pressure  field  displays,  tenperatiare  displays,  gridding  of  j 

satellite  imagery,  curve-fitting  and  statistical  analysis. 

Upon  ootpletion  of  the  system  reoonfigijration  in  December  1977,  99% 
of  the  system  software  was  conpatible  but  approximately  30%  of  applications  | 

software  was  not  caipatihle  or  possessed  logic  faults  vAien  used  in  APGL's  MsIDAS  ■ 

system.  Included  in  the  December  1977  reconfigiaration  were  two  additional  | 

^plications  software  categories.  The  first  category  is  the  NOWCASTING  group  j; 

v4iich  displays  ocxiventional  weather  observations.  Both  hourly  surface  r^xarts  ^ 

(100  hours  of  more  storage  is  available  for  over  1100  stations)  and  vpper-cdr 
rcdiosonde  data  can  be  analyzed  in  a variety  of  forms.  Hie  second  category  of 
additicmal  ^plications  software  is  the  GRAPHICS  groip,  vdiich  controls  and  i 

processes  data  for  a plot  display  file  or  for  the  electronic  WRRM  plotter/analyser. 
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As  of  October  1978,  all  applications  software  was  ccnpatible  with  the 


current  systan  configuration  although  there  are  sane  logic  faults  in  the 
NCWCASTING  software.  These  faults,  however,  do  not  severely  limit  the  display 
capability  of  the  FAA  line  data. 

In  order  to  correct  applications  software  logic  faults,  as  well  as  to 
trouble-shoot  hardware  problems  effectively  and  efficiently,  an  understanding 
of  the  McIDAS  operating  software  is  essential.  There  are  two  major  software 
modules  that  control  McIDAS'  interactive  capability.  The  first  is  named  RQ-ICCM 
and  is  responsible  for  handling  all  asynchronous  inputs  to  the  systan.  It 
handles  ingestion  of  satellite  imagery  and  teletype  data.  Also,  it  processes 
all  inputs  fran  the  interactive  terminals  as  well  as  inputs  from  the  Acronym 
terminal.  The  second  module,  known  as  TVCTRL,  controls  the  entire  video 
display  systan.  It  recognizes  requests  fran  application  software  and  proceeds 
to  effect  the  display  that  is  desired.  These  two  major  software  modules  are 
new  canpletely  docunented.  There  still  remain  a nianber  of  sub-modules  that 
require  docunentation.  The  documentation  has  taken  the  form  of  caments  in 

the  source  input.  Virtually  every  assanbler  statanent  is  associated  with  a : 

oaiiTfent  to  explain  its  role  in  the  entire  ensanble  of  the  software  module.  ! 

! 

Thus,  if  an  individucil  examines  the  documented  modules,  he  can  understand  ; 

cxxnpletely  the  software  modiile's  function  by  reading  the  detailed  carments.  ' 

' i. 

Table  3.1  lists  the  major  operating  systan  modules  and  sub-modules  that  have  f- 

been  docunented  thus  far.  Table  3.2  lists  most  modules  required  to  have  a ^ 

functional  McIDAS  systan.  It  is  easily  seen  that  much  documentation  work 
remains,  ^^lications  software  were  not  listed  in  these  tables. 

Efforts  have  also  been  devoted  to  developing  new  software  to  support 
the  research  programs.  A McIDAS  prototype  version  of  a cloud-motion  vector 
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algorithm  for  eventual  use  on  a Icirge  scale  computer  hcis  been  written.  It 
provides  printout  of  real-time  analysis  of  cloud  motions  from  consecutive 
satellite  images.  Another  new  software  algorithm  will  yield  probability 
estimates  of  visibility  from  current  observations.  The  initial  version  com- 
putes visibility  probabilities  for  Otis  AFB  and  Scott  AFB.  Plans  are  presently 
being  formoLated  to  compute  prctoability  estimates  for  all  military  bases  in 
the  continental  United  States.  Tlvey  will  be  displayed  in  real  time  on  a 
Mercator  projection  map  for  user  interpretation  and  analysis.  Provisions  have 
been  made  in  this  software  to  allow  manipulation  of  input  parameters  to  study 
the  effects  on  the  final  probability  estimates. 

Other  McII3AS  applications  software,  v^ch  has  been  designed  and  imple- 
mented, relate  to  collecting  imagery  and  teletype  data.  In  all,  two  softwcure 
modules  have  been  written  for  them.  The  first  allows  the  main  frame  conputer 
to  schedule  ingestion  of  GCES  imagery  through  the  Off-line  Data  Ingest  System 
(GDIS) . This  irr^ule  is  necessary  because  the  ODIS  microcatputer  cannot  recog- 
nize any  time  beyond  Greenwich  midnight.  Since  the  McIDAS  main  frame  computer 
does  not  have  this  difficulty,  it  can  set  up  a schedule  for  imagery  ingestion 
into  the  ODIS,  without  any  hunan  intervention.  This  is  a benefit  particularly 
if  data  collection  is  required  during  off-duty  hours.  The  second  applications 
software  module  designed  to  collect  data  is  concerned  with  hourly  teletype  data. 
Until  September  1978  this  was  all  done  hy  hand  from  hard  copy.  McIDAS  now 
performs  this  chore  tpon  oonmand.  Future  plans  include  the  automatic  editing 
of  conventional  ctoservations  in  real-time.  This  will  make  the  latest  accurate 
observations  continuously  available  for  analysis  and  display  by  McIDAS. 
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TABLE  3.1 


DOCUMENTED  SYSTEM  SOFTWARE 

REMOTE  OOMMUNICATICN  CONTTROL 
ASYNC  Asyncronous  interrupt  handler 

CATJC  FAA  ingest  file  manager 

CCMM  Afplication  software  message  handler  (receiving) 

CRBH  Card  reader  handler  for  EOT  to  UW 

OQUE  i^lication  software  message  handler  (sending) 

KBC  Keyboard  handler  (all  terminals) 

ND  & NDAT  McIDAS  cormon  output  handler 

NUMC  FAA  line  current  record  address  handler 

REMPSA  Nfein  loop  of  McIDAS  remote  camunication  control  j 

RTAC  FAA  line  ingestion  system  initializer 

SYSQ  Main  applications  software  initiator  loop  modiile  , 

SVCART  Input  byte  from  FAA  line  | 

j 

TTOH  Input  byte  handler  for  terminals 

WD  & WDAT  McIDAS  odimon  input  handler  ' 

I 

XMIT  Special  CdM  handler  for  KBC  t 

i 

i 

! 

I 

[' 

>■ 

I 
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VIDEO 

AHEAD 

CURJOY 

CURSIZ 

FIX 

FHMBAK 

lADCW 

IVSYS 

KURPOS 

KURSIZ 

LTV 

SCHECK 

TPMIL 

TVCPSA 

WRMIO 


TABI£  3.1  (OONT'D) 

DOOJMEOTED  SYSTEM  SQFTVJARE 

DISPIAY  CCNTROI. 

Step  forward  one  frame  on  video  disk 

Determines  cursor  position  as  TV  line  and  elenent  frcm  McIDAS  ocmnon 
Dei>vtmines  cursor  size  as  TV  line  and  elonent  frcm  McIDAS  ocnmon 
Build  video  canmand  half  words  frcm  two  full  words 
Step  backward  one  frame  on  video  disk 
Builds  video  disk  medn  control  word 

Determines  video  system  for  termincd  j 

Hardware  control  word  processor  for  cursor  position  j 

Hardvrare  control  bits  processor  for  cursor  size  word 
TV  space  to  time  mapping  routine 

Auto-swe^  checking  routine  1 

Determines  requesting  terminal  ' 

Main  Video  controlling  routine  j 

WRRRM  outqxit  routine 

f 

I 

t- 
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TABIE  3.2 


McIDftS  System  Software 


REMOTE  OCMUNICATION  (OTNTRDL 


ASYNC 


AsynchrcaTous  intern^  hardier 


CALLC 


FAA  ingest  file  manager 


i^lications  software  message  hardier  (receiving) 


ISBIAN 


IPCfiE 


ND  & NDAT 


REMPSA 


RTSVDC 


SVCART 


Card  reader  "B"  hardier  for  ranote  jc±»  entry  to/fron  UW 
^^lications  software  message  handler  (sending) 

CXitput  line  righthand  deblanker  for  CRT  display 
Fax  data  ingestion  handler 
Arconym  temincil  interrupt  subsystem 
Keyboard  hardier  (all  terminals) 

Operator  Ccrmunication  handler 

Line  prints  "E"  hardier  for  ranote  jcto  entry  to/fron  UW 

MoIDAS  common  outfit  routine 

FAA  line  current  record  address  routine 

Main  loop  of  McIDAS  ranote  camtunication  control 

FAA  line  ingestion  system  initializer  routine 

Main  FAA  line  ingest  loop 

Main  afplications  software  initiater  loop  module 

Input  byte  fron  FAA  line 

Terminal  "A"  intempt  subsystem 

Terminal  "B"  interrupt  subsystem 

Terminal  "C"  interrupt  subsystem 

Irput  byte  handler  for  video  termincLls 
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I 

I 

1 

■m 

TWAITS 
WD  & WDAT 

‘ XMTT 

VIEEP  DISPLAY 

AHEAD 

CURJOy 

CURSIZ 

FIX 

FFMBAK 

FPMSYN 

lADCW 

IVSYS 

IWRM 

IWFMEO 

KURPOS 

KURSIZ 

LOOP 

LTV 

OVRIAy 

t 

SCHECK 
SWEEP 
TCT 
TFMIL 
TVCPSA 
VELOUR 


TABIE  3.2  (CONT'D) 

Key-in  processor  for  applications  software 
Software  executicxi  suspension  routine 
McIDAS  oatinon  input  handler 
Application  software  initiator 

SYSIFM 

St^)s  forward  one  frame  on  video  disk 

Determines  cursor  position  as  TV  line  and  element  from  McIDAS  ootiton 

Determines  cursor  size  eis  TV  line  and  element  fran  McIDAS  cannon 

Builds  video  cotmand  half  words  from  tsio  full  words 

Step  back  one  frame  on  video  disk 

Delays  1/30  sec  for  frame  synchronization 

Builds  video  disk  main  control  word  j 

Determines  video  system  for  requesting  terminal  , 

Builds  WRRRM  control  word  j 

Controls  all  output  to  WRRRM 

Hardware  control  word  processor  for  cursor  position 
Hcirdweure  coitrol  word  processor  for  cursor  size 

Auto-loop  frate  routine  ' 

TV  space  to  time  mapping  routine  i 

Overlays  two  video  output  signads  r 

Check  for  auto-loop  time-out  j 

Controls  auto-loop  B 

Obtains  and  prepares  all  video  control  words  > 

Determines  requesting  terminal  f 

I 

Efeun  video  controlling  routine 
Controls  velocity  cursor 
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3.3  Vertical  Soxmding  Evaluation 

Infrared  sensors  capable  of  detecting  energy  over  snail  spectral 
intervals  in  the  wings  of  the  IS-ndcrcmeter  carbon  dioxide  absorption  band, 
can  be  conbined  into  a system  to  sense  the  vertical  tarperature  distribution 
in  the  atmosphere.  The  principle  is  based  on  the  fact  the  amount  of  energy 
anitted  by  the  atmosphere  in  narrcw  spectral  intervals  on  the  wings  of  the 
IS-microneter  CD2  absorption  band  is  related  to  the  vertical  distribution  of 
terperature  in  the  atmosphere.  Different  spectral  intervals  anit  itost  of 
their  energy  fron  different  altitude  layers  of  the  atmosphere.  Figure  3.2 
illustrates  the  weighting  functions  used  to  ocrpute  the  energy  reaching  space 
fron  a selected  set  of  spectral  intervals  (channels) . 

Ihe  problem  facing  the  satellite  meteorologist,  however,  is  not  how 
much  energy  is  reaching  space.  That  quantity  is  measured  at  the  satellite. 
Ihe  meteorologist  needs  to  knew  vhat  taiperature  distribution  created  the 
amounts  of  energy  sensed  by  the  various  channels.  This  is  analogous  to 
asking  that  given  $2.37  in  change,  how  many  pennies,  nickels,  dimes,  etc. 
are  there  in  my  pocket?  Ihere  are  any  nimber  of  answers  to  the  radiation 
problon  al.  '>.  A catplex  mathematical  method  called  the  inversion  technique 
can  be  used,  in  theory,  to  recover  the  vertical  temperature  profile  fron  the 
radiances.  Statistical  methods,  building  on  past  data,  have  also  been 
developed.  All  preposed  solutions  for  this  prc*>lem  of  temperatxire  retrieval 
are  based  on  assinptions  which  are  more  valid  at  one  time  or  place  than  at 
another.  The  selection  of  the  best  technique  is  therefore  dependent  on  a 
ociiparison  with  direct  observations. 

Such  a oerparison  would  seem  to  be  ooiparatively  simple  and  straight- 


I 
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forward  except  for  one  further  ocxnplication  vhich  arises  fron  the  unoertainty 
of  the  endssivity  of  the  ground  at  any  given  time  or  place.  That  ccrplication 
is  that  sounding  data  are  valid  cxily  over  water  vhere  the  surface  endssivity  is 
accurately  known.  Evaluaticn  of  sounding  data  have  been  made  over  the  oceans 
against  islard  and  occasional  ship  reports.  These  occasional  r^xDrts  were 
usually  separated  in  time  and/or  space  from  the  satellite  r^»rt.  In  February 
1975,  the  Glc±al  AtmosEheric  Research  Program  conducted  the  Air  Meiss  Transforma- 
tion E>q)eriment  (AMIEX)  in  the  P2x:ific  Ocean  south  of  Japan.  Five  weather  ships 
were  stationed  in  a snail  curea  cind  this  presented  a unique  oppoirtunity  to  make 
direct  cotrparisons  between  radiosonde  soundings  and  CMSP  (Defense  Meteorological 
Satellite  Program)  satellite  data.  Figxare  3.3  shows  the  area,  locaticm  of 
radiosonde  stations,  ground  track  of  a typical  satellite  pass  over  the  area,  and 
the  locations  of  the  soundings  along  one  scan  of  the  satellite  system. 

IXjring  the  period  Pdd  14  through  FdD  28,  1975,  there  were  50  CMSP 
soundings  that  were  ^spropriate  for  evaluation.  The  ship  data,  however,  were 
not  of  the  best  quality  and  there  vras  doubt  about  the  value  of  direct  oonparison 
between  the  two  systems.  Accordingly,  in  place  of  the  discrete  ship  radio- 
sondes, ta  eratures  extracted  for  the  c^^jropriate  levels  from  U.S.  Navy  upper- 
air  analyses  of  the  Pacific  ware  used.  These  data  are  called  IT  PAR  and 
contain  a oert2un  amount  of  horizontcil  and  vertical  smoothing.  This  is 
probably  not  detrimental  to  the  ocnparison  since  horizcntal  smoothing  is 
inherent  in  the  observed  satellite  data  because  the  energy  is  integrated  over 
a 33-km  diameter  circle  at  the  satellite's  subpoint  and  verticed.  smoothing 
is  inherent  in  the  conversion  of  radiance  to  sounding. 

The  50  satellite  daservations  were  taken  from  2 CMSP  satellites,  one 
with  ascending  and  descending  nodes  near  dawn  and  dusk,  and  the  other  with 
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nodes  near  noc3n  and  midnight.  Satellite  derived  temperatures  within  ^ 4 hours 
of  OOZ  and  12Z  were  cxxfpctred  with  the  corresponding  upper-air  temperatures. 

Tables  3.3  through  3.7  are  histograms  that  show,  by  pressure  level  in 
the  atmosphere,  the  frequency  that  the  difference  between  the  EMSP  and  TT  PAK 
occurred  for  the  nuiber  of  Celsius  degrees  indicated  along  the  top  of  the  table. 

The  data  are  separated  by  satellite  in  Tables  3.3  and  3.4.  Tables  3.5  and  3.6 

list  the  distributions  obtcdned  from  day  and  night  data  respectively.  All  the  j 

I 

data  are  ocmbined  in  Table  3.7. 

I 

This  information  is  svnmarized  in  Figiares  3.4  and  3.5  v*iich  show  the 
vertical  distributic»is  of  mean  differences  between  the  DMSP  and  TT  PAK  [ 

taiperatures  by  time,  satellite  and  the  entire  data  sample.  These  results  are 
very  disoourciging.  They  shew  largest  differences  near  the  surface  and  at  the  j 

I 

trepopause.  Figure  3.4  shews  the  inconsistency  between  satellites  and  Figure  ' 

3.5,  that  between  day  and  night.  The  standard  deviation  very  seldcm  falls  | 

below  2Pc  and  averages  over  3*^.  Even  in  the  middle  trcposEiiere  at  500  nto,  I 

not  only  is  the  range  between  the  various  average  values  large,  the  standard 
deviation  is  also  large. 

I 

In  & Tmary,  ciside  fromi  the  most  general  description  of  the  atraosE*ere, 

I 

these  satellite  sounding  temperatures  are  unreliable  to  the  point  of  useless-  ! 

ness.  The  study  did  serve  its  purpose,  however,  in  bringing  to  the  attention  ^ ^ 

of  the  Air  Weather  Service  (W4S) , the  need  to  reevaluate  the  v^le  sounder  ^ 

program  fran  sensors  to  the  uses  to  vhich  such  data  could  be  reliably  applied. 

I, 

•jr 


113 


Table  3.3  Frequency, mean  and  standard  deviation  of  ten|)erature  differences  at  stcuidard  levels 
satellite  No.  8531  amd  oorre^xxiding  IT  PAK  data.  (22  soundings) 


Table  3.4  Seme  as  Table  3.3  exo^Jt  for  CMSP  satellite  No.  9532.  (28 


3.5  Sene  as  Table  3.3  except  for  OOZ  (day)  observations  only  of  both  satellites.  (31 


T^le  3.6  Same  as  Table  3.3  except  for  12Z  (night)  observations  of  both  satellites.  (19  soundings) 


I^le  3.7  Sene  as  Table  3.3  except  all  observations  frem  both  satellites.  (50 


Fig  3.5  Same  as  figure  3.4  except  dashed  line  OOZ  (day)  data;  dot-dashed  line  is 
12Z  (night)  data  and  solid  line  is  all  data. 


3-4  Autonated  Cloud  Identification  by  Second  Order  Statistical  Methjds 
3.4.1  Introduction 

Meteorology  is  characterized  by  enormous  quantities  of  conventional 
data  continuously  pouring  into  weather  centrals  fran  all  over  the  world.  Tbe 
addition  of  satellite  data  has  augmented  this  data  glut  orders  of  magnitude. 
For  example,  in  one  day  10^^  data  bits  are  collected  by  only  one  Geosynchronous 
Earth  Orbiting  Satellite  (GOES) . This  high  volume  of  data  has  a low  inforna- 
tion  density.  Methods  are  needed  to  reduce  these  data  into  a set  with  high 
information  density. 

Many  cor^juter  systems  aid  algorithms  have  been  developed  to  assist  in 
the  efficient  extraction  of  information  from  raw  data.  Interactive  systems 
such  as  the  Man-Ooiputer  Interacti  ve  Data  Access  System  (McIDAS) , have  been 
developed  to  combine  the  electronic  speed  of  a ccnputer  with  the  decision 
capability  of  a skilled  analyst.  At  Air  Force  Global  vgieather  Center  (APGWC) 
a ccrputer  data  reduction  algorithm  has  been  developed.  The  algorithm 
produces  a global  cloud  analysis,  called  3-D  Nephanalysis  (3D  NEPH) , that  is 
routinely  generated  for  15  levels  in  the  atmosphere  and  for  25-mile  square 
boxes.  The  cloud  field  is  defined  as  clearly  as  possible  using  surface, 
radiosonde,  aircraft,  and  satellite  data  as  available.  Observational  modes 
are  assigned  priority  in  the  corputations  by  means  of  a decision  tree  process. 
Since  global  satellite  data  is  routinely  available,  it  is  possible  to  generate 
a 3-D  NEPH  from  satellite  data  alone.  Ffcwever,  because  of  ocirputational 
limitations,  the  data  that  are  used  have  been  degraded  in  resolution  to  reduce 
the  quantity  of  the  irput.  Consequently,  there  is  a reduction  in  the  quality 
of  the  output.  APCJC's  operational  ccninitments  require  the  development  of  a 
new  satellite  processor,  capable  of  handling  increased  amounts  of  satellite 
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data  with  a higher  level  of  sc^Aiistication  and  greater  efficiency.  This 
section  briefly  reviews  the  pressed  method  in  the  development  of  the  new 
satellite  processor  for  AFX3WC. 


3.4.2  Nature  of  the  problem 

Despite  APCSC's  leadership  in  the  applicaticm  of  satellite  data, 
requiran^ents  placed  on  them  by  their  custoners  exceed  their  processing 
capability.  Smoothed  mode  satellite  data  at  ^jproodmately  3-nn  resolution  are 
presently  used  as  input  to  3-D  NEPH.  Although  fine  mode  data  at  1/3  rm 
resolution  are  sensed  by  the  satellite,  they  have  not  been  used  until  recently. 
APGWC  cenputers  are  inadequate  for  processing  satellite  data  in  real  time 
with  the  rates  and  volume  of  data  involved.  The  fine  mode  sensor  of  the  DMSP 
Block  5-D  system  generates  over  100  times  as  much  data  as  the  high  resolution 
sensor  in  the  Block  5-C  system. 

Improvement  in  operatiaial  cloud  ^)ecificatic»i  at  AFGWC  implies  the 
utili2sation  of  fine  mode  data.  Thus,  the  new  satellite  processor's  task 
would  be  data  reduction  of  an  enlarged  data  base  ocxiteiining  1/3  rm  imagery 
vAiile  providing  improved  information  content.  Considering  the  quantity  of 
data  involved,  sopiiisticated  programming  alone  cannot  cate  close  to  allevia- 
ting the  prctolan.  Tie  approach  selected  by  APGWC  and  APGL  is  to  utilize 
scfi^sticated  analytical  techniques  that  are  adaptable  to  hard  wired  oatpu- 
tationcil  devices.  Specifically,  utilization  of  Fast  Fourier  Transform  (FFT) 
coefficients  as  indices  of  cloud  types  has  been  pursued. 

3.4.3  The  Fast  Fourier  Transform  (FFT) 

The  Fourier  Transform  in  the  past  has  been  used  for  characterizing 
linecu:  systems  and  for  identifying  frequency  coiponents  making  up  a continuous 
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I waveform.  Hcwever,  v^ien  the  waveform  is  sanpled,  or  the  system  is  to  be 

analyzed  on  a digital  ocmputer,  it  is  the  finite  discrete  version  of  the 
Fourier  Transform  that  must  be  used.  Hie  Feist  Fourier  Transform  (FFT) , as 
developed  by  Cooley  and  Tukey,  is  sinply  an  efficient  method  of  coniJuting  the 
discrete  Fourier  Transform. 

The  discrete  Fourier  Transform  pair  can  be  written  as: 

, N-1 

la)  Y(j)  = 1 z X (k)  exp  (-i2iTik/N) 

N k=0 

N-1 

lb)  X(k)  = Z Y(j)  esq)  (i2TTj)c/N) 
j=0 

where 

Y(j)  is  the  frequency  danain  function 
X(k)  is  the  time  dcmain  function 

i = /:r 

and  N is  '-he  nunber  of  points  sanpled.  The  discrete  Fcxarier  Transform  can 
be  oorputed  fron  equations  la  and  lb  direcrtly,  but  it  is  slower  than  using  the 
FFT.  For  an  N-point  transformaticai,  v^iere  N is  seme  power  of  2,  (N  = 2'*',  m 
being  positive) , the  niiriber  of  ccnplex  additions  and  subtractions  is  given  by 
N^.  For  the  PTT  the  niirter  of  oonplex  additions  and  subtractions  is  reduced 
to  (N/2)  log2N.  For  N = 1024  this  represents  a ocnputational  savings  of  200 
to  1. 

The  output  of  an  FFT  is  sinply  the  coefficients  of  a Fexirier  series 
representation  of  a function.  A one-dimensional  Fourier  series  representation 
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of  a function  reduces  to  a sijiple  hemnnonic  expressed  by; 


2)  F(x)  - Aq  + E (A4  Cos  + B^Sin  ) 

j=l  ^ 

vAiere, 

T = period  of  the  function 
X = satpled  dcroain 
j = jAiase 

and  Ao»  Aj,  and  Bj  are  the  coefficients.  For  a one-dimensional  FPT  of  an 
8-pixel  (picture  element)  scan  line  the  output  consists  of  the  17  coefficients 
of  Aj,  and  Bj,  with  j ranging  frcm  1 to  8.  In  a two-dui^isional  FPT  of  an 
8x8  pixel  area  the  output  would  consist  of  129  coefficients  Aq,  Aj,^  and 
Bj^j^  with  j and  i ranging  from  1 to  8. 

Ihe  developnent  of  a new  satellite  processor  for  APQWC  involves  the 
oaiputation  of  the  power  spectnin  from  two-dimensional  FFT  coefficients.  It 
is  anticipated  that  proper  processing  of  spectra  will  extract  information 
from  the  satellite  imagery  field  vAiile  functioning  as  a da-ha  reduction 
mechanism. 

3.4.4  Ccrputation  of  power  spectra 

Ccrputation  of  power  ^sectra  from  two-dimensional  Fourier  Transforms 
, yields  a new  set  of  parameters  in  response  to  spatial  distribution  and 

frequency  ciiaracteristics  of  the  cloud  field.  If  some  cloud  types  generate 

0 

reproducible  power  ^lectra  parameters,  one  cam  classi:^  clouds  by  using  a 
discriminant  function.  ^Jectrad.  signatures  (reproducible  cloud  spectra) 
indicate  size,  shape  and  texture  of  a particular  cloud  type.  Thus,  power 
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spectra  are  a measure  of  the  amount  of  variance  in  an  imagery  field  due  to 
clouds  of  various  sizes. 


f 


Power  is  defined  as 


3) 


‘’1,1  = ■'''"j-i  * 


B2.  . 


^,i  ^ (cosine  term  coefficient)  and  is  the  imaginary 

part  (Sine  term  coefficient) . Poster  spectra  were  initially  confuted  by 
annxolarly  integrating  the  FFT  coefficients.  Coiputation  of  power  ^jectra  in 
annular  rings  is  permissible  since  the  BTT  retains  the  prxsperty  of  a periodic 
function  allcxidng  dr5)lication  of  FET  output  as  a half-range  expansion.  llTe 
periodicity  allows  drplication  of  the  FFT  output  in  quadrant  I of  the  frequency 
plane,  into  the  II,  III  and  IV  quadrants.  Ihe  dvplication  of  the  FFT  outp»ut 
in  the  ranaining  quadrants  sinplifies  the  annulai'  integration  of  the  power 
spectra. 

Poer  spectra  that  has  been  integrated  annularly  may  be  biased  as  a 
result  of  >'‘*T  aliasing.  Aliasing  is  the  result  of  high-frequency  ocmponents 
inpersonating  low-frequency  ccmponents.  It  is  minimized  by  ocrnputing  the 
half  of  the  FFT  output  by  synmetric  oonplex  ccaijugates.  Coefficients 
with  frequencies  higher  than  the  linear  Nyquist  frequency  (fg/2,  v*>ere  fg  is 
the  sarrple  dcntain)  are  ocnputed  ty  the  symmetric  ocnplex  conjugates  of 
frequencies  Icwer  than  the  linear  Nyquist.  Since  the  spectra  are  integrated 
on  a t>o-dimensional  plane,  the  spectra  are  coiputed  up  to  what  is  termed  the 
tvc-dimensioneil  Nyquist  frequency,  given  by: 
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/(sanple  size  of  first  dimensicgi)  ^ 
4 


(sanple  size  of  secxand  dimension)  ^ 
4 


Tlius,  an  area  of  8x8  pixels  vd.ll  have  a two-dimensional  Nyquist  frequency 
of  6.  With  an  (j,i)  Cartesian  coordinate,  a distance  frcm  the  nean  qaectral 
coefficient  (Pq^q)  can  then  be  ocnputed.  The  ratiaining  terms  (Pj  vd.th  the 
same  distance  from  are  sunned.  Terms  above  the  twD-din«nsional  Nyquist 
(distance  greater  than  6 for  an  8x8  pixel  area)  are  folded  back  into  lower 
distances.  Integration  of  the  FPT  coefficients  in  each  annular  band,  using 
this  technique,  is  done  by  running  through  the  FTT  output  array  only  once. 
Iherefore,  duplication  of  FPT  output  in  Quadrants  II,  III,  IV  is  not  required. 


3.4.5  Normalized  averaged  power  spectra  (NAA)  cotputation 

The  oonputation  that  is  made  is  not  pcwsr  in  the  strictest  definition 
of  the  term,  but  vihat  may  be  called  normalized  averaged  amplitudes  (NAA) . 

The  ooiputations  of  NAA  result  by  dividing  the  sutmed  pcwer  of  eadi  annulcir 
band  by  the  total  number  of  terms  entered  into  the  sun.  The  square  rtxrt: 
of  the  quotient  is  then  taken  and  normalization  is  achieved  by  dividing  the 
averaged  amplitudes  by  the  nunber  of  pixels  entered  into  the  FPT.  NAA  can 
be  expressed  as: 


NAA(f) 


r~\  i j ^ 

Jxl  /otHT  ill  jli  ^j(f),i(f) 


^j(f)  ,i(f)  power  spectra  that  are  d^jendent  xpon  the  distance  they 
are  located  frcm  Pq^q;  Nr(f)  is  the  nunber  of  Pj^fj  teimc  entered  iiico 

the  sun  that  are  the  sane  distance  frcm  Pq^q;  f is  the  distance  fron  Pq^q,* 
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J and  I are  the  dimensional  orders  of  the  imagery  array.  For  8x8  pixels 
where  J = 8 and  1 = 8 normalization  is  achieved  by  dividing  by  the  prodixrt 
of  J X I = 64. 

3.4.6  A study  of  two  cloud  areas 

NAA's  were  ocnpubed  from  the  coefficients  of  FFT's  of  3 m resolution 
satellite  imagery  for  20  October  1976.  TVro  well  separated  areas,  each  of 
which  contained  a roughly  hcmogeneous  cloud  field  over  at  least  a 200-iin  area, 
were  selected  for  the  data  sample.  Identificaticai  of  the  two  cloud  areas  was 
made  solely  frcm  a CMSP  global  mosaic  transparency.  The  first  area  (area  1)  I 

ves  clcissified  as  a mixture  of  stratus  (St)  and  stratocunulus  (Sc)  clouds  and  ^ 

I 

t 

area  2 was  classified  as  a uniform  declc  of  stratus  or  fog.  FPT's  were 
performed  on  each  200-square  nm  area  as  well  as  on  the  64-25  rm  sq  (8x8 

pixels)  area  that  made  vp  the  larger  areas.  j 

There  vere  no  great  differences  between  means  of  NAA's  for  64  j 

(25  X 25  nn)baxes  of  the  two  areas,  other  than  at  frequency  0.  The  difference  I 

in  frequency  0 indicated  that  the  average  brightness  of  area  1 is  greater 
than  that  r ^ area  2.  The  results  were  the  same  vhen  the  200x200  rm  areeis 
were  examined.  A statistic  was  computed  to  determine  vhether  the  differences  ' 

xn  mean  NAA  of  the  two  areas  as  a function  of  frequerxry  were  significant,  ; 

* ^ 

using  each  25x25-rm  box.  The  statistic,  T'  measured  - U2  for  each 
frequency  where  Uj^  represents  the  mean  NAA's  for  area  1 and  U2  represents  the 

I 

mean  NAA's  for  area  2.  At  the  0.01  level  of  significance,  no  significant 
differences  in  mean  NAA  for  frequencies  1,  2,  5 and  6 were  found.  Significant 
differences  for  frequencies  0,  3 and  4 irtply  that  in  one  of  the  areas,  there  ! 

is  a fluctuatioi  in  brightness  between  pixels  that  does  not  exist  in  the  ! 
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other  ctrea. 


Variation  of  NAA  as  a fvinction  of  frequency  for  the  64  (25x25  rm)  boxes 
of  area  1 and  area  2 was  also  ccmputed.  The  variation  of  NAA  for  area  2 is 
vAat  would  be  expected  of  a cloud  deck  with  uniform  tops  that  have  essentially 
identical  brightness  betveen  pixels.  The  higher  variation  for  area  1 indicates 
a cloud  deck  with  non-uniform  tops  and/or  breaks  above  satellite  sensor 
resoluticxi.  A statistic  was  oonnputed  to  determine  if  the  variation  of  NAA's 

o 2 

between  the  two  areas  vras  significant.  If  and  S2  are  the  variances  in 
NAA  for  area  1 and  area  2 respectively,  then  the  statistic  F (F  = Sj^/S2)  could 
be  computed  for  ecich  frequency.  At  the  0.01  level  of  significance,  frequen- 
cies 0,  4,  5 and  6 showed  significant  differences.  The  differences  in 
Vcuriance  for  frequency  0 means  that  sane  of  the  64  (25x25  nm)  boxes  in 
curea  1 had  a greater  variation  in  cloud  brightivess  from  one  pixel  to  another, 
vhile  area  2 did  not  exhibit  this  variable  cloud  brightness.  The  ratios  (F) 
of  greatest  interest  and  inplication  were  for  frequencies  4,  5 and  6.  The 
reason  is  that,  at  higher  frequencies,  the  sampling  rate  of  the  FET  is  at  a 
higher  cycle  than  it  is  at  lower  frequencies.  This  means  that  the  finer 
variable  brightness  in  the  imagery  is  being  sampled.  Thus,  if  one  area  hcis 
a finer  variable  brightness  in  its  imagery  than  another  area  then  the  vari- 
ance of  NAA  cis  a function  of  frequency  will  most  certainly  differ.  This  is 
v^t  was  seen  for  frequencies  4,  5 and  6.  The  finer  texture  may  be  a basis 
for  cloud  t^pe  discrimination,  and  thus  the  differences  in  varianoe  of  NAA 
may  be  a spectral  signature.  The  important  factor  in  this  initial  stucfy  is 
that  enough  differaioe  was  detected  in  FPT  processing  of  the  cloud-field 
imagery  to  justify  further  investigation  of  the  prcblem. 
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3.4.7  Conclusions 


Ihe  results  discussed  in  the  previous  section  (3.4.6)  vgere  based  on 
3-m  resoluticm  data,  primarily  for  test  purposes.  Nevertheless,  there  is  an 
indication  that  processing  of  FFT  coefficients  will  allow  classification  of 
cloud  types  in  an  automated  scenario.  Use  of  the  1/3-nn  resolution  imagery 
would  undoubtedly  enhance  the  high  frequency  distinctions  of  spatial  cloud 
distribution.  This  would  yield  greater  discriminating  power  among  the  class 
condition  covariances  of  the  spectra,  iitproving  reliability  of  the  classifi- 

I 

cation  process.  i 

3.5  Estimates  of  Cloud  Properties  derived  from  Satellite  Sensors 

3.5.1  Visible  and  infrared  radiation  in  cloudy  atmospAieres  I 

A nimber  of  published  studies  has  suggested  that  a physical  basis  j 

I 

exists  for  the  relation  of  cloud  properties  to  satellite  measurements. 

Satellite  sensors  measure  radiation  which  has  been  emitted,  absorbed,  and 

scattered  by  surfaces,  atmospheric  gases,  and  clouds.  These  processes  can  ' 

be  simulated  by  numerical  models.  The  subject  of  radiative  transfer  in  ! 

( 

atmosfAierif'  gases  above  planetary  surfaces  has  been  treated  in  depth  by  | 

Chandrasekhar^.  The  important  case  of  radiative  transfer  in  clouds  ccmposed 
of  sp*ierical  particles  has  been  sunmarized  by  van  de  Uulst^,  vhile  extensive 
tables  of  calculated  results  have  been  published  by  Deirmendjian^.  Both  ' ^ 

authors  extend  the  work  by  Mie^  vhich  mathematically  describes  the  interaction 
of  a hcmogeneous  sp^iere  with  a field  of  electranagnetic  radiation.  t 

The  theory  of  radiative  transfer  in  cloudy  atmospheres  has  been  applied 

f 

to  spjectral  regions  suitable  for  cloud  sensing  from  satellites.  The  | 

principjal  regions  are  the  visible  spjectrxm  from  0.5  to  1.0  iim,  and  the 
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infrared  (IR)  spectnm  from  8 to  12  jjin.  Both  regions  are  suitable  for  cloud 
sensing  since  they  are  "windcw"  regions  for  vrfiich  geis  transndssion  fron  cloud 
to  satellite  is  very  high  and  since  the  radiation  fields  are  strong  due  to 
reflected  sunlight  (during  daytime)  at  0.5  to  1.0  ym  and  terrestrial  emissicxi 
at  8 to  12  van. 

Cloud  particles  absorb  very  little  energy  at  wavelengths  between  0.5 
and  1.0  ym.  Craisequently,  to  relate  cloud  properties  to  radiation  fields 
consistent  with  satellite  vievdng  gecmetry  one  must  primarily  calculate  the 
amount  of  sunlight  scattered  by  cloud  particles  back  toward  the  sensor.  Such 
calculations  have  been  performed  for  ncminal  distributions  of  cloud  particles  i 
in  a nuitoer  of  studies.  IVaney,  Jac(±icwitz,  and  Howell^  found  that  cloud 
albedo,  i.e.,  the  ratio  of  reflected  to  incident  sunlight,  increased  with  the 
total  cloud  thickness  or  mass  of  cloud  particles.  The  relationship  was  not 
linear,  and  for  clouds  with  an  integrated  mass  greater  than  250  g/in^,  very 
little  increase  in  albedo  was  found.  The  maxiimm  albedo  was  about  90%. 

These  calculations  were  for  stratiform,  ncm-precipitating  water  cloixls.  Plass 
and  Kattawar®  made  ccLLculations  for  both  water  and  ice  clouds  at  0.7  ym.  The 
reflected  radiance  was  found  to  be  less  fron  ice  clouds  than  fron  vater  clouds, 
primarily  due  to  less  efficient  scattering  per  unit  cloud  mass  the  larger 
ice  particles.  Multiple  scattering  was  simulated  by  a Monte  Carlo  technique 
for  oarputing  the  exact  three-dimensional  paths  of  scattered  ^totcns.  The 
Monte  Ccurlo  technique  has  been  adapted  to  ntai-stratiform  clouds  by  MoKee  and 

7 

Oox  , among  others.  The  iiqportcuit  result  here  was  that  clouds  of  finite 
horizontal  extent  were  not  so  bright  as  infinite  clouds  vhen  viewed  by  a 
satellite  due  to  a loss  of  scattered  sunlight  out  the  sides  of  the  clouds.  In 
these  studies,  hcwever,  the  rather  arbitrary  cissimpticai  was  made  that  ice 
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cloud  scattering  could  be  approximated  by  a polydispersion  of  ice  spheres 
since  no  general  theory  exists  for  scattering  of  non-spherical  particles  such 
as  ice  crystals.  An  inportant  class  of  ice  particles,  namely,  long  hexagonal 
prians  vMch  are  frequently  c±>served  in  cirrus  clouds,  was  treated  by  Liou® 
by  approximating  the  prians  as  cylinders  and  applying  a scattering  theory 
first  presented  ty  Lord  Rayleigh^. 

In  the  8 to  12  vim  IR  window,  scattering  is  much  less  irrportant  than 
atjsorption  for  modelling  radiative  transfer  in  cloudy  atmospheres  since  water 
and  ice  particles  intercept  and  reemit  substantial  energies  at  these  wave- 
lengths, The  amount  of  incident  sunlight  at  8 to  12  vim  is  negligible,  so 
that  models  are  not  dependent  on  solar  illumination.  On  the  other  hand,  the 
sources  of  levelling  IR  radiation  such  as  surfaces  or  water  vapor  must  be 
known  or  approximated  since  they  contribute  the  energy  vAiich  is  attenuated 
by  clouds  and  also  form  the  backgrounds  against  which  clouds  are  observed  ty 
satellite  IR  radiometers.  IVo  approximations  are  usiaally  made  in  IR  models. 
First,  underlying  svrrfaoes  are  treated  as  blackbody  emitters  so  that  only 
their  tarperatures  and  not  their  reflective  properides  need  to  be  put  into 
the  models.  Second,  many  clouds  are  also  treated  as  blackbody  anitters, 
cil though  this  c^proximation  is  much  more  difficult  to  justify.  In  the  case 
of  thin  ice  clouds  at  high  eiltitudes,  it  is  a very  poor  approximation  as  shown 
by  Kuhn^®.  For  water  clouds  with  particle  mass  densities  on  the  order  of 
0.1  g/m^  or  greater  most  of  the  radiation  measured  by  a satellite  IR  sensor 
originates  within  the  upper  100  m of  the  cloud.  In  other  words,  the  sensor 
"looks"  into  the  cloud,  but  caily  for  short  distances  so  the  blackbody 
approximation  is  acceptable  if  great  precisicwi  is  not  required. 

Models  of  radiative  transfer  have  been  develcped  for  cloudy  atmospheres 
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in  the  8 to  12  pm  region.  Shifrin^^  was  one  of  the  early  investigators  to 
predict  the  strong  absorption  effects  in  clouds  at  these  wavelengths.  Yamamoto, 

Tanaka,  and  Kamitani^^  did  extensive  calculations  of  radiative  transfer  in 
water  clouds.  Calculations  for  ice  clouds  were  later  done  by  Jacobcwitz^^  and 
by  Liou^*^  for  the  important  case  of  semitransparent  or  "greybody"  cirrus 
clouds.  The  calculations  indicated  that  the  upwelling  surface  radiance  would 
be  so  attenuated  as  to  be  undetectable  frcm  satellites  vhen  the  integrated 
v’ertical  cloud  mass  or  cptical  path  was  on  the  order  of  100  g/m^.  Hcwever, 
due  to  the  low  mass  densities  of  ice  clouds  at  cold  tenaperatures , such  a cloud 
might  have  a gecmetrical  thickness  of  5 km. 

In  summary,  considerable  progress  has  been  made  over  the  years  in 

I 

modelling  radiative  transfer  in  the  windcw  regions  vAu.ch  are  used  to  sense  | 

clouds  from  satellites.  These  studies  are  applicable  to  the  forward  problem 
of  modelling,  that  is,  given  certain  cloud  properties  six:h  as  particle  size, 
phase,  shape,  nvxnber  density,  cloud  thickness,  etc.,  predicting  vrtiat  would  be 
measured  frcm  the  satellite.  The  forward  models  clearly  indicate  that 

i 

f 

satellite  measurements  should  vary  with  cloud  prc^ierties  such  as  cloud  thick-  i 

ness  and  meiss  and,  in  particular,  that  dense  clouds  or  thick  clouds  should  be  ^ 

highly  reflective  in  the  visible  window  and  weak  energy  or  cold  sources  in  the 

IR  window.  Moreover,  the  satellite  measurements  would  be  related  to  cloud  t 

amounts  in  a highly  nonlinear  manner  vAiich  depends  on  cloud  physical  proper-  ^ 

r 

ties,  gas  attenuation,  and  the  viewing  gecmetry  frcm  the  satellite.  I 

I 

Even  though  progress  has  been  made  in  forward  models  for  estimating  | 

the  reflectances  and  radiainces  observed  by  satellites  for  given  clouds,  very  j- 

little  parogress  has  been  made  in  inverse  models  for  estimating  the  cloud  ; 

p»roperties  for  given  satellite  measuranents  and  inverse  models  are  the  chief 
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c±)jective  of  the  present  study.  The  forward  models  are  too  simple  to  be 
readily  edited  for  use  as  inverse  models.  First  of  all,  the  radiation 
models  assvme  hanogeneous  clouds  although  in  situ  sanples  fran  aircraft 
indicate  substantial  variations  in  vertical  and  horizontal  dimensions.  The 
radiation  models  assume  that  ice  particles  have  sane  sinple  geonetric  sh^^e, 
vAiich  is  rarely  observed  in  ice  particle  sampling.  Also,  inverse  models 
vrould  require  sane  means  of  estimating  the  cloud  particle  size  distribution, 
v^ch  is  not  directly  measurable  fran  siitellites. 

In  the  present  study,  it  was  not  necessary  to  develop  an  inverse  model 
by  a theoretical  approach  since  measurements  were  available  fran  a lengthy 
experiment  described  by  Conover  and  Bunting^ ^ and  these  measurements  coiiLd 
be  used  to  construct  and  analyze  enpirical  models.  Three  cloud  prcperties- 
Cloud  Mass  (IWC) , thickness  (AH) , and  Environmental  Severity  Index  (ESI)  were 
selected  for  extended  study.  A variety  of  cloud  conditions  was  sampled  from 
research  aircraft  over  mid-latitudes  of  the  United  States  and  analyzed  together 
with  infrared  (IR)  and  visible  (VIS)  nveasuranents  fran  National  Oceanic  and 
Atmospheric  Aininistration  (NQAA)  satellite  sensors. 

3.5.2  Satt  life  measurements 

TTie  satellite  instninents  used  were  the  Scanning  Radioneters  (SR)  of 
the  Improved  TIROS  Operational  Satellite  (ITOS).  These  broadband  sensors 
simultaneously  detected  both  VIS  (0.52  - 0.72  pm)  and  IR  (10.4  - 12.5  pm) 
radiation.  Resolution  at  the  satellite  subpoint  was  approximately  3 km  for 
the  VIS  channel  and  6 km  for  the  IR  channel.  Data  from  consecutive  polar 
orbits  were  recorded  in  the  form  of  a 2048  by  2048  element  array  for  each 
hemispAiere  each  day.  These  arrays  were  aligned  with  the  conventional 
Numerical  Weather  Prediction  grid  (NWP)  so  that  both  IR  and  VIS  were  equally 
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spaced  on  a polar  sterogr^hic  projection.  The  spacing  of  adjacent  grid 
points  on  the  earth  increases  frcm  about  6 km  at  the  equator  to  about  13  km 
at  the  poles.  SR  data  used  in  this  stixJy  were  thiis  m^jped  and  sairpled  to 
provide  a ncminal  10  by  10  km  horizontal  resolution  for  both  IR  and  VIS. 

The  VIS  data  were  normalized  for  solar  zenith  angle  and  inconsistencies 
in  brightness  corrected  at  the  swath  edges  between  satellite  passes  using  a 
bi-directional  reflectance  mcdel  as  applied  Conover  and  Bunting^ Obser- 
vations by  the  IR  sensor  were  less  sensitive  to  viewing  geanetry  and  therefore 
not  normalized  in  a manner  similar  to  the  VIS.  The  IR  data  did  include  a 
correction  for  limb  darkening  due  to  water  v^xar  absorption.  Also,  since  the 
original  IR  data  set  was  recorded  in  counts,  a conversion  bo  tarperature  in 
degrees  Kelvin  was  made: 

For  counts  - 78,  IR  = 164  + counts  (1) 


For  counts 


IR  = 242 


counts  - 78 


At  locations  vrfiere  the  aircraft  recorded  cloud  information,  IR  and  VIS 
data  v^e  averaged  over  about  50  NWP  grid  points  vAiich  represented  an  area  of 
about  70  by  70  km.  This  area  was  increased  in  size  as  a function  of  the 
satellite  zenith  angle.  A maxiimm  area  thus  averaged  was  70  100  km. 

3.5.3  Aircraft  measurements 

At  the  time  of  the  satellite  pass,  an  cdrcraft  equipped  with  cloud 
jAiysics  instrunentation  wculd  descend  from  an  altitude  of  10  km  in  a 35  km 
diameter  spiral  at  a descent  rate  of  0.3  knv'min.  Many  instrvments  were  used, 
including  a "snow  stick"  for  visual  identification  of  crystcil  habit  and  size, 
Johnson-Williams  liquid  vebcr  sensor,  continuous  formvar  replicator,  cloud 
scene  cameras,  and  particle  measuring  spectronebers . An  onboeird  flight 
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director  ctoserved  altitudes  of  cloud  tops  anl  bases,  estinated  cloud  tops  when 
clouds  were  above  the  maxiiman  aircraft  ceiling,  estimated  the  fractional 
coverage  of  cloud  layers,  noted  optical  effects  such  as  halos,  anl  Hcgr-T-j [yvj 
particles  intercepted  by  the  snow  stick.  This  entire  sounding  schedule  was 
ociipletod  in  about  30  minutes. 

Satellite  measurements  are  sensitive  to  clouds  at  all  altitudes  and 
are  recorded  at  a rapid  rate  - less  than  15  seconds  to  cover  the  entire  70  by 
70  km  array.  If  the  aircraft  were  to  loiter  to  <±)tain  a better  estimate  of 
cloud  tcps,  clouds  at  a ].cwer  altitude  might  advect  or  chance  before  the 
aircraft  was  able  to  descend  for  a saitple.  In  all  cases,  the  highest  clouds 
were  sampled  closest  to  the  satellite  pass  time.  Because  of  the  requiranent 
for  the  aircraft  sample  to  be  rapid,  often  the  nixnber  of  cloud  particles 
obtained  was  snail.  In  this  event  other  available  data  such  as  surface  radar 
and  rainfall  reports  were  included  in  an  effort  to  generate  a profile  of 
cloud  mass  density  as  a furiction  of  altitude.  A total  of  41  cases  were  thus 
recorded  and  Table  3.8  is  a suiriary  of  the  data. 

Ihe  aircraft-measured  cloud  properties  were  calculated  fran  the  ccitplete 
profiles  oi  hydroneteor  mass  density  according  to  the  follcwir>g  definitions: 


I«C  = 

/ p dH 

(3) 

> 

ESI  = 

/ pH  dH 

(4) 

AH  = 

/ 6 dH 

(5) 

1 

I- 

Ihe  symbol  p represents  hydroneteor  mass  density,  H represents  altitude,  ^ 
and  6 equals  one  if  p is  greater  than  0.0001  otherwise  6 equals  zero.  ^ 
1W3  is  the  total  clcud  mass,  integrated  vertically,  including  all  ice  and  I 
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TABIE  3.8 


Simultaneous  Satellite  and  Aircraft  Data  Set 

SAIELLITE  DATA  AIFCRAFT  DATA 

CASE 


NO. 

DAIE 

VIS 

LWC 

ESI 

(°K) 

(ft-I/40) 

(g/m2) 

(g  km^/nti^) 

(Ian) 

1 

lLJan74 

281.3 

161.0 

45. 

.35 

1.2 

2 

14Jan74 

259.0 

91.0 

63. 

.18 

.8 

3 

0lFeb74 

263.9 

57.0 

69. 

.65 

.8 

4 

0lMar74 

230.3 

204.0 

2035. 

7.44 

9.0 

5 

03i^r74 

281.7 

40.0 

0. 

0.00 

0.0 

6 

08Apr74 

229.8 

203.0 

661. 

2.79 

8.4 

7 

247^74 

265.3 

184.0 

56. 

0.00 

1.0 

8 

28^^74 

247.5 

203.0 

731. 

2.29 

8.6 

9 

02May74 

234.2 

176.0 

325. 

.77 

2.8 

10 

15Jan74 

283.7 

137.0 

25. 

0.00 

.8 

11 

12Jan75 

223.0 

211.1 

1067. 

3.65 

8.8 

12 

14Jan75 

261.1 

112.8 

486. 

.16 

2.6 

13 

20Jan75 

217.2 

204.8 

1082. 

4.38 

10.4 

14 

22Jan75 

248.4 

141.2 

178. 

1.09 

5.6 

15 

23Jan75 

248.6 

97.8 

71. 

.59 

3.4 

16 

25Jan75 

248.0 

157.9 

474. 

2.24 

8.6 

17 

27Jan75 

262.9 

123.1 

60. 

.01 

1.6 

18 

29Jan75 

257.9 

191.4 

752. 

.23 

4.6 

19 

04Feb75 

251.1 

175.2 

736. 

2.16 

8.6 

20 

llFeb75 

256.5 

128.2 

33. 

.07 

4.0 

21 

13Feb75 

259.9 

102.3 

11. 

.01 

1.4 

22 

14Feb75 

245.1 

73.4 

27. 

.24 

2.6 

23 

23Feb75 

235.8 

192.5 

362. 

.74 

7.0 

24 

0lMar75 

252.5 

171.5 

336. 

.45 

3.4 

25 

03Mar75 

266.9 

126.6 

141. 

0.00 

1.6 

26 

16Mar75 

223.7 

182.0 

595. 

3.20 

6.6 

27 

20Mar75 

249.0 

189.4 

140. 

.06 

2.6 

28 

26Mar75 

266.0 

104.5 

31. 

0.00 

.8 

29 

27Mar75 

229.5 

192.1 

1314. 

5.53 

10.0 

30 

04^^75 

247.6 

167.0 

370. 

.64 

5.6 

31 

10i^r75 

261.1 

168.5 

751. 

1.26 

5.0 

32 

221^15 

259.7 

135.0 

41. 

.04 

3.6 

33 

23Apr75 

230.5 

200.1 

1967. 

7.30 

9.2 

34 

29Apr75 

260.9 

145.2 

129. 

.05 

5.2 

35 

30May75 

244.1 

220.1 

817. 

2.69 

10.4 

36 

0U\jn75 

248.7 

162.4 

255. 

.98 

3.4 

37 

05Jun75 

271.7 

169.5 

131. 

.04 

3.6 

38 

llJan76 

241.5 

175.2 

1042. 

2.64 

7.4 

39 

02Mar76 

239.6 

155.7 

1127. 

2.64 

10.4 

40 

21Apr76 

233.7 

180.7 

880. 

6.25 

10.2 

41 

08Oct76 

241.4 

155.7 

1532. 

2.68 

10.0 
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water  particles.  ESI  differs  from  EWC  since  altitude  appears  as  a weighting 
factor  and  also  since  water  particles  with  diameter  less  than  50  vim  are 
excluded.  AH  is  the  total  gecrnetric  thickness  of  all  cloud  layers. 

3.5.4  Aneilysis  of  variance 

Analysis  of  variance  is  a statistical  technique  that  examines  the 
relationship  between  a dependent  variable  (predictand  or  criterion)  and 
independent  variables  (predictors) . If  the  points  on  a scatter  diagram 
(e.g. , Figures  3.6-3.11)  follow  a line,  that  regression  line  may  be  considered 
tfie  dynamic  mean. 

Although  linear  regression  is  generally  adequate,  it  is  a matter  of 
observation  that  sane  data  variates  are  not  so  connected.  The  investigation 
of  a precise  description  of  the  oonccmitant  variation  of  two  or  more  quantities 
is  one  of  the  problatis  of  curve  fitting,  known  as  curvilineau:  regression.  The 
objective,  then,  is  to  test  or  discover  a mathematical  law  which  relates  the 
variates. 

The  ctiallenge  of  curve  fitting  the  data  of  Table  3.8  was  to  determine 
a suitable  t luation.  Many  equation  types  were  considered  and  four  candidates 
were  selected  for  analysis  as  slx»m  in  Table  3.9. 

TABLE  3.9 

General  Equations  - Clouds  (Y)  and  SR  (X) 


Linear; 

Y = A + BX 

(6) 

Inverse: 

Y = A + B/(X  + C) 

(7) 

Parabolic : 

Y = A + BX  + CX2 

(8) 

Cubic: 

Y = A + BX  + CX^  + DX^ 

(9) 
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Fig  3.8  Scatter  diagram  relating  ESI  to  IR  temperature 
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CLOUD  THICKNESS 


IR  TEMP  KELVIN  - 41  CASES 

Fig  3.10  Scatter  diagram  relating  total  cloud  thickness  to  IR  temperature 


CLOUD  THICKNESS 


VISIBLE  NORMALIZED  - 41  CASES 


Fig  3.11  Scatter  diagram  relating  total  cloud  thickness  to  reflected  sunlight 
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Fran  the  data  set  of  Table  3.8  and  the  equations  of  Table  3.9,  a 
correlation  analysis  was  conducted  to  determine  the  relationship  between  the 
independent  variables  (IR  and  VIS)  and  the  dependent  variables  (IJWC,  ESI,  and 
AH) . If  all  the  data  points  were  to  lie  on  the  regressicn  line,  then  with  no 
variaticai  about  the  line,  the  correlation  would  be  perfect  with  a value  of  1. 
The  correlation  coefficient  values  thus  cotputed  for  the  cloixi-SR  data  ranged 
fron  0.55868  to  0.79602  and  are  presented  in  Table  3.10. 

TABI£  3.10 


Correlation  Coefficients  - Clouds  (IWC,  ESI,  AH)  and  SR  (IR,  VIS) 


LTNCAR 

INVERSE 

PARABOLIC 

CUBIC 

IR,  EWC 

.67495 

.68129 

.66998 

.70445 

VIS,  EWC 

.60431 

.63928 

.63492 

.64506 

IR,  ESI 

.73253 

.75681 

.72262 

.79000 

VIS,  ESI 

.55868 

.61726 

.60278 

.63317 

IR,  AH 

.77887 

.77853 

.77850 

.79602 

VIS,  AH 

.69518 

.71197 

.70992 

.71131 

Examination  of  the  correlation  values  of  Table  3.10  suggests  that  | 

curvilinear  regression  provides  a better  fit  than  linear  regressicai.  A test  f 

for  significance  Jcnown  as  the  variance  ratio  F or  sirnply  the  F-test  was  ^ 

conducted  to  examine  the  significance  of  d^arture  fron  linear  regression.  | 

Fran  tables  found  in  both  Snedeoor^^  and  Croxton  and  Ccwden^^,  a value  greater 
than  7.29  would  be  an  indication  of  significant  curvilinecu:  regression.  Only  > 

I I 

in  the  IR,  ESI  cubic  case  was  the  value  exceeded.  A further  examination  of  i 
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Kiq  3.12  Three  forms  of  c.irve-f itting  for  ESI  as  a function  of  IR  temperature 
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the  cubic  equation  revealed  that  the  curve  bends  downward  (Figure  3.12) , 
which  contradicts  the  physical  models  discussed  previously.  The  inverse 
curve  is  shewn  on  Figures  3.6  - 3.11.  The  F-test  values  of  significance  of 
departure  fron  linear  regression  are  given  in  Table  3.11. 


TABLE  3.11 

F-test  of  Significance  of  Departure  fron 
Linear  Regression 


INVERSE 

PARABOLIC 

CUBIC 

IR,  LWC 

.64 

-.49 

3.23 

VIS,  LWC 

2.94 

2.53 

3.49 

IR,  ESI 

2.04 

2.04 

9.48 

VIS,  ESI 

4.07 

2.72 

5.76 

IR,  AH 

o 

• 

0 

3.18 

VIS,  AH 

.53 

1.69 

2.04 

Ihe  observation  that  the  departures  fron  linear  regression  are  not,  in 
most  cases,  statistically  significant  is  unexpected  since  physical  models  of  | 

cloud  radiation  piroduce  highly  non-linear  relations  of  cloud  pre^jerties  to  VIS  ^ 

and  IR.  It  should  be  noted  that  the  measurements  of  cloud  pre^aerties  have 
substantial  uncertainties  v^ch  could  easily  be  misleading  in  a snail  data 

» 

sarple  such  as  Table  3.8.  Moreover,  the  linear  equations  are  often  not  the 

• * 

best  for  applications  since  they  underestimate  the  extreme  values  of  cloud 
properties,  which  are  usually  the  most  important  in  satellite  ajplications . 
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3.5.5  Multiple  regression 


Eicactly  the  same  principles  are  involved  in  multiple  regression  as 
siitple  regression,  but  the  principles  are  extended  to  more  than  one  independent 
variable.  With  the  use  of  both  independent  variables  IR  cmd  VIS,  an  estimate 
of  U'JC,  ESI,  and  AH  can  be  made  more  accurately  than  with  one  variable  alone. 

A set  of  normal  equations  can  be  derived  to  estimate  the  coefficients  A,  B,  and 
C in  the  equation  of  the  form 

Y = A + BX^  + CX2 

The  resultant  estimating  equations  fran  multiple  regression  based  on  the 
measurements  in  Table  3.8  are 

( 

I 

(11) 

(12)  ■ 

(13)  I 


IWC  = 3950.851  - 16.157  (IR)  + 3.884  (VIS) 
ESI  = 19.870  - 0.0785  (IR)  + 0.00904  (VIS) 
AH  = 30.597  - 0.119  (IR)  + 0.0283  (VIS) 


(10) 


A sumiary  of  multiple  correlation  and  standard  error  of  the  estimate  is  shewn 

I 

in  Table  3 . 2 . 

I 

t 

Table  3.12 


Multiple  Correlation  Coefficient  and  Standard  Error  of 

the  Estimate  for  41  Cases  of  Cloud  Preperties  ^ 


Multiple 

Correlation  Coef. 

Standard 

Error  of  Estimate 

IWC 

.72071 

369.208 

ESI 

.74842 

1.354 

AH 

.83073 

1.889 
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VISIBLE  NORMALIZED  - 41  CASES 


IR  TEMP  KELVIN  - LWC  CONTOURS 

Fiq  3.13  Scatter  diaqram  relatinq  LWC  to  both  IR  temperature  and  reflected  sunlicht. 
Solutions  of  the  linear  Equation  (11)  are  ai&o  plotted. 


VISIBLE  NORMALIZED  - 41  CASES 


VISIBLE  NORMALIZED  - 4 I CASES 


300 


Fig  3.15  Scatter  diagram  relating  total  cloud  thickness  to  both  IR  temperature  and 
reflected  sunlight.  150 


The  regression  plane  defined  by  the  equation  is  designed  to  pass  among 
the  points  so  that  the  sun  of  squares  of  vertical  distances  fron  plane  to 
points  is  a minimum.  However,  the  statistics  of  multiple  regression  are  subject 
to  sanpling  variation  and  a final  test  of  significance  was  conducted.  An 
F-test  similar  to  that  conducted  in  Section  3.5.4  above  proved  highly  signifi- 
cant. This  means  that  IR  and  VIS  as  joint  predictors  of  LWC,  ESI,  and  AH  are 
almost  certainly  valid  and  offer  improved  estimates  over  IR  or  VIS  alone. 
Solutions  of  equations  (11) , (12) , and  (13)  and  shewn  as  straight  lines  in 
Figiures  3.13,  3.14,  and  3.15.  In  multiple  regression,  as  in  single  regression, 
the  IR  is  a better  predictor  than  VIS.  This  conclusion  could  be  expressed  in 
terms  of  statistical  tests,  but  it  is  readily  apparent  by  oenparing  the 
regression  slope  coefficients  of  IR  and  VIS  in  equations  (11) , (12) , and  (13) . 

Multiple  regression  was  also  done  for  inverse  forms  of  the  VIS  and  IR, 
as  independent  variables.  These  nonlinear  forms  generated  slightly  higher 
multiple  correlation  coefficients  and  slightly  lower  standard  errors  of 
estimate  for  the  cloud  properties  LWC  and  ESI.  There  was  less  inprovanent  in 
the  case  of  AH. 

3.5.6  Conclusions 

Clouds  vdrLch  are  brightest  in  the  VIS  and  coldest  in  the  IR  have  the 
largest  LWC,  ESI,  and  AH.  IR  terperatures  are  better  predictors  than  VIS 
f by  a small  amount.  The  equations  to  predict  cloud  properties  from  satellite 

* IR  and  VIS  have  statistical  significance  according  to  established  F tests. 

Hie  siitple  form  of  the  prediction  equations  lends  itself  readily  to  machine 
cenputation  of  cloud  property  estimates.  Nonlinear  forms  of  prediction 
equations  offer  slightly  better  predictions  than  the  linear  forms. 
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3.6  Ctojective  Prediction  of  Cloud  Cover  and  Precipitation 
3.6.1  Introduction 

Hie  Air  Weather  Service  is  making  significant  changes  in  its  c^jeraticxial 
format  to  meet  demands  for  increasing  services  with  fewer  pecple.  Anong  the 
goals  being  sought  by  the  AWS  are  automation  of  observations,  centralization  of 
product  preparation,  and  tailoring  selected  products  to  be  used  by  the  c^jerator 
with  a minimum  of  interpretation  by  meteorologists.  In  many  cases,  the  technol- 
ogy needed  to  achieve  these  goals  is  readily  available  and  needs  only  be 
modified  for  application  to  a specific  prctolem. 

At  the  present  time  sun-synchronous  and  geosynchronous  satellites  cover 
the  globe  many  times  per  day  providing  large  quantities  of  detailed,  continuous 
information  not  readily  available  frcm  other  sources,  ffcwever,  two  special 
prc±)lems  are  posed  in  the  ^pdication  of  satellite  data  to  aubcmated  processing. 
Ihe  first  problan  is  handling  the  huge  volune  of  data  generated  by  the  imagery 
sensors.  An  investigation  of  one  aspect  of  this  prctolem  is  reported  in 
Section  3.4  of  this  report.  A second  problem  is  translating  the  necessarily 
subjective  satellite  photoanalysis  into  an  automated  ctojective  form  vtoich  can 
be  processed  by  ocitputers.  Ihe  ability  to  achieve  a primitive  or  partial 
solution  to  this  second  problem  is  inplicit  in  the  work  reported  herein. 

3.6.2  Objective 

The  space  and  time  scales  of  geosynchronous  satellite  imagery  are 
consistent  with  the  needs  of  short-range,  mesoscale  forecasting.  Ihe  AWS  has 
a requirement  for  continuously  updated,  autcmated,  short  range  forecasts  of 
weather  events  that  are  important  to  operational  missions.  The  work  in  this 
investigation  is  one  phase  of  a larger  study  concerned  with  the  developnent 
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1.16  Area  of  satellite  imagery  data.  Surface  reports  saved  for  stations  with  bold 
markings . 


1 

I 

I 

I 

I 


of  techniques  to  utilize  real-time  satellite  imagery  data  in  the  field  of 
objective  weather  forecasting.  The  intent  is  eventually  to  use  satellite  data 
in  conjunction  with  conventional  data  for  application  to  terminal  forecasting; 
and  by  itself  for  application  to  target  forecasting  or  for  areas  that  are 
without  weather  stations. 

The  ultimate  goal  of  this  study  is  a system  vAiich  extracts  the  appropri- 
ate data  frciti  the  satellite  transmission,  processes  it  thru  a coiputerized 
statistical  package,  and  generates  a forecast  with  the  probability  of  a criticcil 
event  and  its  duration  to  the  operational  decision  maker. 

The  irmediate  objectives  are  a little  more  straightforward.  For  the 
first  st^  we  wish  to  establish  the  reliability  of  the  data,  and  then  investi- 
gate the  variation  in  surface  brightness  as  a function  of  location,  season, 
and  time  of  day.  This  is  important  because  it  is  against  these  brightness 
values  that  the  approach,  retreat,  and  development  of  clouds  will  be  measured. 

3.6.3  Data  sample 

The  geosynchronous  meteorological  satellite  data  used  are  hourly  tapes 
of  1-mile  resolution  visible  and  4-n\ile  resolution  infrared  data  vdiich  cover 
the  northeastern  area  of  the  U.S.  shown  in  Figure  3.16. 

The  visible  data  represents  scene  brightness:  low  numbers  indicate  water 
and  ground  surfaces;  middle  range  numbers,  low  or  thin  clouds,  and  srrw  covered 
land  and  deserts;  and  large  numbers,  high  or  thick  clouds  and  ice  cap®.  The 
infrared  data  represents  terperatures  at  the  tq?  of  clouds  or,  in  the  absence 
of  clouds,  the  earth's  surface:  very  cold  temperatures  r^resenting  high  clouds; 
warmer  temperatures  representing  lower  clouds  or  the  earth's  surface  in  cold 
climates;  and  the  vrarmest  terperatures  representing  deserts  and  trcpical  oceans. 


I 
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In  order  to  identify  cloud  amounts  and  predict  their  motion,  the 
brightness  signature  of  the  earth's  surface  must  first  be  determined.  This 
in  turn  varies  with  the  land/water  configuration,  the  season  of  the  year, 
i.e.,  color  of  the  ground,  and  the  time  of  day. 

IVo  hundred  fifteen  essentially  clear  observations  frcm  March  to 
December  1977  were  chosen  to  determine  the  ground  signature  for  a given  day, 
time  and  location.  The  locations  chosen  for  this  study  are  Bedford,  Boston, 
Albany,  Pittsburg,  and  EXilles  International  Airport  in  Virginia.  Plans  called 
for  collecting  data  throughout  the  working  day.  Preerption  of  the  system  by 
higher  priority  projects  prevented  the  collection  of  afternoon  data  during  May, 
June,  and  July. 

With  the  given  stations  as  the  center  points,  visible  data  to  cover  120 
rows  by  120  elonents  and  infrared  data  to  cover  30  rcws  by  132  elanents  were 
located  and  extracted  frcm  the  archived  data.  Ranking  these  ntmerical  sets  of 
visible  (VS)  and  infrared  (IR)  data  into  24  classes,  and  assigning  than  the 
letters  A - W,vrtien  printed  out,  gave  a dramatic  visual  map  of  the  chosen  areas 
centered  on  the  above  locations.  Land  and  water  configurations  could  be  easily 
seen  and  idt.itified  especially  along  the  New  England  coastline.  These  printouts 
were  initially  used  to  verify  the  locations  of  landmarks  in  the  data  array. 

3.6.4  Data  corrections 

The  McIDAS  system  is  used  interactively  to  obtain  a correction  factor  to 
fine-tune  the  navigation  to  a 2-mile  or  better  accuracy  necessary  for  the  fine 
time  and  space  scales  required  for  mesoscale  forecasting.  This  correction 
factor  is  determined  by  the  difference  between  the  line  and  elonent  specified 
as  the  geographic  coordinates  of  a landmark  after  the  initial  navigation,  and 
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those  obtained  vdien  the  cursor  is  placed  on  that  landmark  on  the  screen.  A 
2-mile  accuracy  results  vhen  3 good  landmarks  are  visible,  and  3-4  mile 
accuracy  for  2 landmarks.  This  correction  factor  is  valid  for  all  points  shown 
in  Figure  3.16. 

Examination  of  the  averages  of  brightness  for  entire  lines  of  data 
indicated  systematic  variations  frcm  line  to  line.  These  variations  were 
indicative  of  slight  differences  in  the  sensitivity  of  the  sensors.  Although 
the  variation  amovinted  to  only  about  1%,  the  brightness  values  were  adjusted. 

The  average  of  all  picture  elanents  (pixels)  sensed  by  each  of  the  4 scanners  j 

individually  were  made  to  agree  by  adding  or  subtracting  appropriate  constants  ; 

fron  each  brightness  value.  ! 

Randan  RF  noise  was  also  taken  into  consideration  and,  if  detected,  an 

i 

adjustment  was  made  in  the  data.  Noise  was  identified  v^hen  the  average  bright- 
ness for  a line  departed  more  than  approximately  13%  fron  the  average  for  the  | 

previous  line.  When  this  occurred,  each  pixel  in  the  line  was  replaced  by  the  ! 

average  value  of  the  pixel  in  the  non-noise  line  above  and  the  pixel  in  the 

next  non-noise  line  below.  ( 

3.6.5  Data  sampling 

The  spatial  stability  of  the  brightness  data  was  examined  by  selecting  [ 

' ^ 

an  area  of  27  x 27  pixels  surrounding  the  selected  site.  This  27  x 27  (729)  [. 

array  was  subdivided  into  nine  9x9  boxes,  Figvire  3.17.  Each  box  was  further 

I. 

subdivided  into  7x7,  5x5,  and  3x3  boxes.  The  total  27  by  27  box  represents  approx-  ^ 

imately  25  by  35  mile  area  vAiich  is  close  in  size  to  the  3-D  Nephanalysis  lanit.  ^ 

For  each  9,  25,  49,  81  and  total  729-pixel  box,  the  average,  standard 
deviation,  maximum,  miniimin,  and  range  of  brightness  were  cerputed.  Figure  3.18 
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CENTER 


BLOCK 

ONE 

TWO 

THREE 

FOUR 

AV6 

21.528 

20.990 

21.235 

21.228 

SD 

3.495 

3.  158 

3.524 

3.779 

MAX 

31 

31 

36 

38 

MIN 

19 

19 

19 

19 

RANGE 

12 

12 

17 

19 

Fiq  3.18  Sample  of  the  format  in  which  9x9  pixel-box  data,  stored  on  tape,  can  be 

printed.  Each  printout-sheet  contains  this  format  for  the  central  and  8 surrounding 
boxes  arranged  as  in  Figure  3.17. 
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is  a sample  of  the  output  data  for  one  9x9  pixel  box  \A^ch  is  stored  on  tape 
for  each  hour  used  in  the  study. 

A ccrtparison  of  the  average,  standard  deviation,  maximum,  minimun  and 
range  of  brightness  for  the  9,  25,  81,  ard  729  pixel  boxes  centered  on  the 
selected  site  on  clear  days  is  shown  in  Table  3.13.  The  chart  indicated  that 
the  statistics  are  relatively  stable  vrfien  corrputed  fran  9 to  729  pixels  and 
that  9 pixels  can  provide  sufficient  accuracy  for  most  uses. 

Table  3.14  shows  the  brightness  statistics  for  nine  9x9  pixel  areas 
arranged  as  in  Figure  3.17  for  four  clear  observations  at  16Z,  and  an  observa- 

I 

tion  with  scattered  clouds,  shown  in  parentheses,  on  one  of  the  same  days  at  | 

14Z.  The  statistics  of  the  clear  observations  shew  ramrkable  similarity  frem 

box  to  box.  When  scattered  clouds  viere  reported,  the  standard  deviation  was 

as  much  as  5 times  larger  and  the  maximum  values  were  higher  because  of  the 

clouds.  The  minimuns  were  lower  because  the  sun  was  lower  in  the  sky.  j 

This  stability  of  the  ground  brightness  values  over  a range  of  sample  | 

sizes  and  spatial  distributions  bespeaks  an  instrunental  stability  that  should 
make  possible  the  detection  of  subtle  changes  in  clouds  and  cloud  configure-  I 

tions. 

3.6.6  Brightness  analysis 

/\n  analysis  was  done  of  the  average  brightness  for  clear  days  at  each 
of  the  4 key  .'•tations,  Bedford,  Albany,  Pittsburg,  cind  Dulles.  An  exanple 
for  July  is  shown  in  Figixe  3.19.  Low  scattered  clouds,  newly  forming,  are  ^ 

indicated  when  the  standard  deviation  is  less  than  10;  more  advanced  scattered 
cloud  development,  denser  scattered  clouds,  or  bigger  patches  of  cloudiness  f 

are  indicated  when  the  standard  deviation  is  greater  than  10.  As  can  be  seen. 
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TABLE  3.13  Brightness  Statistics  at  Bedford  MA  at  1600Z  as  Function  of  Area  on  April  29,  May  3 
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AVERAGE  BRIGHTNESS 


SOLAR  ELEVATION  ANGLE 

Fig  3.19  Scatter  diagram  of  average  ground  brightness  of  a 9x9  pixel  boxes  centered 
at  the  4 test  stations  during  clear  and  scatter  conditions. 
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when  the  standard  deviation  is  amall  the  points  fall  along  the  same  line  as  the 
clear  cases.  Some  of  the  widely  scattered  cases  are  included  in  the  statistics 
to  enlarge  the  sample  for  this  initial  study. 

High,  thin  scattered  clouds  have  little  effect  on  the  average  bright- 
ness, and  clear  sky  data  may  also  include  sane  of  these  cases.  In  contrast, 
haze  did  have  a noticeable  effect  characterized  by  a brighter  shade  and  a very 
small  standard  deviation.  Haze  cases  were  excluded  and  will  be  handled  later. 

In  order  to  determine  whether  the  minimum  or  average  brightness  is  the 
better  indicator  of  clear  skys,  both  the  nvinimum  and  average  brightness  for  ; 

I 

729  pixels  were  correlated  with  solar  elevation  angle  for  selected  days  and 
hours  at  different  locations. 

The  y intercepts  and  tlie  slopes  of  the  regression  equations  were 
CLTputed,  as  was  the  standard  error  of  estimate  Sy.  Figures  3.20  thru  3.24 
show  that  the  resulting  correlation  coefficients,  standard  errors  of  estimate,  I 

and  the  slopes  are  almost  identical  for  both  the  minimum  and  average  brightness.  I 

3.6.7  Conclusions 

Thar  does  indeed  exist  a stable  relationship  between  background  | 

I 

brightness  and  solar  elevation  angle  from  day  to  day  and  place  to  place.  f 

I 

Variations  result  from  seasonal  changes  in  the  characteristics  of  the  surface  | 

L 

texture  ard  coloration.  In  this  sample,  brightness  did  not  vary  greatly  among 
the  4 key  locations,  but  a larger  sample  of  data  from  this  same  region  might 

i 

change  this  conclusion.  Certainly  it  will  change  as  samples  are  analyzed  from 
other  than  the  rortheastem  United  States. 

A sutimary  of  average  brightness  over  729  pixels  vs.  solar  angle  is 
given  in  Figure  3.25.  Almost  all  correlation  coefficients  were  above  0.90  for 
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SOLAR  ELEVATION  ANGLE 

3.20  Graph  of  the  regression  equation  for  average  (B)and  minimum  ground 

brightness  vs.  solar  elevation  angle  on  clear  days  for  stations  indicated  on  graph 
for  March  1977.  r is  the  correlation  coefficient;  Sy , the  standard  error  of 
estimate-  and  n,  the  number  of  cases. 
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BRIGHTNESS 


Fiq  3.21  Same  as  Figure  3.20  except  for  Apri'' 


BRIGHTNESS 
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Fig  3.22  Same  as  Figure  3.20  except  for  May,  June  and  July. 
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Fig  3.24  Same  as  Figure  3.20  except  for  December 
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these  cxmbined  locations.  In  all  but  one  case,  Decenter,  the  average  bright- 
ness have  higher  correlation  coefficients  and  lower  standard  errors  of  estiinate 
than  the  miniimjn  brightness. 

The  need  new  is  to  enlarge  the  data  base  to  collect  a significant 
sanple  for  each  station  for  each  month.  This  larger  sanple  vd.ll  refine  sane 
of  the  conclusions,  rather  than  significantly  change  any  of  the  results. 

Ihe  next  phase  of  planned  work  is  an  hour-by-hour  daily  analysis  of 
cloud  cases  for  (a)  distance,  (b)  direction,  and  (c)  speed  of  movanent  of 
brightness  areas  outside  the  729  pixel  box  into  the  central  location  to  esti- 
mate the  predictability  of  the  advection  and  developnent  of  objectively  defined  ' 

clouds. 

3.7  Cloud  Pattern  Motion  Vectors  | 

I 

3.7.1  Introduction 

The  purpose  of  this  effort  is  to  formulate  a working  corputer  model  for 

I 

the  determination  of  cloud  pattern  motion  vectors  based  on  an  approach  developed 
primarily  by  the  Stanford  Research  Institute  under  contract  to  the  U.S.  Navy. 

A cloud  pattern  is  represented  by  a two-dimensional  array  of  brightness  values.  t 

The  version  of  the  oaiputer  program  assembled  under  this  effort  for  calculating  1 

motion  vectors  will  be  used  for  stvdying  effects  of  the  various  parameters 
involved  and,  subsequently,  for  making  cemparisons  with  results  from  other  ' 

approaches.  The  ultimate  goal  at  AFGL  is  to  formulate  a fully  autonated  tech-  i 

nique  for,  obtaining  cloud  pattern  motion  vectors  for  use  in  shortrange  sensible  ’ * 

k' 

weather  prediction  stixiies. 

The  requirement  that  the  technique  be  automatic,  being  aocemplished  i 

■ t . 

fran  beginning  to  end  by  oemputer,  greatly  increases  the  difficulty  of  such  a 
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developnent  as  catpared  to  those  requiranents  vdiich  can  settle  for  a semi- 
autonatic  approach.  However,  there  are  a few  candidates  for  the  automatic 
approach  to  notion  vectors  Viiiich  show  premise.  The  one  formulated  here  is 
patterned  after  the  one  developed  at  the  Stanford  Research  Institute  (SRI)  and 

1 p 

described  by  Endlich  et  al. 

Much  of  the  wark  in  recent  years  on  calculating  cloud  notion  vectors 
has  either  explicitly  or  inplicitly  had  as  a goal  the  estimation  of  wind 
velocities.  This  is  not  our  goal  either  explicitly  or  inplicitly.  Cloud 
pattern  notion  vectors  have  a raison  d'etre  apart  from  their  inplication  of 
winds;  namely  as  an  aid  in  generating  cloud  and  other  weather  element  forecasts. 

The  ward  pattern  is  used  to  anphasize  the  fact  that  the  goal  is  that  of 

I 

tracking  arrays  of  brightness  instead  of  individual  elements  of  brightness.  j 

The  object  pursued  here  is  scroevSiat  different  frem  that  vAuch  gave  rise  ' 

to  the  SRI  model.  A single  velocity  vector  vdiich  is  representative  of  an  area 
of  the  order  of  90x90  n.mi.  is  sought.  SRI  seeks  synoptic  or  sub-synoptic  ^ 

motion  vectors  over  large  areas  (vAiere  clouds  permit  of  course) . The  Navy 
(NEPRF,  i’-onterey  CA)  has  recently  begun  an  operational  use  of  an  interactive 
version  of  a lore  recent  SRI  model.  Although  this  more  recent  model  xandoubtedly 
includes  experience  gained  over  tlie  past  few  years,  the  principal  change  has 
been  teward  a sinplification  of  the  clustering  process  for  the  sake  of  caiputer  , 

r 

time  savings  rather  than  for  improved  clustering.  Clustering  techniques  are  ; 

usually  developed  and  applied  to  assist  in  the  analysis  of  large  amounts  of  i 

data.  Iheir  developnent  and  use  has  been  very  rapid  in  recent  years.  Briefly,  • 

'.I 

clustering  techniques  are  oonputational  means  of  grouping  data.  They  are 
usually  thought  of  as  techniques  for  finding  natural  groupings  of  data  since 
most  of  the  research  activity  and  applications  occur  in  this  area. 
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Our  approach  to  motion  vector  analysis  uses  briqhtness  centers  to  represent  ' -I 

cloud  patterns.  Given  two  consecutive  cloud  pictures,  brightness  centers  are  j 

I 

obtained  for  both,  and  then  the  two  sets  of  centers  are  matched.  Fran  this  | 

19 

matching  an  overall  vector  is  obtained.  Ball  ^ al  are  accredited  with  j 

develcping  this  approach  around  their  ISODATA  clustering  algorithm.  No  effort 
has  been  made  here  to  inprove  on  the  technique  of  Endlich  ^ ^ (1971) ; an 
effort  has  been  made  to  interpret  it  as  accurately  as  possible  fron  their  report. 

Differences  most  likely  do  exist  since  the  write-up  does  not  give  enough 

details  to  be  certain  vtot  was  done.  This  applies  especially  to  the  routine 

for  finding  brightness  centers.  Strengths  and  wea3aiesses  of  the  steps  taken 

will  show  up  in  experiments  with  the  model  and  ccrtparisons  of  it  with  other 

techniques . | 

i 

3.7.2  Three  significant  steps  to  motion  vector  analysis 

There  are  three  significant  steps  to  motion  vector  analysis  in  this  i 

approach:  (a)  Simplification  of  cloixl  images  by  truncating  at  seme  brightness  ! 

level  or  by  some  other  process  such  as  spacial  averaging;  (b)  application  of 
a clustering  procedure  to  obtain  representative  brightness  centers;  (c)  Pairing  , 

off  of  members  of  two  sets  of  brightness  centers  to  obtain  motion  vectors.  A , 

I 

discussion  of  each  of  these  is  in  order  before  laying  out  the  details  of  the  * 

caiputer  version  assembled  for  AFGL.  The  flexibility  of  this  approach  is,  as 

i 

might  be  expected  fron  the  above  steps,  considerable.  ! 

(a)  Simplification  of  Cloud  Images  • i 

Simplification  of  the  images  is  necessary  because  of  limitations  ^ 

> 

inposed  by  oorputer  storage  and  running  time.  Even  with  simplification  by 
truncating  to  a one-bit  brightness  level  and  taking,  for  exanple,  every  other 
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elanent  of  every  other  row,  the  problem  of  obtaining  good  sets  of  brightness 
centers  by  clustering  is  no  minor  process  in  terms  of  corpuber  time.  Siirpli- 
fication  is  a necessity.  There  are  many  ways  this  can  be  done  other  than 
truncation.  The  main  idea  here  is  to  obtain  a simplification  that  best  repre- 
sents those  features  that  are  of  interest  for  tracking.  This  step  has  received 
little  cament  in  the  literature.  Hew  best  to  siiiplify?  This  is  a good 
questiort  but,  so  far,  there  have  been  no  good  answers.  The  answer  could  be 
inportant  in  any  efforts  to  inprove  the  model  such  as  separating  cloud  motions 
at  different  levels  or  obtaining  the  motion  most  significant  in  terms  of 
station  forecasts.  In  addition  to  this,  there  are  many  more  unknowns  vhich 
require  extensive  experiments  in  order  to  understand  them  well.  Different 
channels  could  be  used  separately  or  oorposited  in  sane  way.  The  possibilities 
for  this  step  seem  unlimited.  The  first  step,  to  save  time  and  energy,  for  ' 

an  indepth  exploration  would  naturally  be  time  spent  at  an  interactive  graphics  I 

console  to  narrow  dewn  the  field  of  possibilities.  j 

(b)  Clustering  Brightness  Areas  i 

The  clustering  of  brightness  areas  and  use  of  their  centers  to 
represent  cr  jd  patterns  is  the  core  of  this  approach  and  is  at  the  same  time 
its  vreakest  link.  Simplifications  in  the  clustering  process  are  required  at 

1 

the  expense  of  cotpleteness . The  guiding  idea  here  is  to  take  the  sinplified  , 

cloud  images  and  obtain  locations  of  major  brightness  centers.  Then  two  ^ 

consecutive  images  would  be  ejq^ected  to  resemble  each  other  except  for  the  1 

shift  resiUting  frem  motion  of  the  cloud  pattern.  The  procedure  must  contend  ^ 

with  cotplications  such  as  (a)  patterns  moving  partially  off  (or  over)  the 
area  of  interest,  (b)  changes  in  brightness  areas  or  number  of  brightness 
centers  due  to  growth  or  dissipation,  (c)  brightness  centers  splitting  or 


173 


conbining,  and  (d)  coping  with  the  flimsy,  irregularity  of  vi^t  we  refer  to 
as  "clouds"  and  "cloud  patterns". 

(c)  Pairing  the  TW3  Sets  of  Brightness  Centers 

Once  brightness  centers  are  obtained  the  remaining  problan  is  to 
pair,  whejre  possible,  centers  on  the  first  image  with  centers  on  the  second 
image  in  such  a way  that  they  would  be  judged  the  same  points  in  the  pattern  at 
different  times.  The  difficulty  here  lies  in  designing  the  model  to  cover  a 
very  wide  range  of  cases.  First  of  all,  the  number  of  points  in  the  two  sets 
need  not  be  the  same  for  the  reasons  mentioned  in  the  previous  paragraph. 

Inportant  factors  in  this  analysis  is  the  distance  between  brightness  centers, 
the  nunber  of  centers,  and  the  area  associated  with  the  centers.  These  factors 
are  obviously  interrelated.  Another  important  consideration  is  the  interval 
of  time  between  the  two  images.  The  distance  between  centers  relative  to 
motion  distance  is  an  important  factor  in  the  process  of  pairing.  Best  results 

i 

are  expected  vdien  center  separation  distances  are  large  relative  to  the  motion  i 

j 

distance  and  vice  versa.  : 

Each  of  the  three  important  steps  in  the  analysis,  simplification, 
clustering  and  pairing,  can  be  ipigraded  without  serious  impact  on  the  other  two.  ! 

Which  peui:s  are  in  need  of  improvement,  if  any,  can  be  determined  from  rather 

simple  comparisons  with  more  reliable  semi-automatic  methods.  { 

; ' 

L 

3.7.3  Computational  details  of  motion  vector  model 

I 

Details  of  the  coiputational  imodel  will  new  be  given  in  sufficient  • t 

detail  to  permit  dvplication.  As  coipared  to  the  actiial  implementation  by  SRI,  ‘ 

the  model  develcped  here  is  simple  and  straightforward.  Ihis  simplicity  of 
structure  permits  modifications  and  improvements  vrfiere  necessary  without 
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getting  into  ccn^lications  that  rcxjuiie  an  unreasonable  amount  of  experimenta- 
tion. It  seems  to  us  there  is  little  reason  to  enploy  a clustering  routine 
that  goes  to  as  great  a depth  of  precision  as  does  the  SRI  ISODATA  routine 
v^rtien  the  problans  of  accounting  for  cloud  pattern  changes  other  than  transla- 
tion are  a major  source  of  error.  Recent  developments  at  SRI  of  their  model 
indicates  an  agreement  on  this  point.  In  any  event,  the  technique  developexi 
here  will  yield  improvements  in  v\hatever  ccmponents  that  might  be  needed.  It 
will  also  yield  readily  to  adjustments  to  different  images,  image  sizes,  and 
pattern  scales. 

a.  Preparation  of  Arrays 

The  first  step  was  to  develop  a procedure  that  would  permit  computer 

runs  of  the  model  on  many  different  kinds  of  cloud  patterns.  A program  was  I 

written  (COZHA)  vhich  takes  a prescribed  array  from  a data  tape  of  the 

Geostationary  Operational  Environmental  Satellite  (GOES)  and  places  it  into  a 

permanent  file  at  the  AFGL  Ccmputer  Center.  These  arrays  are  kept  in  the  same 

resolution  as  is  recorded  on  the  satellite  data  tapes.  For  the  visible  data 

■^hannel,  which  is  of  current  interest,  there  is  afproximately  one  picture 

elanent  per  lile.  These  images  consist  of  every  other  elanent  and  every  other 

line  of  the  initial  one-half  mile  resolution  GOES  data.  Simplifications  in 

these  data  are  needed  that  ^lre  most  advantageous  to  the  model. 

The  array  size  chosen  for  current  studies  is  90x90  n.mi.  The  selection 
t i 

of  this  size  area  r^resents  an  effort  to  obtain  a balarxre  between  (a)  snail  ' 

. ' 

. enough  to  be  representative  of  snail  scale  changes,  (b)  large  enough  to  include  | 

enough  cloud  pattern  to  obtain  a representative  reading  of  motion. 

, The  oorputer  program  COZHB  takes  arrays  at  two  different  times  and 

obtains  frcm  them  a motion  vector. 
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b.  Sinplification  by  Truncation 

There  are  many  ways  to  simplify  an  image  for  further  processing  and 
little  is  known  so  far  about  advantages  and  disadvantages  of  them.  Simplifi- 
cation is  needed  here  to  bring  the  data  within  a size  that  is  reasonable  in 
terms  of  corputer  storage  capacities  and  times  required  for  execution. 

Fran  the  90x90  picture  elanent  array  every  third  column  and  every  third 
rcw  is  taken.  This  results  in  a 30x30  array.  The  hundred  brightest  points 
are  selected  frcm  this  array  for  further  use.  An  elenent  is  either  in  the 
hundred  brightest  or  it  is  not.  Actually,  since  there  may  be  more  than  one 
point  that  ranks  c.ie-hundred , there  may  be  a few  points  more  than  a hundred  in 
this  brightest  category. 

c.  Calculation  of  Brightness  Centers 

Brightness  centers  are  obtained  fran  the  arrays  of  truncated  bright- 
ness by  a clustering  procedure  using  Euclidean  distance  as  a measure  of 

proximity  of  points.  Passes  are  made  through  the  data  points  to  select  bright-  j 

ness  center  locations.  Each  pass  operates  on  the  data  in  a similar  way.  The 
difference  or  improvanent  fron  one  pass  to  another  lies  in  the  fact  that  sane 

t 

center  locations  (clusters)  which  do  not  have  enough  points  in  them  are  dropped. 

This  provides  a means  for  filtering  out  anall  areas  of  brightness.  We  have 

found  that  two  passes  of  the  routine  for  finding  centers,  are  sufficient  to  i 

1. 

ctotain  stable  results.  A minimum  requirement  of  six  points  per  brightness 
center  gives  results  that  smooth  out  small  irregularities. 

« 

The  procedure  for  locating  centers  rests  on  the  computational 
operations  of  a pass  through  the  list  of  x,  y ooiponents  of  the  data  set 
composed  of  those  hundred -plus  highest  brightness  values. 

A description  of  such  a pass  will  be  given.  The  center  of  gravity 
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of  all  the  points  is  obtained.  This  represents  the  location  of  the  first 
hrightJTess  center.  The  points  assigned  to  this  center  are  those  which  lie 
within  half  the  average  distance  of  all  points  from  the  center  of  gravity. 

This  fraction  (one-half)  is  often  referred  to  as  the  sphere  factor.  The  value 
of  one-half  was  found  to  work  well  fran  experience  by  SRI.  The  gtiall  anount 
of  experimenting  we  have  done  indicates  that  it  is  a good  value  for  this  model. 
The  next  center  is  the  first  x,  y encountered  (reading  left  to  right,  top  to 
bottcm)  that  does  not  lie  in  the  first  cluster.  The  points  assigned  to  this 
center  are  those  which  fall  within  the  average  distance  but  are  not  already 
assigned  to  the  first  cluster.  Subsequent  centers  are  found  in  the  same  way 
until  all  the  points  are  exhausted.  Clusters  having  fewer  than  six  points  in 
them  are  discarded.  New  centers  of  the  clusters,  excluding  the  first,  of 
course,  are  defined  as  the  average  of  the  coordinates  of  the  points  within  the 
clusters.  Splitting  of  large  clusters  or  ccmbining  of  small  ones  occurs  on 
the  second  pass  only  as  a result  of  changes  due  to  the  dropping  of  the  very 
anall  clusters.  It  has  been  found  that  a second  pass  through  the  data  is  all 
that  is  required  to  obtain  good  results.  Additional  passes  through  tlie  data  set 
do  not  chant  results  enough  to  make  them  worthwhile, 
d.  Matching  of  Brightness  Centers 

The  cluster  centers  for  the  two  image  times  require  a means  of 
pairing  off,  vdiere  possible,  to  obtain  measures  of  displacements.  The  number 
of  centers  for  the  two  image  times  may  be  different  and  there  may  be  centers 
without  matches,  i.e.,  those  that  have  moved  on  or  off  the  area  of  the  image. 
The  approach  to  matching  centers  followed  here  relies  heavily  on  the  center  of 
gravity  of  the  brightness  field. 

The  motion  of  the  center  of  gravity  is  taken  as  an  initial 
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displacarient . This  is  always  one  ccrrponent  of  the  final  motion  vector.  CXiier 
corponents  are  obtained  when  pairs  of  points  lie  within  20  grid  intervals  of 
each  other  relative  to  this  displacanent . This  method  follows  SRI ' s approach 
and  rests  on  their  expcirimentation.  It  is  expected  that  this  parameter  (20 
grid  intervals)  will  require  sane  refinanent  by  making  it  a function  of  both 
grid  interval  size  and  time  between  images.  It  is  also  expected  that  the 
center  of  gravity  carries  too  much  weight  and  that  modification  to  lessen  its 
influence  would  improve  the  results, 
e . Exanple 

Figxires  3.26  and  3.27  show  arrays  of  points  (after  two  passes  of  the 
subroutine  for  obtaining  cluster  centers)  obtained  frcm  two  GOES  images  taken 
one  hour  apart.  The  grid  size  is  approximately  6 n.mi.  in  both  the  x and  y 
direction.  The  coordinates,  row  and  colunn  numbers,  appear  on  the  left  and  at 
the  tcp,  respectively.  The  difference  in  the  column  numbers  resulted  fron  a 
relative  gridding  error  between  the  two  imaqes.  This  error  was  found  by  making 
corparisons  of  geographical  features.  The  centers  of  briqhtness  are  marked  by  a 
C and  a subscript  vdiich  is  a cluster  number  for  referencing  purposes.  Table 
3.15  shows  the  distance  of  sepjaration  of  the  two  sets  of  centers  relative  to 
the  motion  of  the  center  of  gravity.  Distances  are  in  grid  intervals.  Pairs 
selected  are  shown  on  the  right  side  of  the  table.  Having  obtained  these  pairs, 
an  overall  motion  vector  is  obtained  from  the  vectors  defined  by  each  pxair  of 
» pxjints.  An  average  of  the  x and  y cotponents  of  these  vectors  produce  con- 

ponents  of  an  average  vector.  A velocity  of  7.4  n.mi. /hr  with  motion  fron  315 
degrees  was  obtained  for  this  case. 
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3.7.4  Sumnary 


'I'ho  model  for  tho  autanatic  calculation  of  motion  vectors  described 
here  was  based  on  £in  early  version  of  the  SRI  model.  Ihere  are  several  param- 
eters which  are  poorly  understood  and  require  experimental  results  for  their 
clarification.  Tho  more  inportant  ones  have  been  mentioned  here.  The  greatest 
asset  of  this  approach  is  its  flexibility.  It  should  be  rather  easy  to  pick 
out  those  areas  vhich  require  more  sophisticated  ideas  to  inprove  the  results. 
Development  of  such  ideas  are  not  expected  to  be  terribly  difficult. 
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